NOTICE 


THIS DOCUMENT HAS BEEN REPRODUCED FROM 
MICROFICHE. ALTHOUGH IT IS RECOGNIZED THAT 
CERTAIN PORTIONS ARE ILLEGIBLE, IT IS BEING RELEASED 
IN THE INTEREST OF MAKING AVAILABLE AS MUCH 
INFORMATION AS POSSIBLE 



M£0-3<f3- 

(NAS A-CR-163208) ANALYSIS AND EVALUATION Ii N80-25782 

THE PRODUCTION PROCESS AMD EQUIPMENT AREA OF 
THE LOB-COST SOLAR AHRAY PROJECT Annual 

Report, Oct. 1978 - Oct. 1979 (Pennsylvania Unclas 

Univ.) 203 P HC A1Q/MF A0 1 CSCL 10A G3/44 22379 



Contract JPL-954796 
Annual Report Oct. 1973 to Oct. 1979 
(DRD Line Item 7) 


October 1979 



UNIVERSITY of PENNSYLVANIA 

The Moore School of Electrical Engineering 

PHILADELPHIA, PENNSYLVANIA 19104 


The JPL Low-Cost Silicon Solar Array Project is sponsored by the 
U.S. Department of Energy and forms part of the Solar Photovoltaic 
Conversion Program to initiate a major effort toward the develop- 
ment of low-cost solar arrays. This work was performed for the 
Jet Propulsion Laboratory, California Institute of Technology by 
agreement between NASA and DoE 



DoE/JPL-954796-79/7 
Distribution Category UC-63 


ANALYSIS AND EVALUATION 

IN THE PRODUCTION PROCESS 
AND EQUIPMENT AREA 

of the 

LOW-COST SOLAR ARRAY PROJECT 


Contract JPL-954796 
Annual Report Oct. 1978 to Oct. 1979 
(DRD Line Item 7) 


October 1979 


H. Goldman and M. Wolf 


The JPL Low-Cost Silicon Solar Array Project is sponsored by the 
U.S. Department of Energy and forms part of the Solar Photovoltaic 
Conversion Program to initiate a major effort toward the develop- 
ment of low-cost solar arrays. This work was performed for the 
Jet Propulsion Laboratory, California Institute of Technology by 
agreement between NASA and DoE. 


This report was prepared as an account of work sponsored by the 
United States Government. Neither the United States nor the 
United States Department of Energy, nor any of their employees, 
nor any of their contractors, subcontractors, or their employees, 
makes any warranty, express or implied, or assumes any legal 
liability or responsibility for the accuracy, completeness or 
usefulness of any information, apparatus, product or process 
disclosed, or represents that its use would not infringe privately 
owned rights. 


ii 


Summary 


This Annual Report summarizes the analyses performed and 
conclusions reached during the year from October 1978 to October 
1979. This work was largely contained in three topical reports 
which were issued as Quarterly Reports. In addition, material 
is presented in this report which was covered in Monthly Reports, 
but which was not included in the Quarterly Reports. Much of 
this material makes up chapter V, "Energy Analysis." 

The report is arranged by topics, with the analysis of 
slicing processes covered in chapter II and of junction formation 
processes in chapter III, with chapter IV containing the des- 
cription of a simple method for evaluation of the relative 
economic merits of competing process options with respect to the 
cost of energy produced by the system, and chapter V the energy 
consumption analysis. 

Subsequent to the analysis of the Czochralski crystal pulling 
process performed in the preceding year, the important companion 
to all ingot processes, slicing, was examined. Progress is 
being made or projected towards reduced kerf and decreased wafer 
thickness, together with higher throughput rates, on all approaches 
to slicing, including the generic types of fixed abrasive (diamond 
ID saw, wire saw) and slurry sawing (reciprocating multi- 
blade or multi-wire) . From their current position of comparable 
add-on price, the projected improvements, if successfully carried 
out, could also yield comparable price reductions. However, all 
these advancements cannot eliminate the kerf losses, with the 
waste of valuable material. 
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The analysis of the junction function process contains an 
accumulation of the pertinent technical and economic data on 
current processes, such as the Spectrolab gaseous diffusion 
process for front junction formation, and ion implantation 
using the Varian-Extrion 200-1000 machine. These recent exper- 
ience data were used to assay the projections to future improved 
process methods, such as diffusion processes proposed by Motorola 
and RCA, and ion implantation advances proposed by Lockheed, 
Motorola, RCA and Spire. 

The analysis of the energy consumption in the solar module 
fabrication process sequence, from the mining of the SiC^ to 
shipping, shows, in the current technology practice, inordinate 
energy use in the purification step, and large wastage of the 
invested energy through losses, particularly poor conversion 
in slicing, as well as inadequate yields throughout. The cell 
process energy expenditures already show a downward trend based 
on increased throughput rates. The large improvement, however, 
depends on the introduction of a more efficient purification 
process and of acceptable ribbon growing techniques. 
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I. INTRODUCTION 


The manufacturing methods for photovoltaic solar energy ilization 
systems consist, in complete generality, of a sequence of individual pro- 
cesses. This process sequence has been, for convenience, logically seg- 
mented into five major "work areas": Reduction and purification of the 
semiconductor material, sheet or film generation, device generation, module 
assembly and encapsulation, and system completion, including installation 
of the array and the other subsystems. For silicon solar arrays, each 
work area has been divided into 10 generalized "processes" in which certain 
required modifications of the work-in-process are performed. In general, 
more than one method is known by which such modifications can be carried 
out. The various methods for each individual process are identified as 
process "options". This system of processes and options forms a two- 
dimensional array, which is here called the "process matrix". 

In the search to achieve improved process sequences for producing 
silicon solar cell modules, numerous options have been proposed and/or 
developed, and will still be proposed and developed in the future. It is 
a near necessity to be able to evaluate such proposals for their technical 
merits relative to other known approaches, for their econc.iic benefits, 
and for other techno-economic attributes such as energy consumption, and 
generation and disposal of waste by-products, etc. Such evaluations have 
to be as objective as possible in light of the available information, or 
the lack thereof, and have to be periodically updated as development 
progresses and new information becomes available. Since each individual 
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process option has to fit into a process sequence, technical interfaces 
between consecutive processes must be compatible. This places emphasis 
on the specifications for the work-in-process entering into and emanating 
from a particular process option. 

The objective of this project is to accumulate the necessary infor- 
mation as input for such evaluations, to develop appropriate methodologies 
for the performance of such techno-economic analyses, and to perform such 
evaluations at various levels. 

In evaluating even current processes substantial gaps or uncertainties 
were found in important information required for both technical and economical 
evaluation of the currently practiced processes. In proceeding to the 
evaluation of processes which are still in the developmental or even con- 
ceptual stage, the gaps in needed information become very large. In these 
cases, it is necessary to fill the gaps more extensively with estimates based 
on extrapolations or analogies. Such estimates always leave some doubt on 
the accuracy of the evaluations, and it will be necessary to also make 
"probable error" estimates to reduce decision mistakes based on early 
evaluations. Nevertheless, collecting the information and carrying out 
evaluations at the earliest possible time provides not only a planning tool, 
but also aids in uncovering the deciding attributes about which information 
ought to be obtained at an early stage of the development process. 

This annual report describes the work completed during the last 12 months, 
and summarizes the work in progress. In the preceeding 12 month period, the basic 
methodologies for performing the comparative analyses of competing process 
options were developed, as well as the formats for accumulating the needed in- 
formation. Also, the processes for the reduction of quartzite to silicon and 
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for Czochral ski crystal growth were studied. The reduction process was investi- 
gated primarily from the energy consumption viewpoint, as the process is 
already carried out cost-effectively on a large industrial scale ( -~ 350,000 t /y 
worldwide, > 140,000 t /y in the U.S.A.). During the current 12 month period, 
the analysis of the slicing and junction formation processes was completed 
and is described in detail in this report. The processes for metallization 
(contact formation) and anti-reflection coating are at the state where the 
information which is available in contract reports and published literature has 
essentially been accumulated, and the analysis has been started, with considerable 
gaps in the information available in the reports becoming apparent. Finally, 
the methodology for the comparative evaluations has been refined and completed, 
based on the cost of the energy produced by the system rather than 

independent cost per unit peak power output value. 

This report is arranged in chapters according to the subject items 
investigated. The detailed process-data formats to the individual 
process options which were included in the quarterly reports, have been 
omitted since they would have added greatly to the volume of paper used. 

It should also be noted that the subject items were studied during a 
certain time period and the Conclusion reached then, with the subsequent 
effort devoted to a different subject area. In this approach, the analyses 
are not continuously updated by inclusion and reeValuation of new develop- 
ments. It may also be noted that throughout this report, costs or prices 
are quoted in 1975 dollars, unless, in rare special cases, it is stated 
otherwise. 
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II . Slicing 


1. Introduction 

The fabrication process sequences for solar modules involves many steps, 
the principal ones of which up to crystal growth have been analyzed previously. 
The next major process following Czochralski crystal pulling, or any other 
process that results in bulk ingots, such as the casting of semi crystal 1 ine 
ingots or the Heat Exchanger Method of crystal growth, has to be the division 
of these ingots into wafers, commonly referred to as "slicing". The first 
application of this developing methodology was made to the Czochralski 's 
crystal pulling process. 

Previously, we had examined the reduction of quartzite to metallur- 
gical grade silicon and did a comparative evaluation of competing Czochralski 
techniques for growing single crystal, cylindrical ingots. The next major 
process step in the sequence for producing single crystal silicon wafers, 
today and in the near future (up to 1982), is the slicing technique. The 
evaluations were started with the current methods of multiblade slurry 
slicing, and inner diameter slicing using a diamond coated blade for which 
a large amount of the needed information is available. 

We have tabulated production experience data obtained from Spectro- 
lab^ for slicing 2-cm rectangular, b.4-cm and 7.6-cm diameter wafers 
using the Varian multi-blade slicing system, and similar data obtained 
from HAMCO^, for ID slicing of 10.16-<m diameter ingots using their 
equipment. Experimental data from OCLI^, Varian^ and TI^ for 
multiblade watering, from OCLI^ and STC^ for ID slicing, and from 
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jPl( 8 ) f 0r the Yasunaga multi-wire slurry slicing system, were also tab- 

(91 

ulated. To complete the analysis, projections made by Varian' ' for multi- 
blade slicing, by STC for ID slicing^ by Crystal Systems^*^ for their 
fixed abrasive multi-wire system, and by Solarex^ 1 ^ for the Yasunaga 
multi-wire slurry system were examined. 
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Brief Descriptions 
of the Slicing Technique 


2 . Brief Descriptions of the Slicing Techniques 
2.1 Multi blade Slicing 

The multiblade slurry sawing method is one of the two tech- 
niques used in current production slicing. In its present configuration 
230-250 blades of 38-cm length of hardened 1095 steel are mounted and 
evenly spaced on a blade head that is, for slicing, reciprocated, at 
frequencies below 2 Hz (normally about 1.6 Hz), across the workpiece using 
approximately a 20-cm stroke. The abrasive slurry is pulsed sprayed or, 
at times, dripped onto the top surface of the workpiece and recirculated 
by a pump. The slurry is a Si C abrasive suspended in PC oil. It is nor- 
mally used for one load before it is discarded. There are no practical 
ways, at present, to re-use the abrasive slurry for more than one load. 

The current multiblade slicing machines can accept blade heads 
up to 18.5-cm wide. However, the number of blades in a blade head, and 
consequently, the number of slices that could be produced per load, is not 
limited by the blade head width per se, but rather by the maximum tension 
force the blade head can exert on the blades. This is about 401,800 N 
for current production blade heads^. An adequate saw force commonly 
called “blade load", is necessary to achieve economically acceptable cut- 

(i,) 

ting rates in the slicing process. A blade load of about 1-2 N/blade ‘ , 
is usually applied. Excessive blade loading, and even normal loading 
after some blade wear, can cause deflection of the blades, often called 
"buckling", which results in inaccurately sliced wafers or even broken 
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wafers. To minimize buckling, the blades need to be stressed as much as 

possible, which, in current practice, is 80% of the yield strength of 

(5) 

1095 steel, or 1.37 GPa' . Therefore, the maximum number of blades per- 
mitted per blade head is 401 .8/1 . 37* A , where A is the cross-sectional 
2 

blade area in mm . For a 6.35 mm high blade, 0.20 mm thick, a size that 

( 4 ) 

is normally used in production' ' , the maximum number of blades thus is 

230. Reducing the blade thickness to 0.15 mm will increase the maximum 

number of 6.35 mm high blades to 307. At present, the thicker 0.20 mm 

blades are used in production because of their better wafer yield, as 

they are less susceptible to buckling which can be caused by vertical 

misalignment at the beginning of the slicing process and by increased 

blade tension, resulting from a reduced crossection because of blade wear 

f 51 

near the end of slicing' '. 

There are two types of blade packages available: the drill -pin 

package and the epoxy package. In the former, the alternately arranged 

blades and spacers which determine the thicknesses of the kerf and wafers are 

held together by four threaded rods. It is the cheaper of the two types of 

package ($50 compared to $175), but often requires additional alignment 

(31 

before mounting on the slicing machine' . In the epoxy package, an 

adhesive is applied between the spacers and the blade ends to hold the 
(41 

package together' . 

The production procedure for multi-blade slicing involves first 
mounting the workpiece, or silicon crystal, with wax, epoxy, or other 
suitable cement on a graphite or ceramic base plate. The workpiece is 
then clamped by the baseplate to the slicing machine. To help increase the 
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yield, ceramic bars are often similarly cemented longitudinally onto the 
cylindrical crystal near its top and bottom horizontal tangents. The bars 
"smooth-out" the slicing by decreasing the variation in kerf length 
and blade load as the blades travel downward through the cylindrical 
crystal. In addition, ceramic bars near the top tangent minimize the 
effect of vertical misalignment by reducing blade buckling by the time 
they enter the silicon crystal. Those bars near the bottom, help to 
smooth the transition of the blades cutting into the base material by 
equalizing the slicing properties above and below the crystal to base 
transition. Some of these benefits are also obtained, in some places, with- 
out the use of ceramic bars by varying the blade load according to the 
changing kerf length during the slicing process. After the slicing is 
finished, the wafers, still attached to the base, are removed from the 
slicing machine and the wafers are then detached from the base. 

The effective linear cutting rate of the multiblade process is 

presently about 550 times smaller than the ID diamond saw. The linear 

cutting rate cannot be increased significantly because of the limit on the 

blade load and because of the blade head mass which limits the reciprocating 

frequency. The blade load cannot be increased much beyond its present 

value without significantly increasing blade buckling since the tensile 

strength of the blades is fixed. Varian found that a blade load of 2.77 N/ 

blade caused severe enough buckling to separate the crystal from its 
(41 

mount' . In another experiment, a reciprocating frequency increase to 

(41 

2 Hz resulted in sufficient vibration to break .ill wati-rs Iheretorc, 

in order to increase the throughput rate, or the wafer area produced in 
the multi-blade slicing process per unit time, either the number of slices 
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in the load, or the area yield per load, has to be increased without 
significantly increasing the time of the run. The area output per load 
can be increased in a combination of several ways: by increasing the num- 
ber of blades per unit blade head width, as can be achieved by decreasing 
the blade and/or spacer thickness; by increasing the width of the blade 
head without changing blade and spacer thicknesses; or by increasing the 

width of the workpiece. 

The blade thickness has a lower bound set by its strength. If 
the blade is too thin, it will buckle under the blade load, or break 
from the blade tension, resulting in broken wafers and low yields. 
Reduction of the spacer thickness is limited by the wafer strength. 

Slicing wafers too thin increases their chance of breakage due to pressure 

from the lateral blade movement, blade vibration, blade buckling etc. As 

the blade and spacer thicknesses are decreased, the increased fragility 

of the blades and the wafers ultimately leads to significantly lowered 

yields. Experimentally, Variant has found that using 0.15 mm thick 

blades with 0.30 mm spacers still results in good yields. Under these 

conditions 0.25 mm thick wafers with 0.20 mm kerf are produced. This 

gives, assuming a wafer yield of 95a, which has been demonstrated by 

2 

Varian, an area conversion ratio of 0.9 m /kg-Si which is a 50!?. improve- 
ment over Spectrolab's recently experienced area conversion ratio in 
slicing 5.4-cm and 7.5-cm diameter wafers. 
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Varian Is also currently experimenting with a larger blade head 
width that can accept 900 to 1000 blades. This blade head weighs approxi- 
mately one ton. Therefore, the workpiece will be reciprocated against 
the stationary blades. The workpiece size is projected to be 12-cm in 
diameter and 40.5-cm long yielding a wafer area of 9.67 m^/load using the 
900-blade machine with the aforementioned blade and spacer thicknesses. This 
area yield is over four times higher than obtained in present commercial 
practice. 


A third method to potentially increase the area yield per load with- 
out increasing the slicing time would be to increase the width of the 

workpiece, or the kerf length, by slicing two or more ingots, placed 

(51 

side-by-side, simultaneously. TI' ' has found that the machine slicing 
time, and, correspondingly, the linear cutting rate, is essentially in- 
dependent of the kerf length. TI has therefore proposed slicing two 12-cm 
diameter ingots at one time to increase the multi-blade slicing productivity. 
The area yield per load, with details of this projection given in Tables 
I to III, can thus be doubled without significantly changing the slicing 
time. 

2.2 Inner Diameter Slicing 

In the process of inner diameter, or ID, slicing, one wafer is 
sliced at a time with a rotating, diamond impregnated blade. The rotation 
speed depends upon the blade size, and is 2,100 rpm for a bl? with a 
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15.25-cm diameter hole, and 1650 rpm for a 20.32-cm diameter, inner dia- 
meter blade. The blade consists of a stainless steel core which is 0.10 
and 0.15 mm thick for 15.24 and 20.32-cm blades, respectively, with dia- 
mond plated edges. The total thickness of the 15.24-cm blade is approxi- 
mately 0.30 mm, and the 20.32-cm blade is about 10% thicker. The blade 
is mounted around its rim in a vise-like holder where hydraulic pressure 
is applied to tension it radially. 

The linear cutting rate, or the rate that the inner diameter blade 
traverses the silicon can be up to 305 cm/h, or almost three orders of 
magnitude higher than for the slurry, multi-blade process. There are 
several reasons for this. First, the inner diameter blade speed is approxi 
mately 1,600 cm/sec as opposed to less than 80 cm/sec for multiblade 
slicing. Therefore, the contact length per unit time between the blade 
and the silicon for ID slicing is twenty times higher than for multiblade 
slicing. Also, fixed abrasive slicing removes more kerf in a unit contact 
length because there are two surfaces moving relative to each other instead 
of three as in slurry slicing. In slurry slicing, the abrasive is pushed 
into the workpiece and is "rolled out". Whereas for fixed abrasive slicing, 
the abrasive cuts into the workpiece to remove the kerf. Finally, the 
diamond plated layer on the ID blade increases the blade's rigidity and 
thickness and allows the application of more force, by the blade, on the 
workpiece than in multiblade slicing. The total thickness of the ID 
blade is 300-330 pm thick while the multiblade is 150-200 pm thick. It 
should be noted that the effective ID cutting rate is about 10-207 lower 
than indicated by the blade's linear cutting rate because of the 18 to 24 
seconds between two consecutive slices, when the blade is returning to its 
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original vertical position and the silicon crystal is being indexed. 

In mounting the ingot, one end is attached to a graphite base with 

epoxy and the ingot is then placed in a box with rubber supports along it's 

length to keep it rigid. The stiffness of the mount will affect the 

vibration level between the blade and workpiece, influencing the wafer 

(31 

thickness and yield' . At present, ID machines can accommodate ingots 
(2 3) 

up to 50-cm long' * . The current practice of slicing 10.16-cm diameter 

2 

wafers, 0.50 mm thick with a 0.33 mm kerf, yields a area of 4.8 m /load 
2 

or 0.50 m /kg, a: a practice, wafer yield of 98%. During slicing, either 
water or water mixed with a small percent of Rust-Lick is sprayed on the 
cutting edge, at a rate of about 2 ms^sec, to cool the blade. The blade 
must be dressed, every 50 slices for the 15.24-cm blade and every 25 slices 
for the 20.32-cm blades for proper slicing, in order to remove dirt and 
expose a fresh cutting surface. The dressing is done with 5 cuts of an 
alumina stick. The lifetime of the blade is dependent on the rate of 
diamond "pull-out" and the degree of metal fatigue and varies quite ex- 
tensively from blade-to-blade. The lifetime median is about 3,000 7.52-cm 
diameter slices for the 15.24-cm blade and 5,000 10.16-cm diameter 
slices for the 20.32-cm blade. 

A method being investigated, to increase the ID saw's productivity 
by a factor of two, is crystal rotation^ 7 '’, The cutting speed is doubled 
using a rotatina crystal since the blade has to traverse only half-way 
through the crystal diameter. The half penetration in rotating crystal 
slicing permits the use of a cheaper, smaller diameter, and thinner inner 
diameter blade, lor slicing 10-cin diameter waters with this technique the 
wafer thickness and kerf are expected to be 225 pin and 210 pm respectively^. 
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This process is expected to be In conroercial use by 1982. 


2 . 3 The Yasunaga YQ-100 Multiwire Saw System 

The Yasunaga multiwire saw is a slurry slicing system which 
uses a single wire (600 to 30,000 m in length) routed around a rocker arm 
tensioning device, a wire guide catridge, and a take-up reel. The con- 
tinuous wire forms up to 250 multiple loops around the three grooved wire 
guides, arranged in an equilateral triangle, that are the key parts of 
the wire guide catridge. During slicing, the wire guide catridge oscillates, 
while the workpiece is raised against the wires with a preset force. An 
abrasive slurry is sprayed on the cutting surface. The procedure for 
mounting the silicon crystal for multiblade slicing is similar to that 
described for multiblade slicing. 

The chief potential benefit of the Yasunaga saw is its high 

area-mass conversion ratio by employing closely-spaced, small diameter 

wires. The current YQ-100 model has a workpiece capacity of 10x10x10 cm 

( 8 ) 

and as demonstrated by experiments,' ' it can slice 215, 212+7 um thick 
wafers with less than 200 pr- i ar f using 0.4 mm pitch guides, 0.16 mm 
diameter wire and 13 pm SiC abrasive. Under those conditions an area to 
unit mass ratio of 1.04 m /kg is obtained, which is about 50t higher than 
what any other current production or experimental slicing system achieves. 
This higher area to mass ratio effectively recjces the consumption of 
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single crystal silicon, to produce a given wafer area, by a third. It 

is projected that the Yasunaga saw can achieve an area-mass ratio of 1.42 
2 

m /kg by employing closer spaced pitch guides (0.3 mm), smaller diameter 
wire (0.08 mm) and a finer abrasive ( 5 urn). This would yield a 200 pm 
thick wafer with 100 ym kerf^ 1 ^. 

It is believed that the narrow lapping band of the wires of 
the Yasunaga saw results in wafers with less subsurface damage than with 
other commercial slicing techniques^ 11 and this is being investigated^. 

Currently, the Yasunaga saw is not used for the production of 
silicon wafers, at least not in the USA, although Solarex hr . 'cently 
obtained a machine for pilot line operation. However, Motorola is using* in 
its semiconductor device production an in-house developed, pi^y* -‘"ory 
wire saw with capabilities which seem to be comp; able to that ci the 
Yasunaga saw. 

2 . 4 The Multiwire Fixed Abrasive Slicing Technique ("FAST") 

This method is similar to multiblade slicing, except that the 
silicon is sliced with d iamond- impregnated wires instead of steel blades 
and an abrasive slurry. In FAST, the diamond impregnated wires are 
mounted and evenly spaced, at a linear density expected to be up to 
25 cm -1 , on a light weight frame that is reciprocated across a rocking 
workpiece^ 1 ^. The wires are coated with 22 to 45 ym diamonds imDedded 
in a metal matrix, and can be coated on their bottom halves only 
to reduce abrasive costs. Development is still proceeding towards finding 
an optimum wire composition, but it has been found that, heat-hardened, 
tungsten core wire, diamond-impregnated, and nickel -plated, has a good 
lifetime, which means it could be used for about 10 loads before signifi- 
cantly losing its cutting ability^* 1 ^. 
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Crystal Systems has conducted most of their experiment^ per- 
taining to FAST, on a modified Varian 686 wafering machine. Consequently, 
the slicing potential of multiwire, fixed abrasive slicing has not been 
fully demonstrated. For example, workpiece size has been, for most of 
the experiments only 4 x 4 cm, and the reciprocating rate lower than re- 
quired for optimum fixed -abrasive slicing. A slicing machine, built to 
Crystal Systems' specifications, have just been delivered to them 
and slicing with this machine has just been initiated. The new slicing 
machine has been designed to provide higher cutting rates and lower wafer 
and kerf thicknesses and operate with a much lighter blade carriage, at higher 
reciprocating frequencies, and reduced vibration than the Varian machine. 

It is expected that this multiwire, fixed abrasive slicing technique could 
have a cutting rate of 0.6 cm/h (twice the value previously achieved with 
good yields), with an area to mass ratio of 1.1 m /kg by producing wafers 
200 ym thick with a 200 ym kerf. 

The add-on prices for "FAST" have been projected for 1986 since the 
state of development of the system and the comparatively small base of 
experimental data available, making it unlikely that this slicing tech- 
nique could be in significant commercial operation by 1982. 

3. TABULATION OF OPERATION. LABOR, MATERIAL AND COST DATA 

Tables II-l to I 1-3 summarize the data provided by various organizations 
for the slicing techniques that are being used or developed. Included in 
these tables are production experience data from Spectrolab^^ for multi- 
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TRBLt H-l* 


Organlrat ion 


SLICING OPERATION DATA fCR NL'ITIBLADE WAl'ERING 


Specttolab OCLI Varian TI 

(Production Experience) lExper rr'cntel ) teeperimontl . ('300 blade rxperimenta) 

2 c» | 5.4 cm * 7.5 cm 10.16 rm no. P-005) I (Pro jectionll projection) inrl. Projection 

Rectangular | Diameter I Diameter Diameter 10cm Diameter! lOcm Diameter] 12cm Diameter 12cm Diameter 


16 cm long 16 cm long 15 cm long 11.7 cm long 13.5 cm Ion? 40.5 cm long I 2,13 cm long 

I inqote 

3 2 1 1 1 1 I 2 


' severe 
| j- slight 


600 grit 
SiC 
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TABLE II-1B 


SLICING OPERATION DATA FOR MULTIWIRE AND INNER DIAMETER WATERING 


1 


Multiwire Wafering 


Inner Diameter Slicing 



Crystal Systems 

( Yasunaaa YQ-100 1 

OCLI ! 

KAME0 

Organization 

Fixed Abrasive 
Method 
(projection) 

Experimental) 
7.6 cm 
diameter 

(Projection) 
10 cm 
diame ter 

{ experimental) 
7,6 cm 
diameter 

(Experimental' 

10,16 cm 
diameter 

(Production exp.) 
10 . 16 cm 
diameter 

i. 

Workpiece size 

30x10x10 cm 

10 cm long 

10 cm long 

50 cm long 

25 n long 

46 cm long 

2. 

No. of workpieces/ 
load 

1 

1 

1 

1 

1 

1 

i. 

Slices/load 

250 

215 

333 

725 

350 

555 

4. 

Wafer thickness 

(n») 

0.1 

0.21 + 0.01 

0.2 

C.36 + 0.02 

0.36 ♦ 0,02 

0.50 

5. 

Kerf thickness 

0.3 

0.2 

0.1 

0.33 

0.35 

0.33 

B 

Practical Wafer 
Yield 

1.00 

1.00 

1.00 

0.95 

1.00 

0.98 

| 

Fraction Silicon 
Incorporated in 
Wafer 

0.25 

0.51 

0.67 

0.50 

0.51 

0.59 

B 

Depth of Surface 
damage (um) 

Fissures ex- 
tend 3 um 

'15 

-6.5 

n.a. 

n.a. 

n.a. 

F 

Abrasive 

none 

GC 1200 
(13 Um) 

5 urn Sic 

none 

none 

none 

1 

Vehicle or 
coolant 

1:1 water: 

ethylene 

glycol 

lapping 

oil 

n.a. 

80:1 water: 
rust lick 

80:1 water: 
rust lick 

water 

9 

Concentration 
(kg /l) 

- 

-1.5 

n.s. 

- 

- 

- 

iTw 

Flow rate 

a/h) 

n.a. 

3600 

3600 

7.2 

8.4 

n.a. 

13. 

Type of blade 
or wire 

Ni plated, 
tungsten wire, 
liaxoond im- 
pregnated 

Steel wire 

Steel wire 

Model STC-16 

p blade, 

diamond 

plated 

Model STC-22, 
ID blade, 
diamond 
plated 

ID blade 

diamond 

plated 

14. 

Blade or wire 
dimensions 

0.125 mm 
;ore 0.25 mm 
total diameter 
15 um diamonds 

0.16 mm dia- 
meter 0.4 m 
pitch 
guides 

0.08 mn 
diameter, 

0.3 mm pitch 
guides 

42.23 cm OD 

15.24 cm ID 

0.10 mm thick 
core, 

0.28-0.30 total 
thickness 

55.88 cm OD, 
20.32 cm ID, 
0.15 mm thick 
core, 0.33- 
0 . 36 tota 1 
thickness 

n.a. 

15. 

Amount on 
machine 

250 wire blade 
package 

-17,000 m 

-35,000 a 

1 

1 

1 

16. 

No. of loads 
before blade 
change 

9 

3 

3 

4.1 

14.3 

1 

I 

Wa^er area/load 
(a 2 ) 

7.50 

0.98 

2.62 

3.14 

2.84 

4.41 

1 

Area yield 
(m 2 /kg) 

1.1 

1.04 

1.42 

0.59 

0.60 

0.505 

1 


0.6 

0.84 

0. 3 

305 

305 

305 

20. 

Slicing time 
segment/load (h) 


9.0 

30.0 

23.9 

14.7 

23.12 

21. 

Load/Unload 
time (h/load) 


n.a. 

n.a. 

1.23 

0.735 

0.083 

22. 

Cutting tool 
change, machine 
service (h/load) 

■■ 

n.a. 

n.a. 

1.02 

0.84 

0.33 

■ 

Machine augment 
time (h/load) 

18.0 

10.0(e) 

3.1(e) . 

26.2 

16.3 

23.5 


Machine product- 
ivity (m 2 /h) 

0.42 

0.098 

0.085 

0.126 

0.176 

0.19 
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( 2 ) 

blade slicing and from HAMCO for ID slicing, and experimental results 
for multiblade slicing, from OCLI^, Varian^ and TI^»for multiwire 
slicing from JPL^ and ID slicing from OCLI^. In addition, projections 
made by Varian for multiblade slicing^ , by Crystal Systems^^ for 
their "FAST" method, and by Solarex^^ for the Yasunaga saw are included. 


The operation data for multiblade slicing are listed in Table II-1A, 
while Table II-1B contains the corresponding data for the fixed abrasive and 
slurry multiwire and the inner diameter slicing processes. These tables 
contain the process attribute of slicing which are summarized on Figure II.l 
The first two lines of Table II-l are the dimensions of the workpiece and 
the number of workpieces per load, the product of which is the slicing 
machine's capacity. The wafer area produced in a load is related to the 
workpiece capacity through the wafer and kerf thicknesses and practical 
wafer yield. This wafer area per load (Table II-l, line 17) can also be 
calculated as the product of the theoretical number of slices cut per 
load (Table II-l, line 3), the "practical wafer yield" (Table II-l, line 5) 
and the area of the single wafers. The "practical wafer yield" fraction 
is the number of acceptable wafers divided by the theoretical number 
sliced per load. The wafer area per unit mass (Table II-l, line 18) is 
calculated by dividing the practical wafer yield by the product of the 
sum of the wafer and kerf thicknesses (lable n-l, lines 4 arid b) and the 
density of silicon, or 
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10 * Jf-1.6 


II-l. 18 


m 2 /kg 

(II-l .4+11-1 .5)*2.34 


where II-l - n represents the value from Table II-l , line n. 

The wafer thickness, kerf and practical wafer yield are necessary 
for finding the division of the input silicon crystal or workpiece into 
the silicon incorporated in the work-in-process wafer (Table II-l, line 7) 
and that silicon lost in kerf and broken wafers. 


The procedures for determining the subsurface damage depths, listed 
in line 8 of Table I, were not consistent between organizations. The 
most accurate method for determining subsurface damage depth is to re- 
move wafer surface material until the cell efficiency becomes independent 
of any further removal. Spectrolab's values reflect this procedure^ K 
The other subsurface damage depths were determined by chemical etching 

to remove surface material followed by Wright etching to reveal defects^, 

(51 

by etching and x-ray topography' , and by. angle lapping and Sirtl etch- 
ing^. 

Indirect material requirements, briefly summarized on Figure 1 1. 2, in 
terms of the abrasive and vehicle, or coolant type, the slurry concentra- 
tion and its flow rate or that of the coolant, are listed in lines 9-12 of 
Table I. Lines 13-16 describe the expendible tooling requirements such 
as the type of blade or wire, its dimensions, the size of the blade pack 
and its life expectancy. These data are necessary for determining the 
expendible tooling and material costs. 

The effective cutting rate (Table II-l, line 19) is defined here as 
the workpiece diameter divided by the slicing time segment, which is the 
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time the machine is actually sawing (Table H-l, line 20). The time 
periods when the machine is not actually slicing and cannot be used for 
slicing because of preparatory or service operations, are listed in lines 
21 and 22. The sum of these lines and the slicing time segment is the 
machine segment time (Table H-l. line 23), or the average time needed for 
slicing a load, including loading, unloading and servicing. The machine 
segment time is needed for calculating the number of loads processed 
annually, and the machine productivity (Table H-l» line 24) which is the 
wafer area sliced in a load divided by the machine segment time. 

The requirements per machine load for labor, included that needed 
for service and repair, for indirect material needs, including electricity 
consumption, for capital expenses, which consists of machine and facility 
components, are included in Tables II-2A and II-2B. These data form the basis 
for calculating of the manufacturing cost components of labor, expendable 
tooling, indirect materials, and capital. Also listed in these tables 
are values necessary for calculating direct material or silicon costs: 
the proportion of silicon lost in grinding the cylindrical ingots to a 
uniform diameter, the unit mass of silicon incorporated in the wafer and 
that lost in kerf and broken wafers. 

The labor times required for each part of the crystal slicing 
operation (see Fig. 2), that is crystal mounting, machine loading and 
machine monitoring are listed in lines 1-3 of Table II-2, with their total 
on line 4. The service labor time, which includes changing the blades or 
wires, is listed in line 5. 
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Expendable tooling and indirect material requirements, in terms of 
the blade or wire set costs and the quantities of vehicle or coolant and 
abrasive consumed during a run, are listed in lines 6-8 of Table 1 1 -2 . The 
electrical consumption for a run (Table II-2, line 10) is considered as an 
indirect material and is obtained by multiplying the slicer's power re- 
quirements by the slicing time segment (Table II-l, line 20). 

In order to calculate the potential number of loads that can be 

sliced annually, shown in line 12, the machine segment time (Table II-l, 

line 23) is divided into 8280. This last value, 8280, is taken from 
(121 

SAMICS and is the number of annual hours the wafer slicing plant 
operates. The plant operation schedule is continuous except for one 
1-week vacation, two 4-day weekends, and one 3-day weekend, and was 
chosen to maximize annual production by minimizing slicer shutdowns 
during a run due to plant closings. 

After dealing with expenses, the sum of the machine and facility 
costs, or the capital cost portion of the manufacturing costs needs to 
be considered. The capital costs are dependent on the factors listed 
on Figure II. 3. The annual machine cost (Table 1 1-2, line 14) is the product 

of the initial cost of the slicing machine, including installation, taken 

from the data sources, and the standardized charge rate of 0.2135 y”"* . 

( 12 ) 

This charge rate was taken from SAMICS' , using a depreciation sched- 
ule of 7 years, a state tax of 2% on one-half the capital, a 4% insurance 
premium, and a 12% interest-on-debt rate on one-twelfth the initial 
capital cost. The low ratio of dept to capital, or the low financial 
leverage, is due to the postulate that the photovoltaic industry would be 
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IABOR AND INDIRECT MATERIAL. SLICING REQUIREMENTS 


LABOR TIMES: 

ATTACH SUPPORT BLOCK TO INGOT 
MACHINE LOAD/UNLOAD 
MACHINE MONITORING 
TOOL CHANGE/MACHINE SERVICING 


INDIRECT MATERIAL COSTS: 

SLURRY (COOLANT) TYPE 

UNIT COST 
USAGE 

TOOL (BLADE) TYPE 
COST 
LIFE 

MACHINE REPLACEMENT PARTS 
PURCHASED MACHINE SERVICING 
MISC. (MOUNTING BLOCKS, ADHESIVE) 
ENERGY 


Figure II. 2. 
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CAPITAL COSTS 


MACHINE COST 
(MACHINE LIFE) 

ALLOCATABLE BUILDING AREA 
(SPECIAL SERVICES) 


Figure II. 3. 
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unable to raise large amounts of debt capital, without large Interest rates, 
because It will be a rapidly evolving industry with appreciable risks^^. 

The second capital cost contribution comes from the building. The 

allocatlble building area, shown in line 15 of Table II-2, was taken, accord- 
(121 

Ing to SAMICS , as twice the machine's operating area. The doubling 
accounts for Indirect and overhead space needed e.g., for functions such 
as maintenance, administration and receiving/inventorying, as well as for 
aisles, washrooms, etc. The initial building cost (Table II-2, line 16) is 
taken as $1506. 95/m^, according to SAMICS^^^ , and is based on the machine 
operating area only. This cost figure includes appropriate cost allocations 
for the additional building space needed as outlined above. The facilities 
charge rate used to calculate the annual building cost (Table H-2, line 17). 
from the initial cost, is 0.117 y' 1 . This value was obtained in the same 
fashion as the equipment c large rate, except that a 40-year life expectancy 
is employed for determining the depreciation rate of the building. Also 
a 31% surcharge on the annual cost of capital is included, in the 0.117 y ^ 
factor, to account for special services which are the "indirect" utility 
consumption, that is for heating, air-conditioning, 1 ighting.etc. for the 
building. 

To properly calculate the direct material cost, that is the cost of 
the cylindrical slicing ingot, the amount of the silicon crystal lost in 
grinding is necessary. The grinding of the cylindrical ingots to a uni- 
form outside diameter, previous to slicing, facilitates the slicing oper- 
ation, as well as tooling and handling of the sliced wafers in subsequent 
device fabrication procedures. In calculating the mass fraction of sili- 
con lost in grinding, shown in line 18 of Table II-?, the average diameter 
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loss Is assumed to be 0.6 cm. With this diameter loss, and the consequent 
loss of mass, the price per unit mass of silicon entering Into the slicing 
operation can be determined. Since the grinding diameter loss stays con- 
stant with crystal diameter, the fraction of lost silicon is inversely pro- 
portional to the diameter of the crystal. 

The difference between the add-on processing cost and the work-in- 
process cost Is the cost of the direct material contained in the wafers. 

The latter value for a unit area can be obtained by multiplying line 19 of 

Table I I - 2 by the unit mass silicon r- st. To obtain the amount of silicon 
contained in a unit wafer area, the incorporated silicon fraction is di- 
vided by the wafer area per unit mass (Table II-l, line 18). The incorporated 
wafer fraction is the product of the yield fraction (taken from Table II-l, 
line 6) and wafer thickness (Table II-l, line 4) divided by the sum of the 
wafer and kerf thicknesses. In equation form, the fraction of silicon con- 
tained in the wafer is, 

I r 2 19 - H-1.6MI-1.4 kg 

' (II-l .4+11-1.5)^11-1 .18 7 ’ 


with the roman numerials dash arabic numbers representing the table numbers 
and the final arabic numbers, the nne numbers for that table. The kerf 
and broken wafer loss, recessary for differentiating the operating add-on 
cost from the specific add-on cost, is calculated in a similar fashion to 
line 19 of Table II, except that the kerf loss is represented by the kerf 
thickness and the broken wafer loss by the broken wafer fraction multiplied 
by the wafer thickness. Therefore 
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(II-1.5+(1-II-1.6)*II-1.4) 
(II-l .4+11-1 ,5)*II-1 .18 
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From the operation data and expenses, listed in the first two tables, 
th? add-on components of the slicing manufacturing slicing cost can be 
calculated. For the most part, the add-on cost components, shown in Table II- 
3, on a per unit area basis, are derived from the data of the proceeding 
tables using the relationships given in that table. The exceptions include 
the unit costs of the indirect materials which were taken from the sources 
footnoted in Table II-3. In addition, the purchased service cost fcr multi- 
blade slicing (Table IIt3, line 4), which includes the cost of machine 
maintenance and overhaul performed on the outside or under contract, used 
was $1529.3 y~^ and was obtained from Spectrolab^. HAMCO^ supplied the 
purchased service cost for an inner diameter slicing as $285.7 y \ The 

total material cost which is the sum of the first four lines of Table 1 1-3 

( 12 ) 

was increased by 5.26%, in accordance to SAMICS charge factors to 
account for handling and other miscellaneous expenses. 

The labor costs were calculated using the labor times, listed in Table 
1 1-2 and the labor rates shown in the Cost Account Catalog of the SAMICS 
Support Study^^. For calculating the direct labor costs which involve 
crystal mounting, machine loading and supervision the wages paid an elec- 
tronics semiconductor assembler, whose duties are described under SAMICS 1 

( 13 ) 

occupation classification no. 726884 and wages under catalog ■>. B3096D V ' 

were employed. The maintenance labor rate of a maintenance mechanic II 
(occupation classification no. 726884, catalog no. 83736D) was used to find 
the labor cost of internal machine service and cutting tool charges. The 
listed labor rates were multiplied by 1.432 to take into consideration fringe 
benefits, such as vacations, medical health plans, social security benefits, etc, 
and miscellaneous expenses. A surcharge of 25% was added to the direct 
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1. 

K* parti lb) a 
tooling 

(11*. 4 i 74. W* 
1A.17) 

lo. as 

15.34 

23.81 

74.31 

22.48 

10.49 

2.04 

1.3) 

J. 

luuiiala 

21. IS (a) 

8.95 (a) 

4.9S (a) 

14.19 (b) 

7.50 (e) 

1.40 (c) 

3. 65 (c) 

0.30 (d) 

1. 

electrical 
anerqry cost 
($0,032 *11*. 10 

I 1A-17) 

0.15 

0.41 

0.1) 

0.42 

0.47 

0.12 

0.14 

0.11 

4. 

Meplacanant 
ports $ pur- 
ctv**4d oonrico 

1.45 

1.01 

2.37 

a.»4 

4.62 

2.90 

0.47 

0.83 

i 

■ 

35.37 

29.19 

35.22 

98.96 

47.86 

18.01 

7.07 
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4.85 
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1* 31 

1.15 

7. 
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■M 
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4.90 

3.59 

2.43 

0.94 

0.94 

1. 

Other indirect 
labor (25% of 
4. ♦ 7.) 

■ 

6.53 

5.14 

2.J.5 

2.10 

1.44 

0.44 

0.52 

H 

Total labor 
(4. ♦ 7. » S.) 

19.32 

12.67 

25.70 

10.73 

10.50 

7.18 

2-33 

2.61 

I 
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4.94 

• .Of 

4.17 

7.90 

12.43 

7.80 

2.71 

2.22 

9 

PaclUtios coat 
(XU. 17 i XU .12 
•1-.17) 

1.14 

1.85 

1.46 

1.81 

2.85 

1.79 

0.41 

0.34 

I 

Capital Cost 
(10. 4 11.) 

6.10 

9.93 

7. S3 

0.71 

15.28 

9.59 

4.90 

2.56 

| 

Overhead 
(0.059 <• (10.1 
♦ 0.1:1 *ftl.» 

0.42 

0.6i 

0.93 

0.64 

1.04 

0.65 

0.20 

0.17 

i 
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i 
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95.93 

159.25 

114.29 

58.10 
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(a) P. Schmid and C.P, ►hat tack (“Heat Exchanger-Ingot Cast: mg/Si icing Process" 
Cryatal Systems, EROA/JPL 954 m-77/3, pp. 78-79 (10/77), 

(b) Calculated using $ 12 . 1 0/kq ' c 1 for the abrasive and 51. 25/t' “ for the PC oil 
and assuming the alurry ia used twice, 

(c) LSSA Project Report, "Multivire Slurry Watering Demonstrations, • Jet 
Propula ion Laboratory, DOE/JPL- 1012-781 7, (2/?8j, 

(d) Estimated f roe materials coat of Yeaunaqa'a 7.6 cm diameter ingot. 

(e) H.I. Yoo, "Aaseemont of Proaent State-of-the-Art Sawing Technology," OCLl, 
OOE/JPL 954BJ0-77/12, p. JB (12/77), 

l f) Estimated from oCLI’a material coat data, 

l g) Assuring total purchased service is $2,000 for the machine's lifetime. 
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labor and maintenance labor costs to account for the cost of supervisory, 
management, and other support personnel. 

The unit area equipment and facility costs, which constitute the 

capital cost, were obtained by dividing the respective annual costs by the 

annual area factory output. The overhead, listed in line 13 of Table II-3, 

is defined as the insurance, state taxes, and interest-on-debt payments on 

(121 

the working capital. As suggested by SAMICS' the working capital was 
taken as 15% of the equipment plus facility cost, or 15% of the capital 
cost. 

The profit and the amortization of one-time costs is represented by 

the return-on-equity (ROE), shown in line 14 of Table II-3. This value is 

equal to the SAMICS 1 return-on-equity (EQR), which is 20% of the equity 

(121 

portion of the book value' , plus the amortization of the start-up costs 
(AOC), minus the income tax investment credit (ITC) on 10% of the annual 
equipment depreciation divided by the product of one minus the federal in- 
come tax credit ( 1 - t) and one minus the miscellaneous expense fraction, 

(1 - x), or 


EQR + AOC - ITC ? 

ROE (II-3.14) = $/m • 

(1-x) * (1-t) 

The add-on cost components described above can be used to calculate 
a unit area wafer price that ignores the cost of the silicon ingots. This 
add-on price shown in line 15 of Table IJ-3, is the sun; of the -material , 
labor, capital, overhead and return-on-equity. To convert this value into 
a wafer price, the unit mass cylindrical crystal price, and the add-on 


11-30 


grinding cost must be added to it. The unground silicon crystal or ingot 
prices shown for 1970, 1982 and 1986 are taken from our previous evaluations. 
For 1978, the ingoc price is based on pulling 7.8-cm diameter ingots with 
a Leybold-Heraeus single charge puller. The silicon ingot prices employed 
for the years 1982 and 1986 are projections for mul ti-pull ing Cz-grown 
10.2-cm and 1 5.2-cn diameter ingots, respectively. 


Previous to slicing, the silicon ingots must be ground to a uniform 

diameter and this cost has to be included in the cost of the direct material. 

The add-on cost of grinding, listed in line 16 of Table H-3, consists of 

two parts: a) the cost of the grinding operation which is projected to be 

$0.20/cm-crystal length, based on industry data^ ; and b) the cost of the 

1 1 -2 18 

silicon lost from grinding, which is equal to (Tq ' q ' ^ -)q * (Si ingot 
price ($/kg)), where II-2.18 is the percentage of material lost in grinding 
Summing the add-on grinding cost to the Si ingot price yields the cost of 
ground silicon prices (Table I 1-2, lines 17,-22, 27) which are used to cal- 
culate silicon wafer prices. 


Also of interest in our analysis is the cost cf the silicon lost in 
kerf and broken wafers. These values, shown in line 18, 23, 28 of Table H-3, 
are the product of the unit area kerf and wafer loss mass (Table H-2, line 
20) and the ground silicon prices. The add-on wafer prices, shown in lines 
19, 24 and 29 of Table I 1-3 are defined, here, as the sum of add-on wafer 
price assuming a zero silicon price (Table 1 1 -3 , line 15) and the cost of 
the lost silicon. 


To arrive at a unit area wafer price listed in lines 20, 25, and 3C 
of Table III, the add-on price and the cost of silicon incorporated in the 


wafer are summed. The latter value is the cost of the ground silicon 
ingot multiplied by unit area silicon mass contained in the wafers (Table 
II-?, line 19). 
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4. Cost Structures of the Slicing Processes 


The more important unit area manufacturing cost components for 
selected current production or experimental slicing capabilities, using 
1978 silicon prices, and projected future capabilities, using 1982 and 
1986 projected silicon prices are surrmarized in Table II-4. These silicon 
prices apply to single crystal ingots ground to a uniform diameter. The 
diameter tolerance for the ground ingots has been given as + 0.125mn 
standard by Si 1 tec, and as + 0-075itm by Spectrolab. Also included in this 
table are the costs of the lost silicon and that contained in the wafer. 

In Table II-4, one can observe the decreases in expendible tooling, 
indirect materials, labor and capital costs that are expected for 1982 
in ID multi-blade and slurry multiwire slicing. Illustrated in Figure II-4 
are the more relevant data of Tables II-3 and II-4, in a bar graph format. 
In Figure II-4, the relative impacts of the material, labor and capital 
costs can be readily compared to each other for the current multiblade 
and ID slicing processes and for the near future (1982) projected multi- 
blade, ID, and slurry r ulti'-ire processes. 

As evidenced in Table II-4, the indirect material costs (primarily) 
slurry) and the costs for expendible tooling (the steel blades or wires) 
are much higher for the slurry sawing processes (multibladc and Yasunaga 
multiwire) than those for the fixed abrasive approaches (ID saw and FAST 
wire saw). This is a consequence of the more effective utilization of the 
abrasive in the fixed abrasive system, coupled with longer tool life. Re- 
ductions of these expendible tooling costs for the multiblade and slurry 

multiwire slicing processes are expected in the future through lower cost 

(911) (Ql 

tool fabrication technioues ’ ' and through improved lifetimes v . The 

lower tool cost fabrication techniques are expected to result from larger 
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scale, automated assembly^ and a simplification of the assembly process^). 
Investigations are currently being conducted into possibilities for lowering the 
slurry costs, for instance by reclyclina the slurry or substituting a 
cheaper vehicle (e.g. mineral oil) for the PC oil. In spite of these 
projected reductions, the indirect material and expendible tooling costs 
for the multiblade and the Yasunaga multiwire techniques remain sizable 
components of the total add-on costs for those processes. In the near- 
term projections, these components are 44% and 73% of the add-on cost for 
the multiblade and slurry multiwire processes, respectively. This compares 
to 20% and 9 % of the 1982 projections for the add-on costs in the ID and 
fixed abrasive multiwire saws, respectively. 

The current prices are essentially equal for production wafers cut 
by either the Varian multiblade or the ID sawing processes, although the ID 
saw has twice the productivity (Table II-l, line 24) and experiences lower 
indirect material and tooling costs. The higher productivity directly 
results in lower labor, capital, and return-on-equity costs, as shown in 
Figure II.4.- There Tower processing costs for the ID slicing are counter- 
balanced, however, by a higher silicon consumption resulting from the 
practice to cut the wafers to greater thickness with higher kerf than 
achieved with the slurry saws. At the current silicon prices, this has 
a considerable cost impact. 

The 1978 wafer prices shown here are somewhat lower than the con- 
temporary commerical wafer and the 1978 values of the LSA Interim Price 

(14) 

Allocation Guidelines . This difference results from two facts: a) the 

data of this report do not include the cleaning, etchinq, or polishing 
process steps usually included in commercially sold wafers; and b) the 
standardized indirect cost model (SAMICS-IPEG) purposely omits several 
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indirect charges on partially processed items such as wafers. Since the 
indirect cost structure models a vertically integrated industry, marketing 
costs for wafers, c.cj. are not incurred. 

5. CONCLUSIONS 

The cost-analysis data, and particularly the projections, which in- 
clude reduced expendible tooling and indirect material cost components, shov/ 
that the dominant influence on the add-on price of sliced wafers is the 
productivity of the slicing machine. The machine productivity (the time 
rate of output unit expressed in wafer area) has a direct, inversely propor- 
tional impact on the capital cost allocation to the wafer area produced of 
the cost components for equipment and facility, and on that part of the labor 
expenditures which are devoted to machine monitoring and maintenance, as 
shown in Figure II. 5. Figure II. 5 shows that the effective linear cutting rate 
(the workpiece diameter divided by the slicing time-segment) is 0.55 + 0.3 
cm/h for the multiblade and multiwire processes. The inner diameter diamond- 
coated blade process has an effective linear cutting rate of approximately 
300 cm/h, a nearly C50 times larger value than that for the other processes. 

To achieve comparable machine productivities, the low linear cutting rates 
have to be compensated by simultaneous multiple slicing. The current efforts 
of Crystal Systems, Solarex, and Varian are therefore directed at increasing 
the number of wafers sliced during a run. Current multiblade packages con- 
tain about 250 blades. Varian has built an experimental slicer incorporating 
a blade pack of over 900 blades. Similarly, the wire package proposed by 
Crystal Systems^ ^ is projected to have 750 cutting wires. Solarex hopes to 
slice^^ 333 wafers at a time with the Yasunaga YQ-100 slicing machine. 

The slicing technology improvements projected for the 1982 produc- 
tion lines are based on the results of recent experimental runs and on 
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developments in progress (Table 11-4, Fig. 1 1. 5). For the multi blade saw, 
the primary advancement will be a nearly four-fold productivity increase 
via Varian's development of a machine using a 900-blade-pack. Simul- 
taneously, a 2S% blade thickness reduction In combination with a 

37.5% wafer thickness decrease, while maintaining a wafer yield of 95%, 

2 

is projected to result in an area yield of 0.9 m / kg-Si crystal, a 50% 
increase from Spectrolab's mass to area conversion ratio in slicing 
round wafers. 

Slice and kerf thickness reductions to values similar to those 
projected for the multi blade slurry process, are also expected for the 
ID-sawing method. Recently acquired data from STC are reflected in a 
1982 projection for 10-cm diameter crystals using ID slicing with ingot 
rotation, as shown in Table II-4. 

The wafers from this process are expected to be 225 ym thick with 210 ym 
kerf. In addition, crystal rotation is expected to double the effective 
cutting rate of the ID process. This essentially doubles the productivity 
of the ID saw, and results in comparable projected productivities for the 
900-blade multiblade and the ID sawing processes. Remaining differences in 
the costs of these two processes are, however, overshadowed by the cost 
of the silicon incorporated into the wafer or lost. At the projected 
1982 price for ground single crystal ingots, the cost of this silicon 
still amounts to nearly 80% of the wafer price. 

One slicing method has been projected to 1986, primarily, because 
only a comparatively small base of experimental data is available, so 
that this method cannot be expected to be in significant commercial 
operation by 1982. This method is Crystal Systems' fixed abrasive 
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multiwire sawing. The current projections are contained in Table II-4, 
while Table II-3B is based on earlier inputs. The difference results 
primarily from a recently communicated reduction in tooling costs based 
on wlrehead fabrication Improvements, and from the use in Table IV of 
a more conservative effective cutting rate corresponding to the experi- 
mentally found rates averaged over the life of the bladehead. The pro- 
cess add-on costs are comparable to those of the two previously discussed 

processes. If the silicon price of 1982 would have been used, an approx- 
2 

imately $1 1/m lower wafer price would have resulted in comparison to the 
ID process. While the fixed abrasive multiwire process currently projects 
the lowest wafer price, it is also the one with the least experience 
data. It is therefore of great importance to gain a significant data 
base through pilot line operation. 

Considering the uncertainties in the projections, the data indi- 
cate no considerable differences in the competitiveness of the three ap- 
proaches, and a reasonable potential for all three to meet the 1986 
guideline goal . 
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III. Junction Formation 


1. Introduction 

In general, the first major step in the process 
sequence from quartzite to complete solar modules, follow- 
ing the generation of the silicon wafers, ribbons, or 
sheet, is that of pn junction formation. Of the primary 
present junction forming processes, gaseous diffusion was 
examined in great detail as the base case. Then, other 
diffusion processes and ion implantation, proposed as methods 
for lower cost junction formation, were analyzed. 

As with our crystal growing and slicing studies, the 
evaluations were started with the current methods of diffu- 
sion and ion implantation for which a large amount of the 
needed information is available. 

For the diffusion process, we have tabulated production 
experience data from Spectroiab^ and projections made by 

Motorola ^ and RCA. ^ In our studies of ion implantation 

( 4 ) 

of pn junctions, experimental data from Spire using a 

modified Varian-Extrion machine, along with material, labor, 

and capital projections made by Lockheed^, Motorola, ^ 

( 4 ) 

and Spire for their proposed machines were examined. 
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2. Diffusion Processes 


The principles of the diffusion processes are ex- 
tensively described in the literature and will not be repeated 
here. These processes are the ones most widely used in the 
semiconductor industry for the formation of pn junctions. 

There are two basic variations: the infinite dopant source 

method, which results in a complementary error function 
distribution of the dopant with depth below the surface, and 
the finite source method, which results in a Gaussian distri- 
bution. The infinite source method represents principally a 
one step process in which the dopant is transferred from an 
impurity carrier material into a surface layer of the silicon 
wafer which becomes dopant-rich up to the solid solubility 
of the respective impurity in silicon (source layer) , and 
from which the diffusion into the wafer takes place, 
simultaneous with the formation of this source layer. In 
this case, the dopant carrier material (the "infinite source") 
remains present throughout the diffusion process. In the 
finite source method, the dopant carrier material is removed 
after a "back-on" cycle for source layer formation, and 
diffusion continues subsequently for a longer time period 
in the "drive-in" cycle, redistributing the dopant 
originally contained in the source layer deeper into the wafer. 

In solar cell production, the infinite source method is 
used exclusively, and the processes vary only in the dopant 
carrier material used and the form of its application. 

Very commonly used is phosphine gas in an "open rube" system, 
generating a phosphorus glass layer on the silicon wafer. 
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This glass layer is formed during the first part of the 
diffusion step, but in the same furnace and at the same 
temperature as the rest of the process. This phosphorus 
glass is then the dopant carrier material from which the 
source layer is formed. The phosphorus glass is removed only 
after completion of the diffusion process. In an alternate 

approach , the dopant carrier material is supplied in semi-liquid 
form and applied to the wafer in a spray-on or spin-on step, 
with subsequent glass formation and diffusion occuring in a 
belt furnace potentially as a continuous flow process. While, 
in the "gaseous" diffusion process, the gas cycling and flow 
control is performed fully automatically, it is still at 
batch process with usual furnace loading and unloading. 

As diffusion is currently the major competitive process, 
we have examined the attributes and costs of present and pro- 
jected future diffusion processes. In the current production 
operation, Spectrolab uses open-tube diffusion with phosphine 
diluted heavily in hydrogen to form a pn junction. Thanks 
to the data supplied generously by R. Oliver and E.L. Ralph 
of Spectrolab, ^ we have been able to make a detailed analysis 
of the present diffusion process as a baseline case. The 
detail data of this process are presented in Tables III-4A 
through 4C and the cost summary in Table III-3. The diffusion 
process takes approximately 35 minutes for a run containing 
75 wafers of 7.62-cm diameter. We have observed that the 
process as performed by Spectrolab is very labor intensive. 

The reason is that only two diffusion furnaces are needed 
to nandle the entire production, but one operator is needed 
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to attend the process. Thus, this operator devotes most 
of his/her time to manually loading and unloading wafers 
onto and from the quartz diffusion boat, which could be 
done mechanically. If one assumes automatic wafer feeding, 
the operator's time could easily be reduced to 10 minutes 
per run, and the processing add-on price would be re- 
duced to approximately $9. 50/m from the present value 
of $12.74/m^ (SAMICS methodology). 

Another significant cost contributor, and one that 
has oeen ignored in most projections for future diffu- 
sion processes, is that for cleaning the quartz furnace 
tubes and boats, which is usually done with a HF-HNO^-f^O 
solution, as often as twice a day. Frequent quartzware 

cleaning has been found instrumental to maintaining high 

2 

cell efficiency, but it contributes $2. 23/m to the 
diffusion add-on price in the Spectrolab process. This 
price contribution was calculated assuming that the 
quartz cleaning operation requires 1 h/work day of labor, and 
a tube cleaning tower which costs $15,000 including installa- 
tion, and which is shared between the two furnaces. About half of 
this cleaning cost contribution is due to equipment costs, 
with the remainder, listed in decreasing magnitude, shared 
between labor, facility, and material costs. 

Future diffusion price projections, such as for 

Motorola's phosphine (PH^) process , ^ also detailed in 

Table III-4A to 4C, are about a factor of four lower 

2 

than present calculated prices ($3. 10/m compared to 
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$12. 74/m 2 ). The Motorola process which has approximately 
the same wafer throughput rate as Spectrolab's current 
process, is applied to 12-cm diameter wafers, rather than 

7.62-cm wafers in the Spectrolab process. 

The 12-cm wafers have an area that is nearly 2.5 
times larger than that of the 7.62-cm wafer, accounting 
for most of the cost difference between Motorola's and 
Spectrolab's diffusion processes. The rest of the 
cost difference can be attributed to the more automated 
nature of the Motorola process, requiring half-as-much 
labor as Spectrolab requires, and the lack of inclusion, 
by Motorola, °f costs for cleaning the quartzware. On the 
other hand, notable are Motorola's projected use of 
significantly more energy and direct material (pnosphine) 
than Spectrolab is consuming now. 

Currently, the PN junction formation process by 
diffusion is not a large cost-contributing factor in cell 
processing. In application of the diffusion process, a 
separate annealing step is not required, at least not beyond 
a somewhat slowed cooling rate from diffusion temperature. 

A separate post-annealing step is, however, required after 
the ion implantation process to reduce the crystal damage 
resulting from implantation, and to activate the impurity 
species. Therefore, the annealing cost must be included in 
any cost analysis of ion implantation. Using a Thermco 
eight-tube diffusion furnace, which has an output rate of 

2 

1,000 12-cm diameter wafers/h, an add-on price of $1.18/m 
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was calculated for the annealinq process steo. 

If ion implantation is to replace diffusion, it may be 
able to become cost competitive only as part of a more 
extended sequence of vacuum processes, or by producing 
cells of significantly higher performance than achievable 
by the diffusion process. 


3 • Principals and Application of Ion Implantation 


Ion implantation is a method for introducing dopant 
material below the semiconductor surface to form PN or 
high/low junctions. In the common type of ion implantation 
machine, the source material , usually a cnemical compound 
containing the dopant, is broken down and ionized under 
electron bombardment in the ionization chamber, the ions 
are extracted from this chamber by an electric field and 
further accelerated and collimated, purified using a mass 
spectrometer, and then scanned either electrically, 
magnetically , or mechanically while impinging on the semi- 
conductor wafer to be implanted. The top portion of 
Figure III.l shows a schematic presentation of such an ion 
implanter with a magnetic analyzer. In simpler machines, 
as shown in the bottom part of Figure III.l, functions such 
as beam collimation, mass analysis, or scanning may be omitted. 

In the machine shown in the top part of Figure III.l, the 
source material can be ionized in a number of ways, 

the principal ones of which are: heating and electron 
bombardment of the source material from a high temperature 
emitter, called the"hot cathode source"; electron dis- 
charge from a low work function emitter, such as barium, 
under the influence of a strong electric field, into 
the vaporized source material to form a plasma (cold 
cathode source); or by microwave discharge. In any of 
the mentioned sources, a magnetic field can be applied 
to concentrate the plasma density and increase the 
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efficiency of ionization. This will also result, 
though, in lower source lifetime and a larger energy 
spread of the ions. 

Three principal types of hot cathode ion sources 
are used in the implanters mentioned in this report. 

In all, the current density from a metal surface at 
temperature T with a work function of $ is principally 
described by: 

2 -e 4>/kT TTT 

n * AT e III- 

J e 


However, at adequately high emission rates, the current 
density j is usually reduced below the value given by 

Eq. III-2.1 because of space-charge effects, in which the 
mutual repulsion of the electrons crowding the space near 

the filament inhibits further emission. The electron density 

then becomes: 

v 2/ 3 

j = \ 1/2 III-2.2 

9Ttd (m/2e) * 


where V is the voltage between the cathode and anode, d is the 
thickness of the electron sheath and m/e is the electron's mass 
to charge ratio. The production of positive ions in the source 
chamber tends to neutralize this "electron cloud" and reduce 
the space charge effects. The cathode current thus increases 
in the presence of positive ions. 

In the "Freeman source", the heated wire cathode has its 
terminals on opposite sides of the "extraction gap" through 
which the ions leave. In the "Chavet source", the filament wire 
is looped so that its electrodes are on the same side of 
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the extraction gap. The thavet filament configuration was 
designed to increase the filament's lifetime by decreasing 
its exposure to the back-streaming ions and thus reduce the 
sputtering caused by them. Another thermionic source is the 

hollow cathode in which the interior of a cylindrical 
cavity is coated with a low work function material, such as 
barium oxide. Upon introduction of the vaporized source 
material, an arc discharge takes place between the cathode 
and anode so that the source material is ionized. As a result 
of applied high voltage, the ions are extracted through 
a hole in the cathode. Vaporized atoms also pass through this 
aperture. They are subsequently ionized by the accelerated 
electrons. One configuration of a cold cathode source known 
as the "Penning source", has an anode that is also cylindri- 
cal in shape with the end plates forming the cathode. In 
addition, a magnetic field is applied parallel to the cylindri- 
cal axis of the "Penning source" to force electrons from the 
cathode to form helical trajectories, thus increasing their 
path length and enhancing the ionization efficiency. 

After the ion beam is extracted from the source chamber, 
it is accelerated through a potential drop. For small 
acceleration energies (< 30 keV) , a single gap electrode 
could be used. The accelerated ion beam is then subjected 
to a magnetic field for mass separation. A singly charged 
ion of atomic mass M (AMU's) moving through a magnetic field 
with strength B (in gauss) will be deflected into a circular 
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path with the radius of curvature equal to 

„ 143.95 ...... 1/2 ^ , 

R » — g (MV) ' cm, III-2.3 

where V is the acceleration voltage. The dispersion between 
ions of two different masses is 

D M - ^ R cm . m-2.4 

M M 

In order to achieve good mass resolution, power supplies 
to the acceleration and magnet regions must have stabilities 
of 1 part in 10,000. 

To form the junction, the analyzed beam is then 
scanned, with one of the techniques mentioned previously, 
on the silicon substrate. Overscanning is necessary be- 
cause of the tails in the Gaussion distribution of the 
ion concentration in the beam. 

Junction formation using ion implantation offers 
several potential advantages over the diffusion process. It is 
a dry, vacuum process, thus avoiding potential con- 
tamination from impurities contained in spin-on or gaseous 
vehicles for the dopants used in some varieties of tne 
diffusion process. Where selective introduction of the 
dopant is wanted, this may be accomplished without application 
of masking and subsequent stripping, and without back-surface 
etching because of double-sided impurity penetrations. 

Thus, ion implantation can involve fewer handling or 
transferring operations than the diffusion process, and 
consequently can result in labor savings and increased 

Iii-ll 
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thickness of the electron sheath and m/e is the electron's mass 
to charge ratio. The production of positive ions in the source 
chamber tends to neutralize this "electron cloud" and reduce 
the space charge effects. The cathode current thus increases 
in the presence of positive ions. 

In the "Freeman source", the heated wire cathode has its 
terminals on opposite sides of the "extraction gap" through 
which the ions leave. In the "Chavet source", the filament wire 
is looped so that its electrodes are on the same side of 
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the extraction gap. The thavet filament configuration was 
designed to increase the filament's lifetime by decreasing 
its exposure to the back-streaming ions and thus reduce the 
sputtering caused by them. Another thermionic source is the 

hollow cathode in which the interior of a cylindrical 
cavity is coated with a low work function material, such as 
barium oxide. Upon introduction of the vaporized source 
material, an arc discharge takes place between the cathode 
and anode so that the source material is ionized. As a result 
of applied high voltage, the ions are extracted through 
a hole in the cathode. Vaporized atoms also pass through this 
aperture. They are subsequently ionized by the accelerated 
electrons. One configuration of a cold cathode source known 
as the "Penning source", has an anode that is also cylindri- 
cal in shape with the end plates forming the cathode. In 
addition, a magnetic field is applied parallel to the cylindri- 
cal axis of the "Penning source" to force electrons from the 
cathode to form helical trajectories, thus increasing their 
path length and enhancing the ionization efficiency. 

After the ion beam is extracted from the source chamber, 
it is accelerated through a potential drop. For small 
acceleration energies (< 30 keV) , a single gap electrode 
could be used. The accelerated ion beam is then subjected 
to a magnetic field for mass separation. A singly charged 
ion of atomic mass M (AMU's) moving through a magnetic field 
with strength B (in gauss) will be deflected into a circular 
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path with the radius of curvature equal to 

„ 143.95 ...... 1/2 ^ , 

R » — g (MV) ' cm, III-2.3 

where V is the acceleration voltage. The dispersion between 
ions of two different masses is 

D M - ^ R cm . m-2.4 

M M 

In order to achieve good mass resolution, power supplies 
to the acceleration and magnet regions must have stabilities 
of 1 part in 10,000. 

To form the junction, the analyzed beam is then 
scanned, with one of the techniques mentioned previously, 
on the silicon substrate. Overscanning is necessary be- 
cause of the tails in the Gaussion distribution of the 
ion concentration in the beam. 

Junction formation using ion implantation offers 
several potential advantages over the diffusion process. It is 
a dry, vacuum process, thus avoiding potential con- 
tamination from impurities contained in spin-on or gaseous 
vehicles for the dopants used in some varieties of tne 
diffusion process. Where selective introduction of the 
dopant is wanted, this may be accomplished without application 
of masking and subsequent stripping, and without back-surface 
etching because of double-sided impurity penetrations. 

Thus, ion implantation can involve fewer handling or 
transferring operations than the diffusion process, and 
consequently can result in labor savings and increased 
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yields. However, ion . .^plantation requires an annealing 

step, which will be further discussed later on. It has 
( 8 ) 

been suggested to use ion implantation as an integral 
part of a total vacuum process sequence for fabricating 
solar cells after wafer or sheet generation. Such a sequence, 
although high in capital costs, could result in labor savings 
and high yields. 

The charge on the dopant ions allows for mass-spectro- 
scopic separation using magnetic fields, and for accurate 
measurement of the ion flux entering the deposition region, 
as lonq as nejtral and doubly charged particles are handled 
correctly. The ion beam currents can be readily measured 
by placing a Faraday cup in the beam's path, but this requires 
a preceding calibration to determine the fraction of uncharged 
and doubly charged ions. The mass analysis and ion current 
measurement features of the ion implantation process can pro- 
vide better control over the quantity and quality of the 
dopant than other processes, and can therefore be applied 

to obtain better process uniformity and repeatability. Dose 

(9) 

uniformities of + 5% (2o) are achievable 

Since ion implantation can be performed at or near 

room temperature, low energy implantations can result in 

o 

original dopant penetration of less than 100A. This is shallower 

than can be achieved in most high temperature source deposition 

steps in the diffusion process. 

Upon entering the substrate, the dominant interactions 

of the ion are wi th the electrons of the lattice, which 

slow the ion down through kinetic energy transfer. After 
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slow the ion down through kinetic energy transfer. After 
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this initial slow-down to sufficiently low energies, i.e., 

2 

ion velocity less than Z^ /fi, collisions of the ion take 
place with the nuclei which completely stop the ion. In 
most cases, the stopped ion rests interstially in the 
crystal lattice. The largest impurity concentration is 
thus found at a penetration distance, "x ", from the 
surface. As a first approximation in the region where 
nuclear collisions dominate, the penetration depth is 
proportional to the square root of the ion beam energy. 
This penetration depth is described by: 


0.7 (Zj 273 + Z si 2/3 ) 


z z . 

I Si 


*1 + M si 


M, 


E i (X). III-2.5 


Ej is the energy of the ion beam in eV, Z and M refer to 
the atomic number and atomic weight, respectively, while 
the subscripts I and Si refer to the ion and to Si, re- 
spectively. The concentration varies from the penetration 
distance approximately according to a Gaussian distribution 
and the impurity distribution can be described by the 
empirical relationship. 


C (x) - C p exp (- (x-Xp) 2 /2 a R 2 ) . HI-2. 
C is the concentration at the penetration distance, and 
o- is called the standard deviation of the concentration 
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function, or the distance from x a t which the concentra- 
tion is equal to C p //e7 The peak concentration depends 
upon the ion beam current (i) and the implantation time (t) 
or 

10 16 i • t 

C D cm III-2.7 

' 4 °R 

The unit of i is nft, that of t is seconds, and o is 
is given in nm. 

The penetration distance can be calculated from 

electron and nuclear ionic collisions only if "channeling" 
does not occur. Channeling is the name given to the con- 
siderably enhanced penetration distance of ions which are 
aligned with low index crystallographic directions, and 
therefore travel parallel to and in between high atomic density 
crystal planes. Since ions travelling this path experience 
relatively fewer collisions with silicon atoms, they can travel 
further into the silicon. To avoid channclinn, the beam must 
be oriented at a slight angle (‘7°) from the orientation of 
the low index crystallograph ' c arcs. This increases the 
apparent distribution of atoms in the crystal plane 
normal to the ion beam's path, and thus increases the pro- 
bability of ion-nuclei collisions. 

The implantation process results in the displacement of 
the silicon atoms from their normal lattice sites by the ion 
collisions, thus creating "vacancies" and "interstitials". 

The implanted impurity atoms, predominately located at intersti- 
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tial sites, are not electrically active. Thus few impurity 
atoms which take up substitutional positions, tend to 
compensate the originally present impurity atoms of oppo- 
site dopant type, and to shift the Fermi level of the silicon 
towards the center of the energy gap. Annealing of the 
ion implanted wafers is required be . 1 ^ to reduce the mentioned 
crystal structure damage resulting from the implantation 
process, some of which is electrically active (recombination 
and trapping centers) , and to electrically "activate" the 
dopant impurity by moving its implanted atoms from inter- 
stitial to substitutional sites. This annealing is usually 
accomplished by a high temperature soak, called thermal 
annealing . 

Thermal annealing broadens the impurity profile, usually 
to a junction depth as great or greater than obtained by use 
of relatively low temperature, short tiir^ diffusions as they 
are normally used for solar cell production. Nevertheless, 
ion implantation followed by thermal annealing is capable of 
producing solar cells with efficiencies equivalent to those 
prepared using diffusion. A thermal annealing cycle of lh at 

450°C and 0.5h at 859°C has been repeatedly found to yield per- 
formance-wise competitive silicon solar cells. This, in part 

negates the potential advantage of being able to control the 
dopant profile at will by varying the implantation energy and 
dosage. Such "designed profiles" might lead to higher effic- 
iency solar cells than obtained so far. Electron and laser beam 
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annealing, as short transient annealing methods, have therefore 
been and are being investigated because of the potential 
profile maintenance as well as several other anticipated 
advantages. These advantages have, however, so far not 
been realized, and cells with efficiencies comparable to 
those obtained by the oven annealing process have so far not 
been reported. However, electron beam annealing, followed by 
a low temperature soak, at approximately 500°C, is now said 
to produce cells with efficiencies comparable to those using 
the "thermal annealing" process. The necessity of the low 
temperature soak seems to indicate that the electron beam 
annealing process does not reduce crystal damage or permit 
gettering as well as thermal annealing. Although the 

electron beam pulse anneal, followed by a low temperature soak, 
will be a more costly process than a higher temperature 
activation/annealing soak one, it appears attractive because 
of the greater freedom in selection and control of the impurity 
profile. If this pulse anncal/soak process, or some other, 
simpler process for activation/anneal i ng could be developed, 
so that ion implanted solar cells might attain higher 
efficiencies than cells prepared by diffusion processes, then 
ion implantation would become a most interesting process 
option, even at a possibly somewhat higher process cost than 
diffusion . 

With the attainability of a potential efficiency advantage 
of the ion implanted solar cells over the diffusion produced 
cells not demonstrated, the usefulness of ion implantation, 
as part of an LSA solar cell sequence will be determined by 
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the potentially achievable cost reduction. Currently, the 
high capital costs, the low reliability, and the low through- 
put rate of ion implantation machines, make junction formation 
with them too costly to be used for large scale solar cell 
production. Large cost reductions are, however, expected to 
be accomplished in the future (1986) by several approaches. 
Approaches to this end include the introduction of large 
throughput machines with high current, hot cathode ion beam 
sources incorporating an analyzer and more automated operation 

*7 \ 

through computer control' , and the development of ion im- 
planters with unanalyzed or roughly analyzed ion beams 
using hollow cathode sources. Some current and future 
applications of ion implantation are listed in Table I along 
with the conditions contingent to the two potential advantages 
of lower cost and higher efficiency. 

In consequence of this discussion, it has to be observed, 
particularly in reading the following sections, that ion 
implantation and activation/annealing are inseparably con- 
nected, and that, even though only one of the process steps 
may be mentioned, it cannot be usefully carried out without 
the other. But there are several ways of carrying out ion 
implantation as well as annealing, with various combinations 
of the two, so that one really deals with an ion implantation/ 
annealing submatrix of options. Consequently, both the pro- 
cess cost and the performance of the finished solar cells 
depend on the option selection within this submatrix. 
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efficiency of ionization. This will also result, 
though, in lower source lifetime and a larger energy 
spread of the ions. 


Three principal types of hot cathode ion sources 
are used in the implanters mentioned in this report. 

In all, the current density from a metal surface at 
temperature T with a work function of $ is principally 
described by: 


. _2 -e 4>/kT 
n * AT e 
J e 


III-2.1 


However, at adequately high emission rates, the current 
density j is usually reduced below the value given by 

Eq. III-2.1 because of space-charge effects, in which the 
mutual repulsion of the electrons crowding the space near 

the filament inhibits further emission. The electron density 

then becomes: 


le = 


V 


2/3 


9Ttd 2 (m/2e) 1/2 


III-2.2 


where V is the voltage between the cathode and anode, d is the 
thickness of the electron sheath and m/e is the electron's mass 
to charge ratio. The production of positive ions in the source 
chamber tends to neutralize this "electron cloud" and reduce 
the space charge effects. The cathode current thus increases 
in the presence of positive ions. 

In the "Freeman source", the heated wire cathode has its 
terminals on opposite sides of the "extraction gap" through 
which the ions leave. In the "Chavet source", the filament wire 

is looped so that its electrodes are on the same side of 
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the space charge effects. The cathode current thus increases 
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terminals on opposite sides of the "extraction gap" through 
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is looped so that its electrodes are on the same side of 
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the extraction gap. The thavet filament configuration was 
designed to increase the filament's lifetime by decreasing 
its exposure to the back-streaming ions and thus reduce the 
sputtering caused by them. Another thermionic source is the 

hollow cathode in which the interior of a cylindrical 
cavity is coated with a low work function material, such as 
barium oxide. Upon introduction of the vaporized source 
material, an arc discharge takes place between the cathode 
and anode so that the source material is ionized. As a result 
of applied high voltage, the ions are extracted through 
a hole in the cathode. Vaporized atoms also pass through this 
aperture. They are subsequently ionized by the accelerated 
electrons. One configuration of a cold cathode source known 
as the "Penning source", has an anode that is also cylindri- 
cal in shape with the end plates forming the cathode. In 
addition, a magnetic field is applied parallel to the cylindri- 
cal axis of the "Penning source" to force electrons from the 
cathode to form helical trajectories, thus increasing their 
path length and enhancing the ionization efficiency. 

After the ion beam is extracted from the source chamber, 
it is accelerated through a potential drop. For small 
acceleration energies (< 30 keV) , a single gap electrode 
could be used. The accelerated ion beam is then subjected 
to a magnetic field for mass separation. A singly charged 
ion of atomic mass M (AMU's) moving through a magnetic field 
with strength B (in gauss) will be deflected into a circular 


ITI-10 


path with the radius of curvature equal to 

„ 143.95 ...... 1/2 ^ , 

R » — g (MV) ' cm, III-2.3 

where V is the acceleration voltage. The dispersion between 
ions of two different masses is 

D M - ^ R cm . m-2.4 

M M 

In order to achieve good mass resolution, power supplies 
to the acceleration and magnet regions must have stabilities 
of 1 part in 10,000. 

To form the junction, the analyzed beam is then 
scanned, with one of the techniques mentioned previously, 
on the silicon substrate. Overscanning is necessary be- 
cause of the tails in the Gaussion distribution of the 
ion concentration in the beam. 

Junction formation using ion implantation offers 
several potential advantages over the diffusion process. It is 
a dry, vacuum process, thus avoiding potential con- 
tamination from impurities contained in spin-on or gaseous 
vehicles for the dopants used in some varieties of tne 
diffusion process. Where selective introduction of the 
dopant is wanted, this may be accomplished without application 
of masking and subsequent stripping, and without back-surface 
etching because of double-sided impurity penetrations. 

Thus, ion implantation can involve fewer handling or 
transferring operations than the diffusion process, and 
consequently can result in labor savings and increased 
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yields. However, ion . .^plantation requires an annealing 

step, which will be further discussed later on. It has 
( 8 ) 

been suggested to use ion implantation as an integral 
part of a total vacuum process sequence for fabricating 
solar cells after wafer or sheet generation. Such a sequence, 
although high in capital costs, could result in labor savings 
and high yields. 

The charge on the dopant ions allows for mass-spectro- 
scopic separation using magnetic fields, and for accurate 
measurement of the ion flux entering the deposition region, 
as lonq as nejtral and doubly charged particles are handled 
correctly. The ion beam currents can be readily measured 
by placing a Faraday cup in the beam's path, but this requires 
a preceding calibration to determine the fraction of uncharged 
and doubly charged ions. The mass analysis and ion current 
measurement features of the ion implantation process can pro- 
vide better control over the quantity and quality of the 
dopant than other processes, and can therefore be applied 

to obtain better process uniformity and repeatability. Dose 

(9) 

uniformities of + 5% (2o) are achievable 

Since ion implantation can be performed at or near 

room temperature, low energy implantations can result in 

o 

original dopant penetration of less than 100A. This is shallower 

than can be achieved in most high temperature source deposition 

steps in the diffusion process. 

Upon entering the substrate, the dominant interactions 

of the ion are wi th the electrons of the lattice, which 

slow the ion down through kinetic energy transfer. After 
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this initial slow-down to sufficiently low energies, i.e., 

2 

ion velocity less than Z^ /fi, collisions of the ion take 
place with the nuclei which completely stop the ion. In 
most cases, the stopped ion rests interstially in the 
crystal lattice. The largest impurity concentration is 
thus found at a penetration distance, "x ", from the 
surface. As a first approximation in the region where 
nuclear collisions dominate, the penetration depth is 
proportional to the square root of the ion beam energy. 
This penetration depth is described by: 


0.7 (Zj 273 + Z si 2/3 ) 


z z . 

I Si 


*1 + M si 


M, 


E i (X). III-2.5 


Ej is the energy of the ion beam in eV, Z and M refer to 
the atomic number and atomic weight, respectively, while 
the subscripts I and Si refer to the ion and to Si, re- 
spectively. The concentration varies from the penetration 
distance approximately according to a Gaussian distribution 
and the impurity distribution can be described by the 
empirical relationship. 


C (x) - Cp exp (- (x-Xp) 2 /2 a R 2 ) . HI-2. 
C is the concentration at the penetration distance, and 
o- is called the standard deviation of the concentration 
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function, or the distance from x a t which the concentra- 
tion is equal to C p //e7 The peak concentration depends 
upon the ion beam current (i) and the implantation time (t) 
or 

10 16 i • t 

C D cm III-2.7 

' 4 °R 

The unit of i is n&, that of t is seconds, and o is 
is given in nm. 

The penetration distance can be calculated from 

electron and nuclear ionic collisions only if "channeling" 
does not occur. Channeling is the name given to the con- 
siderably enhanced penetration distance of ions which are 
aligned with low index crystallographic directions, and 
therefore travel parallel to and in between high atomic density 
crystal planes. Since ions travelling this path experience 
relatively fewer collisions with silicon atoms, they can travel 
further into the silicon. To avoid channclinn, the beam must 
be oriented at a slight angle (‘7°) from the orientation of 
the low index crystallograph ' c arcs. This increases the 
apparent distribution of atoms in the crystal plane 
normal to the ion beam's path, and thus increases the pro- 
bability of ion-nuclei collisions. 

The implantation process results in the displacement of 
the silicon atoms from their normal lattice sites by the ion 
collisions, thus creating "vacancies" and "interstitials". 

The implanted impurity atoms, predominately located at intersti- 
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tial sites, are not electrically active. Thus few impurity 
atoms which take up substitutional positions, tend to 
compensate the originally present impurity atoms of oppo- 
site dopant type, and to shift the Fermi level of the silicon 
towards the center of the energy gap. Annealing of the 
ion implanted wafers is required be . 1 ^ to reduce the mentioned 
crystal structure damage resulting from the implantation 
process, some of which is electrically active (recombination 
and trapping centers) , and to electrically "activate" the 
dopant impurity by moving its implanted atoms from inter- 
stitial to substitutional sites. This annealing is usually 
accomplished by a high temperature soak, called thermal 
annealing . 

Thermal annealing broadens the impurity profile, usually 
to a junction depth as great or greater than obtained by use 
of relatively low temperature, short tiir^ diffusions as they 
are normally used for solar cell production. Nevertheless, 
ion implantation followed by thermal annealing is capable of 
producing solar cells with efficiencies equivalent to those 
prepared using diffusion. A thermal annealing cycle of lh at 

450°C and 0.5h at 859°C has been repeatedly found to yield per- 
formance-wise competitive silicon solar cells. This, in part 

negates the potential advantage of being able to control the 
dopant profile at will by varying the implantation energy and 
dosage. Such "designed profiles" might lead to higher effic- 
iency solar cells than obtained so far. Electron and laser beam 
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annealing, as short transient annealing methods, have therefore 
been and are being investigated because of the potential 
profile maintenance as well as several other anticipated 
advantages. These advantages have, however, so far not 
been realized, and cells with efficiencies comparable to 
those obtained by the oven annealing process have so far not 
been reported. However, electron beam annealing, followed by 
a low temperature soak, at approximately 500°C, is now said 
to produce cells with efficiencies comparable to those using 
the "thermal annealing" process. The necessity of the low 
temperature soak seems to indicate that the electron beam 
annealing process does not reduce crystal damage or permit 
gettering as well as thermal annealing. Although the 

electron beam pulse anneal, followed by a low temperature soak, 
will be a more costly process than a higher temperature 
activation/annealing soak one, it appears attractive because 
of the greater freedom in selection and control of the impurity 
profile. If this pulse anncal/soak process, or some other, 
simpler process for activation/anneal i ng could be developed, 
so that ion implanted solar cells might attain higher 
efficiencies than cells prepared by diffusion processes, then 
ion implantation would become a most interesting process 
option, even at a possibly somewhat higher process cost than 
diffusion . 

With the attainability of a potential efficiency advantage 
of the ion implanted solar cells over the diffusion produced 
cells not demonstrated, the usefulness of ion implantation, 
as part of an LSA solar cell sequence will be determined by 
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the potentially achievable cost reduction. Currently, the 
high capital costs, the low reliability, and the low through- 
put rate of ion implantation machines, make junction formation 
with them too costly to be used for large scale solar cell 
production. Large cost reductions are, however, expected to 
be accomplished in the future (1986) by several approaches. 
Approaches to this end include the introduction of large 
throughput machines with high current, hot cathode ion beam 
sources incorporating an analyzer and more automated operation 

*7 \ 

through computer control' , and the development of ion im- 
planters with unanalyzed or roughly analyzed ion beams 
using hollow cathode sources. Some current and future 
applications of ion implantation are listed in Table I along 
with the conditions contingent to the two potential advantages 
of lower cost and higher efficiency. 

In consequence of this discussion, it has to be observed, 
particularly in reading the following sections, that ion 
implantation and activation/annealing are inseparably con- 
nected, and that, even though only one of the process steps 
may be mentioned, it cannot be usefully carried out without 
the other. But there are several ways of carrying out ion 
implantation as well as annealing, with various combinations 
of the two, so that one really deals with an ion implantation/ 
annealing submatrix of options. Consequently, both the pro- 
cess cost and the performance of the finished solar cells 
depend on the option selection within this submatrix. 
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efficiency of ionization. This will also result, 
though, in lower source lifetime and a larger energy 
spread of the ions. 

Three principal types of hot cathode ion sources 
are used in the implanters mentioned in this report. 

In all, the current density from a metal surface at 
temperature T with a work function of $ is principally 
described by: 

2 -e 4>/kT TTT 

n * AT e III- 

J e 


However, at adequately high emission rates, the current 
density j is usually reduced below the value given by 

Eq. III-2.1 because of space-charge effects, in which the 
mutual repulsion of the electrons crowding the space near 

the filament inhibits further emission. The electron density 

then becomes: 

v 2/ 3 

j = \ 1/2 III-2.2 

9Ttd (m/2e) * 


where V is the voltage between the cathode and anode, d is the 
thickness of the electron sheath and m/e is the electron's mass 
to charge ratio. The production of positive ions in the source 
chamber tends to neutralize this "electron cloud" and reduce 
the space charge effects. The cathode current thus increases 
in the presence of positive ions. 

In the "Freeman source", the heated wire cathode has its 
terminals on opposite sides of the "extraction gap" through 
which the ions leave. In the "Chavet source", the filament wire 
is looped so that its electrodes are on the same side of 
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the extraction gap. The thavet filament configuration was 
designed to increase the filament's lifetime by decreasing 
its exposure to the back-streaming ions and thus reduce the 
sputtering caused by them. Another thermionic source is the 

hollow cathode in which the interior of a cylindrical 
cavity is coated with a low work function material, such as 
barium oxide. Upon introduction of the vaporized source 
material, an arc discharge takes place between the cathode 
and anode so that the source material is ionized. As a result 
of applied high voltage, the ions are extracted through 
a hole in the cathode. Vaporized atoms also pass through this 
aperture. They are subsequently ionized by the accelerated 
electrons. One configuration of a cold cathode source known 
as the "Penning source", has an anode that is also cylindri- 
cal in shape with the end plates forming the cathode. In 
addition, a magnetic field is applied parallel to the cylindri- 
cal axis of the "Penning source" to force electrons from the 
cathode to form helical trajectories, thus increasing their 
path length and enhancing the ionization efficiency. 

After the ion beam is extracted from the source chamber, 
it is accelerated through a potential drop. For small 
acceleration energies (< 30 keV) , a single gap electrode 
could be used. The accelerated ion beam is then subjected 
to a magnetic field for mass separation. A singly charged 
ion of atomic mass M (AMU's) moving through a magnetic field 
with strength B (in gauss) will be deflected into a circular 
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path with the radius of curvature equal to 

„ 143.95 ...... 1/2 ^ , 

R » — g (MV) ' cm, III-2.3 

where V is the acceleration voltage. The dispersion between 
ions of two different masses is 

D M - ^ R cm . m-2.4 

M M 

In order to achieve good mass resolution, power supplies 
to the acceleration and magnet regions must have stabilities 
of 1 part in 10,000. 

To form the junction, the analyzed beam is then 
scanned, with one of the techniques mentioned previously, 
on the silicon substrate. Overscanning is necessary be- 
cause of the tails in the Gaussion distribution of the 
ion concentration in the beam. 

Junction formation using ion implantation offers 
several potential advantages over the diffusion process. It is 
a dry, vacuum process, thus avoiding potential con- 
tamination from impurities contained in spin-on or gaseous 
vehicles for the dopants used in some varieties of tne 
diffusion process. Where selective introduction of the 
dopant is wanted, this may be accomplished without application 
of masking and subsequent stripping, and without back-surface 
etching because of double-sided impurity penetrations. 

Thus, ion implantation can involve fewer handling or 
transferring operations than the diffusion process, and 
consequently can result in labor savings and increased 
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yields. However, ion . .^plantation requires an annealing 

step, which will be further discussed later on. It has 
( 8 ) 

been suggested to use ion implantation as an integral 
part of a total vacuum process sequence for fabricating 
solar cells after wafer or sheet generation. Such a sequence, 
although high in capital costs, could result in labor savings 
and high yields. 

The charge on the dopant ions allows for mass-spectro- 
scopic separation using magnetic fields, and for accurate 
measurement of the ion flux entering the deposition region, 
as lonq as nejtral and doubly charged particles are handled 
correctly. The ion beam currents can be readily measured 
by placing a Faraday cup in the beam's path, but this requires 
a preceding calibration to determine the fraction of uncharged 
and doubly charged ions. The mass analysis and ion current 
measurement features of the ion implantation process can pro- 
vide better control over the quantity and quality of the 
dopant than other processes, and can therefore be applied 

to obtain better process uniformity and repeatability. Dose 

(9) 

uniformities of + 5% (2o) are achievable 

Since ion implantation can be performed at or near 

room temperature, low energy implantations can result in 

o 

original dopant penetration of less than 100A. This is shallower 

than can be achieved in most high temperature source deposition 

steps in the diffusion process. 

Upon entering the substrate, the dominant interactions 

of the ion are wi th the electrons of the lattice, which 

slow the ion down through kinetic energy transfer. After 
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this initial slow-down to sufficiently low energies, i.e., 

2 

ion velocity less than Z^ /fi, collisions of the ion take 
place with the nuclei which completely stop the ion. In 
most cases, the stopped ion rests interstially in the 
crystal lattice. The largest impurity concentration is 
thus found at a penetration distance, "x ", from the 
surface. As a first approximation in the region where 
nuclear collisions dominate, the penetration depth is 
proportional to the square root of the ion beam energy. 
This penetration depth is described by: 


0.7 (Zj 273 + Z si 2/3 ) 


z z . 

I Si 


*1 + M si 


M, 


E i (X). III-2.5 


Ej is the energy of the ion beam in eV, Z and M refer to 
the atomic number and atomic weight, respectively, while 
the subscripts I and Si refer to the ion and to Si, re- 
spectively. The concentration varies from the penetration 
distance approximately according to a Gaussian distribution 
and the impurity distribution can be described by the 
empirical relationship. 


C (x) - C p exp (- (x-Xp) 2 /2 a R 2 ) . HI-2. 
C is the concentration at the penetration distance, and 
o- is called the standard deviation of the concentration 
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function, or the distance from x a t which the concentra- 
tion is equal to C p //e7 The peak concentration depends 
upon the ion beam current (i) and the implantation time (t) 
or 

10 16 i • t 

C D cm III-2.7 

' 4 °R 

The unit of i is nft, that of t is seconds, and o is 
is given in nm. 

The penetration distance can be calculated from 

electron and nuclear ionic collisions only if "channeling" 
does not occur. Channeling is the name given to the con- 
siderably enhanced penetration distance of ions which are 
aligned with low index crystallographic directions, and 
therefore travel parallel to and in between high atomic density 
crystal planes. Since ions travelling this path experience 
relatively fewer collisions with silicon atoms, they can travel 
further into the silicon. To avoid channclinn, the beam must 
be oriented at a slight angle (‘7°) from the orientation of 
the low index crystallograph ' c arcs. This increases the 
apparent distribution of atoms in the crystal plane 
normal to the ion beam's path, and thus increases the pro- 
bability of ion-nuclei collisions. 

The implantation process results in the displacement of 
the silicon atoms from their normal lattice sites by the ion 
collisions, thus creating "vacancies" and "interstitials". 

The implanted impurity atoms, predominately located at intersti- 
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tial sites, are not electrically active. Thus few impurity 
atoms which take up substitutional positions, tend to 
compensate the originally present impurity atoms of oppo- 
site dopant type, and to shift the Fermi level of the silicon 
towards the center of the energy gap. Annealing of the 
ion implanted wafers is required be . 1 ^ to reduce the mentioned 
crystal structure damage resulting from the implantation 
process, some of which is electrically active (recombination 
and trapping centers) , and to electrically "activate" the 
dopant impurity by moving its implanted atoms from inter- 
stitial to substitutional sites. This annealing is usually 
accomplished by a high temperature soak, called thermal 
annealing . 

Thermal annealing broadens the impurity profile, usually 
to a junction depth as great or greater than obtained by use 
of relatively low temperature, short tiir^ diffusions as they 
are normally used for solar cell production. Nevertheless, 
ion implantation followed by thermal annealing is capable of 
producing solar cells with efficiencies equivalent to those 
prepared using diffusion. A thermal annealing cycle of lh at 

450°C and 0.5h at 859°C has been repeatedly found to yield per- 
formance-wise competitive silicon solar cells. This, in part 

negates the potential advantage of being able to control the 
dopant profile at will by varying the implantation energy and 
dosage. Such "designed profiles" might lead to higher effic- 
iency solar cells than obtained so far. Electron and laser beam 
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annealing, as short transient annealing methods, have therefore 
been and are being investigated because of the potential 
profile maintenance as well as several other anticipated 
advantages. These advantages have, however, so far not 
been realized, and cells with efficiencies comparable to 
those obtained by the oven annealing process have so far not 
been reported. However, electron beam annealing, followed by 
a low temperature soak, at approximately 500°C, is now said 
to produce cells with efficiencies comparable to those using 
the "thermal annealing" process. The necessity of the low 
temperature soak seems to indicate that the electron beam 
annealing process does not reduce crystal damage or permit 
gettering as well as thermal annealing. Although the 

electron beam pulse anneal, followed by a low temperature soak, 
will be a more costly process than a higher temperature 
activation/annealing soak one, it appears attractive because 
of the greater freedom in selection and control of the impurity 
profile. If this pulse anncal/soak process, or some other, 
simpler process for activation/anneal i ng could be developed, 
so that ion implanted solar cells might attain higher 
efficiencies than cells prepared by diffusion processes, then 
ion implantation would become a most interesting process 
option, even at a possibly somewhat higher process cost than 
diffusion . 

With the attainability of a potential efficiency advantage 
of the ion implanted solar cells over the diffusion produced 
cells not demonstrated, the usefulness of ion implantation, 
as part of an LSA solar cell sequence will be determined by 
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the potentially achievable cost reduction. Currently, the 
high capital costs, the low reliability, and the low through- 
put rate of ion implantation machines, make junction formation 
with them too costly to be used for large scale solar cell 
production. Large cost reductions are, however, expected to 
be accomplished in the future (1986) by several approaches. 
Approaches to this end include the introduction of large 
throughput machines with high current, hot cathode ion beam 
sources incorporating an analyzer and more automated operation 

*7 \ 

through computer control' , and the development of ion im- 
planters with unanalyzed or roughly analyzed ion beams 
using hollow cathode sources. Some current and future 
applications of ion implantation are listed in Table I along 
with the conditions contingent to the two potential advantages 
of lower cost and higher efficiency. 

In consequence of this discussion, it has to be observed, 
particularly in reading the following sections, that ion 
implantation and activation/annealing are inseparably con- 
nected, and that, even though only one of the process steps 
may be mentioned, it cannot be usefully carried out without 
the other. But there are several ways of carrying out ion 
implantation as well as annealing, with various combinations 
of the two, so that one really deals with an ion implantation/ 
annealing submatrix of options. Consequently, both the pro- 
cess cost and the performance of the finished solar cells 
depend on the option selection within this submatrix. 
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EVALUATION OF ION IMPLANTATION 
FOR LSA PRODUCTION 


PRINCIPAL 


APPLICATION 

STATUS 

ALTERNATE PROCESSES 

pn JUNCTION FORMATION 

PROVEN: PERFORMANCE 

DIFFUSION 


EGUAL DIFF'D JCTN. 

CVD/ EPI 

BSF OR BACK HI/LO JCTN. 

CONCEPTUAL 

THICK FILM/ALLOYING . 
DIFFUSION; CVD/ EPI 

FSF OR FRONT Hl/LO JCTN . 

EFFECTIVENESS NOT 

DIFUSSION 


YET PROVEN 

CVD/EPI 

CONTACT METALLIZATION 

CONCEPTUAL 

THICK FILM 
ELECTROLESS PlATING 
VACUUM EVAPOR'N 
SPUTTERING 


ION IMPLANTATION FOR pn- JUNCTION FORMATION 

CONCEIVED 

ADVANTAGES CONDITIONS STATUS 


LOWER COST 


HIGHER CELL PER- 
FORMANCE THAN 
ACHIEVABLE BY 
ALTERNATE PRO- 
CESSES 


LIKELY ONLY IN SEQUENCE 
WITH OTHER VACUUM PRO- 
CESSES 

DEPENDS ON SUPERIOR 
IMPURITY PROFILE/ 

FEWER CRYSTAL DEFECTS 


TECHNOLOGY ADVANCEMENTS 
REQUIRED. 


STILL TO EE DEMON- 
STRATED 


Table III.l 
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4. Appraisal of Present-day Ion Implanters 

The application of ion implantation for pn or high/ 
low junction formation in process sequences for future 
large scale LSA manufacture depends on the fulfillment of 
either of two conditions: 1.) its costs are equal to 

or lower than those for pn junction formation using diffu- 
sion or high/low junction formation using alloying or 
diffusion, possibly in combination with each other or with 
other process steps; or 2.) the performance of the solar 
cells fabricated by use of ion implantation is adequately 
higher than that of cells prepared by other processes 
so as to justify a higher price. 

Ion implantation is currently used in semi-conductor 
industry production activities for implanting* in solid state 
devices, impurities of low dosage and relatively deep penetra- 
tion (high energy) . In order to gather information on the cur- 
rent state of production line ion implantation, we visited, 
among others, RCA-Somerville, where a Varian-Extrion 200-1000 

ion implanter is used for integrated circuit manufacture, as 

( 12 ) . 

well as for solar cell fabrication in pilot operations. 
Implantations are routinely performed at beam currents rang- 
ing from 0.1 pA to 1.5 mA, at voltages up to 100 kV, alternating- 

+ + 4 - 

ly with P , B and As ions, in a 24 hour-a-day schedule of 
5 to 7 days-a-week. 

The Varian-Extrion 200-1000 ion implanter is available with 
a semi-automatic cassette wafer feeding mechanism that allows 
continuous processing, ii. creasing its output rate to 300, 7.62-cm 
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round wafers per hour. In order to achieve this output rate, 

the ion implanter also has to be modified to operate in a high 

currer.L (4 mA), low voltage (<25 keV) mode. These options are 

included in a Varian -Extrion 200-1000 iinplanter in operation at 
( 2 ) . . 

Spire. Additional options provide an off-axis beam tilt 

to minimize channeling. To achieve dose uniformity and avoid 
shadowing from a tilted beam on a texture-etched surface, 
the wafer is rotated about its axis at 1 rev/sec. 

The cost of such a machine is approximately $315,000 and 
it requires one full-time operator. To achieve acceptable 
machine operation, the RCA personnel have found it necessary 
to have a skilled technician stationed within the immediate vi- 
cinity of their ion implanter at all times, and to make adjust- 
ments in the machine operating parameters quasi-continuously . 

They believe that computer controlled functions, similar to 
those proposed by Spire ^ in their ion implanter design, 
could considerably reduce the need for continuous skilled 
attendance. They mentioned, however, that designing adequate 
computer controls might be difficult since, so far, adequate 
sensing of the status of all parts of the machine and of the 
parameters affecting its operation does not exist. Thus, 
correctly operating the ion implanter is still more of an art 
than a science and requires the adjustment of many functionally 
interrelated controls. Similar statements were variously 
heard, summarized ny Varian-Extrion personnel, in the remark that 
successful machine operation depends very heavily on the operator, 
and that wide variations are experienced among the various users. 
RCA personnel has found that leakage from the high voltage 
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machine elements, in part due to condensed source material, 
tends to interfere with the sensitive dose rate measurements and 
the machine control. Other problems resulted from persistent 
leakage of cooling fluid which could be reduced oy the use of 
freon in lieu of the more common deionized water, albeit at 
significantly higher costs for the make-up fluid. 

One of the major problems mentioned at RCA and elsewhere, 
is the deposition on many parts of the machine of atoms of the 
implanted species as well as of material sputtered off the 
various parts of the source. Arsenic is especially troublesome 
in this respect because of its relatively low vapor pressure 
compared to other implanted species. This deposition occasion- 
ally results in electrical malfunctions, such as shorting of 
insulators and arcing, which occasionally has led to power 
supply or logic board damage. The machine, therefore, requires 
frequent thorough cleaning of the affected regions. Phosphorus 
also condenses on the machine's interior, and we have heard 
of short phosphorous fires upon opening the machine. 

Much of the unscheduled maintenance is performed under 
service contracts. RCA personnel mentioned that such a 
service contract with Varian-Extrion has an annual cost of $13,000. 
This contract provides the so-far extensive on-location servicing 
by Varian-Extrion personnel and replacement of failed parts, 
frequently circuit boards. RCA personnel estimates that about 
two-thirds of this money covers time and expenses of the service 
personnel, and the remainder replacement parts. RCA has re- 
cently introduced regular scheduled maintenance of their Varian- 
Extrion 200-1000 "high current" implanter for which 4 hours per 
week are allocated. During these maintenance periods, the 
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machine interior is cleaned, filaments, if needed, are replaced, 
vacuum pump oils are changed, the machine inspected, and poten- 
tially unreliable parts identified and replaced. Since this 
institution of preventive maintenance, the previously frequent 
machine breakdowns have decreased to a tolerable level. At RCA, 
the experienced filament lifetime, as plotted on Figure III. 2, 
is in the 60 to 120 h range for an average ion beam current of 

around 0.75 mA, although much implantating is done with a 1 mA 
( 12 ) 

beam current. 

Because of the relatively frequent machine breakdowns 
of ion implanters, RCA's personnel have found it necessary 
to keep an extensive spare parts inventory, so that bad or 
suspicious parts can be replaced with minimal machine down time, 
in order to maintain production schedules and to reduce the 
impact on operating costs which are heavily influenced by the 
high cost of the equipment. 

An ion implanter has been in operation at 
Western Electric since 1974 in a production line, high through- 
put mode. This implanter, called the Pk- 30, has an output 
rate of 450, 7.62-cm diameter wafeis/h at a dosage of 
1x10^ ions/cm 2 . The machine can accomodate either a hot 
cathode. Freeman- type source, or a cold cathode (Penning) 
source. It operates in a low voltage (30keV) mode. In the 

case of the cold cathode source, a phosphorus current of 

(14) 

5 mA is obtained, with a source lifetime of 40 h. The wafers 

(7.62-cm diameter) are placed on a disk, 30 at a time. The 
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FIGURE in. 


disk is mounted horizontally in the ion implanter. After pump 
down to 10 5 torr, which takes approximately 3 minutes, the 
disk is rotated at 900 rpm while the undersi.de of the wafers is 
exposed to the fixed ion beam. The total time of each run is 
approximately 4 min. 

The PR- 30 is physically, a relatively small machine. The 

implantation unit, without controls, occupies a floor area of 

1.8m x 2.1m. Two standard instrument racks house the control 

units. The PR-30 Is used only in Western Electric factories, 

and it is not sold on the open market. We have been given an estimated 

price for this machine, if it could be marketed, of less than 
(15) 

$300,000. 

A high current (10 mA) and low voltage (10-50 keV) ion 
implanter, designated NV-10, is currently being readied by 
Nova Associates for introduction into the marketplace. The 
machine uses a Freeman, hot cathode source, with an expected 
lifetime of 16 h at 10 mA. The machine costs approximately 
$410,000. Its output rate for a 2 x 10^ icns/cm^ dosage 

of 270 wafers per hour of 7.62-cm diameter is limited by the 
wafer feed mechanism. If a Tester feed mechanism could be in- 
stalled, the output could be increased to 3-4 times the 
present one, to take better advantage of the machine's high beam 
current. The wafers are mounted, 18 at a time, in a disk that 
is rotated in a near vertical plane during implantation. The 
stationary beam is approximately 1cm x 2-3cm. As with the 
Western Electric implanter, wafer rotation eliminates the need 
for magnetic or electrical beam scanning. 


1 11-24 


2 

An add-on process price of $38. 96/m for implanting 

I *» * 

phosphorus with a 1-2 x 10 ions/cm' dose was calculated lot 

the modified Varian-Extrion 200-1000 machine. This price 

includes the cost of the silicon sheet lost-in-process. The 

sheet price used applies to silicon wafers which have been 

texture etched on one side. The slicing cost was taken from our 
(17) 

previous study of current production slicing costs (HAMCO 

ID data) . 

The add-on process price for ion- implantation using the modified 
Varion-Extrion 200-1000 is low compared to other prices calculat- 
ed for currently used ion implanters. For instance, the calculated 
add-on process price for the Varian-Extrion 200-20 A machine 
is $303. 42/m 2 This high price is due to the machines 's 
low throughput rate as it was designed for high voltage, 

low current (under 0.2 mA) operation. Its hourly output rate 
therefore is only 10, 12-cm diameter cells. 

It should be noted that the given add-on price calculation 
for the modified Varian-Extrion 200-1000 implanter is based on 
experimental, not production line data. Therefore, this value 
does not reflect the breakdown or maintenance problems ex- 
perienced by ion implanters in production operations. However, 
reliable detail data are not yet available for the cost 
components of regular production ion implantation, since 
this process was only rather recently introduced as a production 
process. Still, if such data would be available, they would 
not represent the ultimately achievable costs, after machine 
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and process maturity have been attained. While efforts are 
in progress to adapt the ion implanters better to production 
line operation by increasing their throughput rate , mechanizing 
their operation and improving their reliability, it will 
be some time before the process will be a mature production • 
operation with similar costs experienced by the various users. 
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5, Technology Development for Future Ion Implantation Machines 

The realization of the 1986 cost projections for ion im- 
plantation is contingent on several improvements in the tech- 
nology of ion implantation machines. For one, the ion beam 
current has to be increased significantly to achieve econo- 
mically acceptable throughput rates. Also needed to be in- 
creased is the lifetime of tne source, in terms of mAh's, 
to avoid excessive costs from changing and rebuilding the 
sources, as well as macnine downtime. To reduce skilled labor 
requirements, the implanter's controls should be as automatic 
as possible. In addition, continuous or semicontinuous wafer 
feed, along with appropriate vacuum pumping mechanisms have 
to be employed. Also, care has to be taken in the mass analyses 
and the control of large current/small voltage ion-beams needed 
for solar cell fabrication, because space charge effects make 
those operations difficult. In some LSA process sequences, 
ribbon material is planned to be the substrate. Since rotation 
of elongated rectangular workpieces about their axis is impracti- 
cal, other procedures to achieve uniform deposition have to 
be utilized in the future implanter, e.y., magnetic or mechani- 
cal beam scanning. 

As mentioned previously, at present, PN junction 

formation using open tube diffusion is a small cost contributor 

( 18 ) 

to the solar cell module cost, constituting approximately 1%. 

A replacement process for diffusions in future LSA process se- 
quences would require lower costs, or yield hiyner performing 
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cells, or offer a simplified fabrication sequence. Implanta- 
tion costs are expected to be lowered dramatically by increasing 

the ion implanter's throughput rate from about 2 m 2 /h to 
2 

nearly 200 m /h. To accomplish this, the total ion beam 
current, flux rate of ions impinging on the silicon, is ex- 
pected to be increased from 4 mA to 100 mA. If multiple sources 
are used, then the ion beam current per source needs to be 
increased by a factor of 4 to 5. Increasing the beam current 
will, in general, increase the implanter’s output rate in 
the same ratio. But, as shown on Figure III. 3, the increase in 
the machine's cost per unit beam current decreases with beam 
current. In Figure III. 3, the experienced machine cost per unit 
beam current is plotted as a function of the beam current to- 
gether with an extrapolation to tne future. The first four 
open circles reflect the costs of ion implantation machines 
that are in operation and the solid circles reflect projected 

data from the listed organizations. 

In addition to larger ion sources, future implanters 

would have to be more reliable than current ones. The hi b 
capital cost of ion implanters necessitates their utiliza- 
tion rate to be as high as possible. Proposed future machines 
(Lockheed, RCA, Spire) have been projected to h*ve utiliza- 
tion rates between 85-95?) as opposed to today'. . For 

Motorola's unanalyzed ion beam implanter, the uptime fraction 
is not as significant because of" its relatively low cost . 

The Motorola machine is expected to cost $85,000 as opposed 
to at least $500,000 for any of the other three proposed 
machines which employ analyzing maqnets. One reliability 
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improvement is expected from increasing the source lifetime 
in terms of mAh's with beam current. However, although the 
source life expectancy de-creases as beam current output in- 
creases, as shown in Figure III. 2, the product of source current 
and lifetime increases with increasing output. Therefore 

more silicon can be processed between filament changes. In 

, ( 7 ) 

one proposal , multiple and spare sources are employed 
so that they could be replaced while the machine is operating. 

As listed on Table III-2, the source lifetime (mAh) is expected 
to increase in the future by a factor of ten. 

Another projected improvement is the reduced dependence 
of the ion implanter's performance on operator skill. At 
present a skilled operator needs to monitor the operating 
ion implanter continuously to achieve optimum output rates. 

These skilled labor requirements are expected to be decreased, 
in future implanters, by simplifying the machine's operation, ^ 
by larger batch loads, or by using microprocessors. 

It is thus hoped that future implanters could be operated with 
unskilled labor, with skilled labor called upon only occasionally 
for mechanical and electrical servicing. 

Since annealing is an integral part of the implantation 

process, studies are being conducted in the -LSA pro- 

(19) 

gram on an optimum process. Processes studied include 

thermal, electron pulse, and laser annealing with only thermal 
annealing yielding solar cells of comparable efficiency 
to those prod- ;d from diffusion. Thermal annealing costs, 
as mentioned previously are significant compared to those 
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TECHNOLOGY PROBLEMS TO BE RESOLVED FOR 
SUCCESSFUL LOW COST ION-IMPLANTATION 


5-FOLD INCREASE OF BEAM CURRENT PER SOURCE (A mA -» 20 mA) 

10-FOLD INCREASE OF SOURCE FILAMENT LIFE (25 mAh -» 280 mAh) 
REDUCED RELIANCE ON OPERATOR INFLUENCE FOR EFFICIENT MACHINE 
PERFORMANCE 

REDUCED FREON LOSSES FROM COOLING SYSTEM (HIGH VOLTAGE.) 

EASIER CLEANING OF SPURIOUS MATERIAL DEPOSITED IN SYSTEM. 

(DEPOSITION PROBABLY NOT AVOIDABLE.) 

UNIFORM DEPOSITION W/0 WORKPIECE ROTATION 
REDUCED CAPITAL COST (CURRENT SINGLE SOURCE. 2 mA MACHINES 
CAPABLE OF 200 KAFERS/H COST ~$0.5 Mill.,) 

IMPROVED ANNEALING METHODS (PULSE ANNEALING?) 

REDUCED ENERGY CONSUMPTION, 

ORIGINAL PAGE IS 

OF POOR QUALITY 

ALTERNATE TECHNOLOGI ES: 

COLD CATHODE/ SOLID SOURCES (simpler system, but: faster source 

EROSION? MORE SPURIOUS DEPOSITION 
IN SYSTEM THAN FROM GASEOUS SOURCES?) 

OMISSION OF ANALYZING MAGNET (capital cost and energy savings. 

but: impurity profile acceptable? 

SPURIOUS IMPURITIES CONTROLLABLE?) 


Table III-2 
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for future ion implantation processes. Also thermal annealing 
decreases the potential efficiency of ion implanted cells. 

The shallow implanted PH junction depth that can be obtained 
from implanting with low energy ions has the potential 
of yielding better performing cells than those from the dif- 
fusion process, because of greater UV-response. However, 
the thermal anneal cycle broadens the shallow implanted 
junctions depth, making it comparable to that obtained 
using gaseous diffusion. 

For an effective use of the ion implantation process, 
an extended, automated, vacuum, production sequence has been 
proposed by Spire. For this sequence to be practical., the 
annealing process has to be performed in a short time interval. 
Since the conveyor belt, in the Spire sequence, moves at a rate 
of 30 cm/sec, a thermal annealing cycle of only 5 minutes 
would require an effective furnace length of 90 m. Electron 
or laser beam annealing would be compatible with a rapid 
production line, since they can be performed in fractions 
of seconds . ^ 19 ' However, solar cells annealed with either of 
these two techniques show a decreased performance. A summary 
of some other technical problems that need to be solved for 
the successful implementation of ion implantation for future 
solar cell manufacturing processes is listed in Table III. 2 
These problems include uniform deposition of ribbon-shaped 
wafers, more effective coolant usage and convenient removal 
of deposited source material. 

The importance of beam current size to implantation output 
is shown by the expression for the unit area ion implantation time 
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[1.602 x 10 (a*sec/ionJ • ion flux (ions/cm ) - 

tp ■ k — sec/cm . 

Ion beam (amps) 

ill-2. 8 

The proportionality constant, k, is < 1 and depends on the 
degree of overscan and the beam utilization. Therefore, 
as a first approximation, the throughput rate of an ion 
implanter is proportional to its beam current. Because the 
implantation process is capital intensive, lowering the machine 
cost per unit beam current will lower the implantation cost 
in about the same ratio. As can be observed in Figure III. 3 
the machine costs norma li zee. to their beam current are ex- 
pected to decrease approximately proportionally with increased 
beam current. For future ion implanters, a large capital cost 
decrease per unit of output is anticipated by increasing the 
beam current without proportional increases in machine costs. 

There are several approaches for increasing the ion 
beam current. One approach, proposed by Spire, is to increase 
the size and number of the hot cathode sources to 20 mA and 
10, respectively. The source current lifetime is increased by 

changing from a Freeman to a Chavet type filament. Higher 
currents are tolerable in the latter source, because the 
Chavet filament is looped and therefore is not as heavily 

degraded by the back ion bombardment. Although source 
lifetime does decrease with increasing currents, as shown 
in Figure III. 2, this decrease is less than the increase in 
current. In another approach, a hollow cathode source, similar 
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to that used in ion beam thrusters, is proposed. v ' 
source is expected to yield a current of 100 mA, but be- 
cause of the non- col lima ted, large crossection nature of the 

beam it cannot be mass-analyzed. In an ion implanter proposed 

by RCA, two 10 mA ion beams are used simultaneously. One 

is used to implant the front of the wafer with phosphorus 
15 -2 

at a 1 x 10 cm dosage while the other implants boron with 

a dose of 5 x 10"^ cm ^ The Lockheed proposal has ore 

10 mA beam that can process about one 7.62-cm diameter wafer/sec- 
(5) 

ond. The wafers are loaded and unloaded to and from 

4 side chambers which surround the central implant chamber. 

In the proposed Spire machine, 7 of the 10 sources are 
operated simultaneously with six running at a current of 
16 mA, and the seventh at 4 mA. The ion beam from each source 

passes through a collimator with a slit geometry of 2 x 75 mm 

to provide mass analysis. The larger six sources are broken 
into two sets with an analyzing magnet for each set. Three 
ion beams strike the moving silicon wafers at +15° to the normal 
and three at -15°. The wafers are transported on 20x20 cm carriers 
on a belt moving at a rate o r 30 cm/sec. The seventh and 
smaller ion beam is used for a final dose control. The three 
remaining sources are used as spares. As plotted in Figure m. 2 , 
it is expected, by Spire, that the average source lifetime 
can be increased to 24h, or approximately 400 mAh. This 
would mean, on average, a source replacement every 4 h with 
each replacement requiring 10-15 min. labor. A "dead" source 
is expected to be ready for replacement within 24 h. The 
implantation energy is designed to be 10 keV, dose uniformity 
to be + 10%, and analysis ;o + 0.5 AMU. In order not to 
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enhance the space charge effect of the large beams, electric 
fields after the extraction gap are avoided in the Spire machine 
design. The scanning deflector, shown in Fig. III.l, is 
operated magnetically. 

The narrow width of the 16 mA ion beams makes them analy- 
able since the radius of curvature of the ion beams caused 
by the magnetic field can be made larger than the beam's 
width. The radius of curvature is given by: 

r = ( m / e) x v / B III-2.9 

where (m/e) is the ion's mass to charge ratio, v is the ion 
velocity, and B is the magnetic field strength (in gauss) . 

If a linear magnetic field is assumed, then the deflection 
angle is sin ^ |Te/m) x (B/v) "£J, where i is the length of the 
magnet. Therefore, the angle of deflection depends on the (e/n) 

ratio. A slir. in the ion path placed preceding the beam selects 
the desired ion, as seen in the top drawing of Figure III.l. 

A large, transient temperature increase ('~'800°C) can 
cause considerable stress in the silicon wafer, and make sub- 
strate movement and handling difficult . The energy flux den- 

15 - 2 

sity J, of a 10 keV ion beam at a density of 1 x 10 ions/cm , 

is 1.602 j/cm 2 . With the proposed output of the Spire implanter 
of 180 m 2 /h, the implantation time is 0.002 sec/cm 2 . The 
temperature rise of implanted wafers is given by, 

AT = J/C p p(2D t p ) i/2 °C, III-2.10 

where C p is the specific heat of silicon (0.71 j/g°C)^ 21 ^ at RT, 
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and p it. the silicon density (2.34 g/cm 3 ) . D is known as the 

heat diffusivity which is equal to k/C p, where k is the thermal 

conductivity of silicon (1.47 j/sec*cm °C) . The expression 

/2Dt is the thermal diffusion length and cannot exceed 
P 

the wafer's thickness. Equation III-2 . 10 is valid when the im- 
planted junction depth is small compared to /2Dt~ or the 
thickness of wafer. This condition is satisfied for NP 
solar cell NP junction formation. The junction depth is normally 
approximately 0.2 urn and the implantation time is sufficiently 
long to make the diffusion length several hundred microns. 

For implanting a 200 pm thick wafer, with the porposed Spire 
machine, the temperature increase over the environment is ex- 
pected to be 48°C. 

In Motorola's proposal, shown in the bottom drawing of 
Figure III.l, a large ion current beam (100 mA) is obtained from 
a hollow-cathode source derived from ion thruster technology. 

Ion thrusters, using ionization of mercury, have very large 
beam currents (several amps) , and lifetimes of thousands of 
hours. It is thought not to be difficult to modify the 
thruster to ionize phosphorus or other suitable dopants. 

However, the ion thruster beam can not be mass-analyzed 
because of its circular cross-section and large diameter. 

The dispersion caused by a magnetic field would be less than 
the beam's diameter. In addition, the energy spread of ions 
emitted from a ion thruster type source hinders good magnetic 
separation, since the curvature radius of ions under the influ- 
ence of a magnetic field is directly proportional to its 
velocity. The effect on solar cell efficiency of implanting 
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with an unanalyzed or a "roughly" analyzed beam is not yet 

(23) 

known and investigations have just been initiated. The 

proposed Motorola ion implanter is fairly simple in design; 

the wafers are transported (past the ion beam) by a belt 

through differentially pumped vacuum chambers. Dose uniformity 

might be a problem, because of the Gaussian distribution of 

the beam's intensity and an individual wafer might be exposed 

to only a selected portion of the ion beam. It takes less 

than 0.75 sec. to implant a 12-cm diameter wafer with a 
15 -2 

2 x 10 cm dosage of phosphorus with a 100 mA beam. The 
low capital cost of this implanter, makes the Motorola pro- 
posed ion implanted process the lowest cost one studied in 
this report. 

In the RCA and Lockheed proposed machines, hot cathode 

(3) 

ion sources are employed. In the RCA-proposed machine, both 

the PN and PP+ junctions are formed simultaneously by using 

two separate 10 mA beams. One beam is used for phosphorus and the 

other for boron. This machine can process approximately 
2 

100 cm /sec, and allowing time for beam scanning and beam 
loss at edges, the machine's throughput is 2000, 7.62-cm 
diameter wafers per hour. The wafers are transferred auto- 
matically from 500 wafer cartridges to 50 wafer cassettes from 
which they are then removed to a holder for implantation. 

The high capital cost of the RCA implanter relative to its 
output, makes the RCA process the most expensive of the future 
implantation process projections. 

The Lockheed proposed machine uses a 10 mA beam, and can 
implant 3000 wafers/h (7.62-cm diameter) . The wafers, 
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which are batch-loaded, are held in 1200 ring-shaped trays 
or carousels (50 wafers/tray) that are stacked and distributed 
among 4 cylindrical vacuum chambers adjacent to the implanta- 
tion chamber. During the implantation process, the trays are 
transferred to the central chamber where they are rotated 
such that each wafer is scanned on its underside by the ion 
beam. This is repeated 4 times for each tray to assure dose 
uniformity. The ion beam is kept constant at 7° to the normal 
while the wafers are rotated. This eliminates the need for 
electrical or magnetic beam scanning. After all the wafers 
in the machine have been scanned, vacuum in the implantation 
system is broken and the wafer loading cylindrical chambers 
are replaced. It takes approximately 20 hours for the com- 
pletion of one run: 2 hours for loading, 16 hours for pro- 
cessing, and 2 hours for unloading. The Lockheed process 
employs phosphorous pentaf louride (PF^) as the source gas, 
instead of PH 3 or P. Phosphorous pentaf louride is very 
expensive and is a large cost contributor (about 16%) to 
the add-on process price. 
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6 . 


Junction Formation Material, Labor, and Capital Require - 


ments Cost Structures 

The costs of present and future junction formation pro- 
cesses, broken up into their material, labor, capital, over- 
head, and return-on-equity components, are summarized in 
Table II. 3. Also listed in Table III. 3 is the throughput rate, 
in terms of number of wafers processed per hour and their dia- 
meter. The cost calculations are based on ■ ® SAMICS method- 

( 24 ) 

ology . The detail process parameters, including direct 
and indirect material, labor, and equipment and facility re- 
quirements are presented in Tables III-4A to III-4C for the 
diffusion processes, and in Tables III-5A to III-5C for the 
ion implantation processes. These data are given both in 
terms of unit area of cells processed and of net plant 
operating time. The hourly consumption rates were converted 
to unit area data by use of the "effective output rates" 
shown in line 1 of Tables III-4A and III-5A. The effective 
output rate is the produce* of the operating machine's through- 
put rate and the usage fraction, (process "up-time" divided by 
plant operating time) . The plant operating hours were taken 
to be 8280 h per year. For calculation of the material costs, 
the unit prices shown in Table III. 6 were used. The labor 
costs were obtained by employing a labor rate of $3.894/h for 
the semiconductor assembler (SAMICS B3096D and of $5.29 
for the electronics technician (SAMICS B3704D). The total 
labor costs include an indirect labor charge of 25%, and a 
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* Refers to two ■achlnes operating to 


5.26% addition for overhead expenses, both according to 

SAMICS. The equipment and facility costs were similarly 

(24) 

calculated consistent with the SAMICS methodology. 

The cost of the wafers, which are reflected in the lost- 

in-process cost (Table III. 3, line 18) are taken from our previous 

(17) 

studies of slicing processes and the 1986 silicon and sheet value 

(25) 

goals listed in JPL-LSA's price allocation guidelines. In 
addition to slicing, the cost of one-sided texture etching 
is included in the current and future wafer prices. The 
etching is performed by applying etch stop in the form of 
wax on one surface, texture etching with 30% NaOH at 90°C, 
and removing the wax with plasma etching. The etching step 
costs have been derived from information published by 
Motorola, and add up to approximately $3.09/m^. The cal- 
culated prepared wafer prices are $350.98/m2 and $41.21/m2 
for 1978 and 1986, respectively. The specific process for 
the current wafer price is slicing 10.16-cm diameter wafers 
with a HAMCO ID saw. 

The first two columns of Table III. 3 refer to current 
implantation and diffusion techniques, while the other columns 

detail the costs of proposed processes. Two multi-step se- 
quences for producing front and BSF cells are also shown on 
Table III. 3. The 5-step Motorola diffusion process, which 
is detailed in Table III. 7, consists of protecting the front 
surface by spinning-on-silica, diffusion of the BSF using 
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Table 1 1 1 - 6 


Material Prices Used for 
Calculating Costs 


Item 

Unit Price ($) 

Source 

Direct needs 
500 ppm of PH^ in N 2 

r 

42.86/m 3 

Spectrolab 

PH 3 gas 

0.991/* 

Motorola 

POCI 3 

0.0204/g 

SAMICS E1504D 

5% PH 3 in H 2 

0 . 02 9/H 

SAMICS E1472D 


0.0051/ml 

Matheson 

SeG-phosphorus 

2.76/g 

Alfa Products 

Indirect Mater 

Lai 

(Ventron) 

Liquid nitrogen 

0 . 202/1 

SAMICS C1080D 

Nitrogen gas 

0 . 0004/1 

Obtained from LN 2 

Argon 

0 . 005/1 

SAMICS E1112D 

Oxygen gas 

0. 000184/j, 

SAMICS E1448D 

Compressed air 

— 

SAMICS C2032D 

Cooling water 
Electricity 

jO. 566/kWh 
U. 996/m 2 
0. 0319/kWh 

SAMICS C112BD 
SAMICS C1016B 
SAMICS C1032B 
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BCljf a spin-on silica protection of the back surface, 
phosphine diffusion, and stripping of silica from both surfaces 
with a 4:1 NH^OHjHF solution. The result is an N+PP+ wafer 
with no silica coating, ready for metallization or AR-coating. 

The other multi-step process consists of RCA's double-sided 
ion implantation followed by thermal annealing. The RCA 
2-step process yields wafers equivalent to Motorola's 5-step 
wet chemical sequence. 

The cost components for activation annealing are shown 
in Table III. 3 because it is presently a necessary step after 
ion implantation to achieve state-of-the-art performing cells. 
Annealing costs are significant compared to those derived using 

the high throughput implanters proposed by Motorola and Spire. 

The major cost components from Table III. 3 are graphically 
represented on Figure III. 4. In addition. Figure III. 4 includes the 

cost of RCA's proposed gaseous diffusion using POCl^. This 
diffusion process takes approximately one hour and has an 
output rate of 2,000 7.62-cm diameter wafers per hour. 


The prices for the proposed PH 3 and POCI3 diffusion processes 
are $3.01/m^ and $3.86/m^, respectively. These two processes 
should be available for near term production sequences; no major 
technical problems need to be solved for their applicability. 

Ihe cost decreases for the diffusion processes are about a 
factor of four lower from current ones, and for the most part, 
depend upon throughput increases. The higher output rate 
for Motorola's diffusion process, as compared to Spectrolab * s , 
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Yielded Add-on Cost Components for Motorola's 5-Step 
Wet Chemical Front Junction and BSF Sequence '($/«*) 
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is based on processing larger wafers (12-cm) at the same 
rate as smaller ones (7.62-cm). The larger area (2.48x) 
of the 12-cm diameter wafers accounts for the higher through- 
put rates of Motorola's process. RCA proposes to automatically 
transfer wafers from cassettes to furnace boats, and load/ 
unload the boats from the furnace, to increase output rates. 

The loading and transferring machines add to the capital cost 
of the RCA process, but increase output sufficiently to 
lower unit area costs. 

The RCA 2-sided implantation process, which is included 
on Figure III. 4 with annealing, and the Lockheed implantation 
proposal, could be ready for near-term production ( 1982— 

1984). Both these machines have 10 mA ion beams. The RCA 
implanter actually has two 10 mA beams but only one is used for 
the front junction formation. This beam size is only twice 
as large as some machines in operation and a 10 mA machine, 
the NV-10 , by Nova Associates, should be introduced into the 
marketplace shortly. 

The processing costs from employing the high current (100 mA) 
machines by Motorola and Spire are the lowest ones listed 
on Table III. 3 for junction formation. However, a longer time 
than for the other options discussed above will be needed 
before these machines are suitable for production use, because 
of larger extrapolations of ion currents and throughput 
rates. In addition, the Morotola 100 mA proposal has redactions 
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in labor, and capital costs because of its greatly simpli- 
fied operation, it employs an unanalyzed beam from a hollow 
cathode source, thus eliminating the need for any acceleration, 
magnetic, and scanning capabilities. The hollow cathode 
source originally designed for space propulsion use in ion 
thrusters, should give the high currents and lifetimes necessary 
for a low-cost, high throughput operation, needed in solar cell 
manufacturing. But, the effects on cell performance by 

implanting with an unanalyzed beam are unknown, although 

(23) 

investigations have recently been initiated. Spire expects 
their implanter to have a 100-fold increase in output rate 
over current machines. This is to accomplished, for the 
most part, by increasing the beam current to 100 mA, by 
having a continuously pumped, belt system, feed mechanism, and 
by increases in machine reliability by extensive use 
of microprocessors and redundant beam sources. 
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7 . Conclusions 


In order for the front junction formation processes 
involving diffusion and ion implantation to fit into 
future (1986) low-cost solar cell fabrication sequences 
their costs will have to decrease by factors of approximately 
four and ten, respectively . At present, the phosphorus diffu- 
sion process cost is $12.74/m2 while the ion implantation 
of phosphorus costs $37. 86/m . It is anticipated that the 
future cost contribution for front junction formation would 
be less than $3. 20/m 2 . 

The costs in the long term ion implantation projections 
by themselves seem significantly lower than those of the diffu- 
sion processes, but adding the cost of the necessary activa- 
tion annealing makes the costs comparable. For combined front 
and BSF sequences, the cost difference between a wet- 
chemical process (the 5-step Motorola sequence) and an equi- 
valent multi-step process employing ion implantation, is about 
$l/m 2 . The closeness of these two projections makes it diffi- 
cult to judge which would be economically advantageous in 
1986. From our calculations, it would appear that ion im- 
plantation and diffusion could be competitive. 

Future junction formation processes will have to fit 
well into high volume process sequences. Even though 
currently, gaseous diffusion is an inexpensive step in 
manufacturing solar cells, its costs have to be reduced 
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even more to fit into the future LSA framework. Cost 
reductions depend upon larger throughput rates to be achieved 
by processing larger diameter wafers, and by automatic wafer 
transfer. Wafer transferring could be accomplish* using 
specifically designed machines, which would increase the 
capital cost of gaseous diffusion. Another cost reduction for 
diffusion is related to quartzware (boats and tube liners) 
cleaning using mild chemical etching. The cleaning is necessary 
to minimize wafer contamination and is currently a significant 
cost contributor to the diffusion process. The required cleaning 
frequency and alternative cleaning procedures should be investigated. 

Ion implantation has recently been introduced into 
production activities, and its state-of-the-art performance 
is rapidly changing. During the last decade, ion beam current 
(and consequently the throughput rate) has increased by a 
factor of 1,000 - from a few microamps to a soon to be intro- 
duced 10 mA. For low-cost solar cell junction formation, 
the implanter's beam current would have to increase by an 
additional order of magnitude and its cost reduced by 
approximately a factor of 20. The feasibility of achieving 
these goals cannot, at this time, be assured. But certainly, 
activity in this area should be continued. 

If ion implantation's cost reductions could be accomplished 
through larger throughput rates and greater reliability, and 
if a compatible annealing process could be perfected, then ion 
implantation would be a strong candidate for junction formation 
in future LSA process sequences. 

The cost reductions required for gaseous diffusion to meet 
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the LSA future price goals are not as dramatic as those 
needed for ion implantation, and are not as dependent on 
technical development. However, studies should be continued 
in automatic wafer handling and in quartzware cleaning methods, 
because these are potential add-on price reduction areas for 
diffusion. 
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IV. Methodology for Energy-Cost Effectiveness Evaluation of Sybsystem 
Design and Manufacturing Process Options 

1,.-* Introduction 

One of the important attributes of a photovoltaic solar energy conversion 

system is its economic viability. The evaluation of this attribute is regularly 

performed in decision making about the use of such a system in a particular 

application) as well as in comparing the merits of one particular system design 

or solar cell production process against another. The key aspect in such an 

evaluation is the comparison of the cost of electrical energy produced by the 

photovoltaic system with the cost of competitively available electrical energy. 

The unit cost c_ of the electrical energy delivered from the photovoltaic 
Etl 

system can be expressed, following ref. (1), as: 

c = _Vp C £ a P . . ($ kWh” 1 ] (IV. 1) 

En ^d 

where C are the annual operating costs [ $ y 1 ] , C is the capital spent 
op cap 

in acquiring the system [$] and is the equivalent annual cost of capital 

[y 1 ], This equivalent annual cost of capital may, outside of the usual com- 
ponents of interest, taxes, depreciation, etc., include such considerations as 
desired profit, present or discounted value of life cycle costs, inf lation .ad- 
justments, etc. As the "fuel" in a solar energy utilization system is "free", 
the operating costs are essentially reduced to the maintenance costs, at least 
for the smaller distributed systems. And since it is generally assumed, in the 
absence of information to the contrary, that the systems will be designed and 
built for high reliability and thus require little maintenance, the maintenance 
costs are usually neglected in comparison to the costs of the capital, is 

the electrical energy usefully delivered during a year from the photovoltaic 
system to the load. 
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c = y ~ Cm 

En 'cap 


Y • r ; |$ kWh' 1 ] 
cap 1 


(IV. 2) 
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En E Ld 


As Y is a constant for a particular company at a given time, but will 
cap 

differ from company to company, the system dependent energy cost determinator 
is really the quantity: 


c i 

r = “ E ;[$kWh"y] (IV. 3) 

Ld 

which is the ratio of the required investment to the energy delivered per year. 

The evaluation and optimization of this quantity is therefore of primary interest. 


2. The Energy Delivered to the Load 

The energy E^ delivered during the year to the load is clearly related 
to, although different from, the energy E q delivered by the photovoltaic array 

itself to the remainder of the system. For a photovoltaic array of total 
2 

exposed area A^m ], E q is given by: 

8760h 

E o = A Ar/ H(t > n Ar dt; [kWh*y -1 ] (IV. 4) 

0 

where n (H(t), T(t)) is the effective array efficiency in the time interval 
Ar 

dt around time t, with n. being dependent on the temperature T(t) of the array 

Ar 

-2 

and on the irradiance H(t) [kW m ) during that time interval, as well as on 
the varying spectral distribution and the angle of incidence of the light. <J>(t) 
is a factor of magnitude between zero and one, which describes whether, or 
how much, energy can be delivered by an array for transfer to the load or to 
storage, depending on the existence of load and on the stains of the storage 
system during the respective time interval. Eq. (IV. 4), being a definite integral, 
can be expressed as: 
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4 *rV V. .td f U • ^OilkWh-y' 1 ] 


(IV. 5) 


following the custom of referring the output to "nameplate rating", or peak 

power output capability, which is, for the solar array, expressed by the 

product of the expected peak irradiance H pk and the array efficiency n Ar std 

measured under standardized conditions (including the peak irradiance H^)* 

The connection to eq. (IV. 4) is made via the "load factor" which is the ratio of the 

output actually delivered during the year to the "nameplate rating." f ^ is 

usually determined from the results of a system simulation computer run for 

a one-year period, which includes the solar energy availability statistics - 

normally weather bureau data for a selected year - and the expected load 

statistics. Ideal would be a simulation run over the system life to determine 

an f. , value which represents the average over the system life. However, 

Ld 

forward looking solar energy availability data do not exist, and even forward 
looking load statistics will be of doubtful validity. A compromise could be 
a backward looking simulation over a period equal in duration to the system 
life, using real data. The limited gain in confidence, however, generally does 
not justify the additional expense. A one-year run is usually felt necessary 
to properly include the seasonal changes and the short term meteorological 
variations . 


The total number of hours in the year (8760 h), 
multiplied by the load factor f^ represent in "equivalent time" t during 
which the array could have operated at peak power capability to produce the 
same amount of energy as actually delivered. It is additionally useful to 


define the quantity Pp^* the peak power output capability per unit area of the 


array, which is simply; 


p . =• II . • i|, |kW nf 2 | 

pK pk A r , s t d 1 1 


(IV. 6) 


IV - i 


The energy E^ delivered from the array directly to the load will 
generally be less than E q , being reduced by the power conditioning subsystem 
efficiency n pc » and by the fraction f St of the annual array output which is, 
In the average, transferred into the storage subsystem 

E Ld, dir ■ K o <» - f st> - W < lv - 7 > 


In addition, the energy St is delivered from the storage subsystem, to the 
load: 


"Ld,St 


E o f St 


'St 


n p(; ; fkWh’y -1 ] 


(IV. 8) 


where n Ct . i8 the efficiency of the storage subsystem. 

In the relationship of eq. (IV. 7), the assumption is made that all power conditioning 
occurs after storage. Otherwise, the efficiency n pc would have to be broken into 
several terms. 

Summing eq. (IV. 7) and (IV. 8) yields then the total energy E^ delivered 
to the load: 


E ld = E o 11 " f St ° " r| St )1 M PC ; IkW'-y" 1 I (IV. 9) 

The expression in the brackets, which is a function of the Joad curve relative 
to the solar energy availability curve, as well as of system design, incl. type 
and capacity of the storage device, could be represented by a "storage transfer 
factor" Tj, , so that eq. (IV. 9) can be written as: 

E Ld = E o * T St * n PC ;|kWh y_1) (IV ‘ 9a > 

Tt has to be noted that the load factor f , included in E o , is also dependent 
on the same variables as x , and generally increases with Increasing f and 

w t DC 
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The system power delivery capability 1’ which is usually limited by the 

sy 

power conditioning subsystem and/or the storage subsystem capacities* can be 
related through the factor fp Q to the array peak power capability: 


Sy 


A. P , 
Ar*pk 


St 


'Pc 


Po’ 


IkWj 


(IV. 10) 


The factor fp Q may be smaller or greater than unity. 

The storage subsystem capacity can illustratively be expressed by the 
"equivalent storage time" t^ t which is the time interval for which the storage 
device, when originally fully charged, could provide energy at the peak system 
power delivery rate, until discharged to a predetermined minimum charge state: 


E 


st 



t St ; [kWh] 


(IV.11) 


3. Evaluation of the Energy-Cost Effectiveness of Competing Subsystem Options 

The entire photovoltaic solar energy conversion system is composed of a 
network of subsystems, basically connected in series according to the energy 
flow, as indicated in Fig. IV. 1. The individual subsystems may be defined in 
any way which facilitates the system analysis or t ho cost determination. Tlius, 

« foundation for the solar array may be considered a subsystem, as a circuit 
breaker for system protection, or a battery for energy storage may be. Clearly, 
the entire system cost is the sum of all the individual subsystem costs C^: 

N 

C ca P “ 1 C i ;t$1 (IV. 12) 

i-1 
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and the system performance is a function of the performance of all the sub- 
systems. Frequently, some subsystems are not directly in the line of energy 
flow, as Indicated in Fig. IV. 1 by subsystems 3.1 to 3.M. For an evaluation as 
discussed here, it is best to combine these into a "subsystem group" which, as 
a whole, is in the line of energy flow. The cost of the subsystem group is 
then the sum of the costs of the subsystems within the group. The evaluation 
of the cost effectiveness of an individual subsystem of the group can be per- 
formed by expansion of the methodology outlined here*. 

In general, both cost and performance of n subsystem are the result 
of an engineering design trade-off in which the performance characteristics 
of available devices and their commercial prices are considered, as well as 
the subsystem complexity arc assembly cost. It is the purpose of this section 
to outline a methodology for assessing the cost effectiveness of such trade- 
offs from the viewpoint of the cost of the energy produced by the system. 

Since the entire system (Fig. IV. 1) can be viewed as a series connection 
of subsystems i or groups of subsystems, its efficiency can be expressed as 

the product of the efficiencies nj of the individual subsystems or groups of 
* subsystems: 


or 


N 

n syst “ n Ar,std * T St* n PC * f Ld = ^ ''V 


N 


n „ * f. , n n,; 

syst Ld . i’ 
i=*l 


(IV. 13) 


(IV. 13a) 
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Representation (IV. 13) !• a generalization to the subsystem level of the ex- 

preasion contained in eq.IV.9a where the n’* include all the contributions 

contained in the t'ficienciea n. ... n u „, and in the quasi-efficiencies 

Ar.stn rt. 

i St and f w . In the second version (eq. IV. 13a), the H* factors contain all 

the direct efficiency-like influences of each of the subsystems, while all 

indirect, or second-order influences are relegated i-. the "reduced load 

factor" f’,.. The application and practicality of this approach will be 
La 

recogniz'd later in this paper. 

These subsystem efficiencies have an impact on the dimensioning of the 
individual subsystems, and consequently on their costs, since the system has 
to be designed to satisfy a given load by supplying a certain E^. Thus, sub- 
systems placed nearer the beginning of the series connected subsystem chain 
have to be dimensioned relatively larger to account for the losses of subsequent 
subsystems. This principle is recognizable in eq. IV. 9a where the array output 
K is larger by the inverse of the product fpf; • than the ener-y which 

is delivered to the load. 

The division into subsystems can be practically pursued to the smallest, 
separately identifiable, functional urits will, their individual efficiencies. 

This shall be illustrated by example of the photovoltaic array, with its 
array efficiency n which is frequently considered as composed of 

AIT f S CO | 

"suharravs" which are made up of "modules". The module contains a group of 

.solar ceils (a subsystem) which have an average efficiency U In 

Ce, std 

series/parallel connecting these cells of slightly di'lering characteristics 
Into a "matrix" (a subsvst on) , a small loss in polonl i.il power output is 
incurred, expressed in the " matrixing ellleieney" p . The interconnect 
wiring in this matrix is another separately idem i I i ah I <• subsystem with its 
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Joule losses, which are accounted for in the "wiring efficiency*' n^. The 
encapsulation forms two functional subsystems. The first is the window, in- 
cluding adhesive or pottant, with its optical transmission losses, leading 
to the encapsulation efficiency The second performance influencing 

attribute of the encapsulation is the heat transfer to the environment which 
determines the operating temperature of the array which controls the instantaneous 
operating efficiency of the module. This effect produces an "average annual 
cooling effectiveness factor" f^ o , a quasi-efficiency which usually is included in the 
load factor f,,. At the subarray (subscript SA) and the array (subscript Ar) 
levels, matrixing and wiring losses are again incurred, so that the cell energy 
output will have to be: 


; ; [ kWhy” 1 1 (IV, 14a) 

Ce n Ma Hyi n En * n Ma,SA fi Wi,SA ’ ri Ma,Ar n Wi,Ar 


The installation of the subarrays forms another subsystem which influences 
system performance twofold: through the subarray orientation, which may include 

one-or-two-dimensional tracking, and through the cooling effectiveness. Both 
of these attributes form quaai-cf f icicncy factors which are part of the load 
factor f, ,. 


Since the orientation/tracking effect is a direct influence which can, 
under exclusion of variable atmospheric effects, bo analytically evaluated, it 
can be beneficial to eliminate this performance factor from the (reduced) load 
factor and attach it as an efficiency factor to tho installation (or tracking) 
subsystem. 

Formally applying these principles by combining eq.IV.S, IV. 9a and IV. 13, 
yields an expression for the energy delivered to the load, K , in terms of the 


subsystem efficiencies and quas i-ol I ic ienc io 


i ; i i : : 

i 
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Introducing this expression together with cq.IV.12 Into eq. IV. 3 |ivu the 
energy cost deterninntor T in terms of subsystem cost and performance data, and 
constants, only: 
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(IV. 16) 


Following the approach used by Rodfield in his "cost/Watt optimization" 
(2), the parameters of a single subsystem or subsystem group k of interest can 
be isolated in eq. IV. 16: 


or: 
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The expressions 7. and Z stand for the sum or product, respectively, 
l#k i*k 

over all values of i from l to N, except for tin- value k. ThlR form of r 
permits the evaluation of various design options for the same subsystem, or 
group of subsystems, with differing costs and efficiencies, with respect to 
their influence on the cost of the energy produced. Such an evaluation is 
particularly simple, if only and are variables of the design options. 

Then, a first order Taylor expansion yields: 


ah ■ b • 



(] + 


Y C, 

i*k 1 


) An 



t$kWh‘ 1 yl 


(IV. 18) 


where AC^ and Ar^ arc - positive or negative - differences *ea«nst the base 
case in subsystem cost and efficiency, respecti”"l v , which ret from the 
change in design of subsystem k. The constant I*, in cq.IV.18 is the product of 
the first two of the three terms on the right hand .side of eq.IV.17a. A negative 
AT indicates a reduction in energy cost, and consoi.u ^nlly a design improvement. 
It is readily apparent from rq. IV. 18 that cost reductions and efficiency 
decreases counteract rich other. 


The condition imposed for the drr ivat ion of oq . IV. 18, that only and n'^ 
are variables of the design options, is in apparent conflict with several state- 
ments made in the preced^rp discussion. Thus, load factor can he affected 

by changes in tfie system efficiency, particularly • changes in tfie storage 
transfer factor t , To make the evaluations tractable, it Is practical to 

i 

proceed iteratively by considering tfie reduced load factor f ^ as temporarily 
constant and re-evaluating it only after several changes in the efficiencies. 

Tills procedure illuminates the need for the definition of a "reduced load factor" 
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according to eq. IV. 13a which contains only second order effects of the 
efficiencies and quasi-efficiencies. The iteration is frequently speeded 
by reinforcing properties of the second order effects. For instance, efficiency 
improvements tend, at constant E^, to result in /i. creased load factors. 

The condition for the validity of eq.IV.18 further requires that 
and are independent of the designs of the other subsystems, and particularly 
that the design choice of subsystem k does not influence costs and efficiencies 
of the remaining subsystems. Th*s condition can, in principle, always be 
fulfilled by judicious choice of the designation "subsystem k", so that inter- 
dependent parts of the system are included in the same subsystem. 

A change in the efficiency of one subsystem affects, however, the system 
as a whole. While the resulting change in output energy E^ ^ is appropriately 
accounted for in T , one or more of the subsystems subsequent to the changed 
subsystem in the chain may now be over- or underdimensioned, and the load may 
no longer be supplied as desired. This problem requires considering the system 
of concern in somewhat more detail. 

The majority of the functional subsystems of a photovoltaic solar energy 
conversion system are basically modular and thus essentially without economics 
of scale, at least within the range of concern in an individual design trade- 
off study. The costs of these subsystems can therefore be expressed as a unit 
cost times a quantity factor. Such quantity factors are the array area 
the power handling capacity P of some subsystems, and, for some energy storage 
related subsystems, the energy capacity E. Generalizing the usage in ref (1) 
and (2), the system cost can then he expressed as: 
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cap 
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(IV. 19) 



iv ' 12 



The area-based unit costs C A k apply to the array related subsystems, Ln- 


f eluding its installation and land costs. The power-based unit costs Cp^ 

are connected with the power conditioning and other power handling equipment, 
although a part of the costs of the energy storage subsystems can also be 
proportional to power, for instance through the charge or discharge rates. 

The energy-based unit costs c_ are concentrated in the storage subsystem. 

The remaining costs, including the system— status sensors and the control logic, 
represent the "fixed" costs, Cp 

Using this expression IV. 19 for the capital costs in the energy cost 

determinator eq. IV. 16, and simultaneously extracting the iteratively constant 

reduced load factor f' , from the efficiency product H n. , yields the form: 

Ld j *- 
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[$ kWh _1 y ] (IV. 20) 


It will be observed that, in general, for every index i in the sum, 
only one of the unit cost factors c ., c ., , nr C . will be unequal 

to zero. An exception to this rule is known to exist in certain advanced 
storage batteries whose price is based on a combination of energy and power 
rating. Also, power conditioning subsystem groups may contain fixed cost sub- 
systems, such as control elements. 

In considering the quantity factors A^, P., and E. , several possible 
simplifications are immediately noticeable. First, all area based unit costs 
are commonly related either to the array area or to the solar cell area. The 
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latter is connected to the array area through the packing factor f pg < 1: 
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(IV. 21) 


Second, cq. IV. 11 relates to P Sy through the equivalent storage time, and 
thus permits combined treatment of the second and third terras of eq. IV. 20. 
Third, the dimensioning of each subsystem i of power dependent cost in the 
chain is determined by the system output specifications and the efficiencies 
of the subsystems subsequent to i in the direction of energy flow, so that: 



i = i+1 


(kW] 


(IV. 22) 


This permits expressing eq. IV. 20, under use of eq. IV. 6, IV. 10, IV. 11, IV. 13a 
and IV. 21, and under application of the subscript Ce to the cell area based costs, 
Ar to the array area based costs, as: 
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where t , is zero or t , depending on the existence of an energy based 

O t I 1 St 

cost contribution in subsystem i. The fixed costs, shown in the third term 
in the brackets of ccj. IV. 23, contribute to the energy costs independently of 
the subsystem's or the system's performance. They are also the only ones ex- 
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hibiting any direct economics of scale. 

The first term in the bracket of eq. IV. 23 can be evaluated for the 
impact of design options for an individual subsystem k in complete analogy 
to eq. IV. 17 to IV. 18. The second term, however, requires a slightly different 
treatment: 
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Consequently, first order Taylor expansion of eq. IV. 23 yields the total cost- 
effectiveness criterion A F, for a design change of subsystem k: 
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are the respective "investment per (unit energy per year)" ratios for all 
subsystems, except subsystem k, with area based unit costs, combined; for all 
subsystems preceding subsystem k in the chain, with power or energy based unit 
costs, combined; and for all subsystems, except subsystem k, with fixed sub- 
system costs, combined, correspondingly defined are the total subsystem in- 
vestments: 
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which represent the combined normalized costs of all subsystems, except sub- 
system k, with area based unit costs; of all subsystems preceeding subsystem 
k in the chain, with power or energy based unit costs; and of all subsystems, 
except subsystem k, with fixed subsystem costs; respectively. Examples of the 
application of eq, IV. 25 are shown in the section entitled: "Examples of 

Application of the Methodology." 

It is interesting to note that the three terms in the "cost effectiveness 
criterion" AF^ contain the "investment per (energy per year)" ratios for the 
remainder of the system, multiplied by the difference between two terms which 
are based on the relative cost change and the relative efficiency change, 
respectively. It is to be noted, however, that the relative cost change 
is based on the cost of the remainder of the system, while the relaHvp 
efficiency change is based on the efficiency of t lie subsystem under evaluation. 
The expression "remainder of the system" refers here to the subsystems with 
equally based unit costs, and, in the case of power or energy based unit costs, 
only to the subsystems preceding in the chain the subsystem being evaluated. 

For the "fixed cost subsystems", there is no efficiency influence. 

It may also be noted that the cost-effectiveness criterion (eq. IV. 25) 
contains the terms 
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where the relative efficiency change of subsystem k can refer to a subsystem 
of power based unit cost, but influence the cost-effectiveness through the 
subsystems of area based cost structure, or vice versa. The latter, inverse 
case is, however, not likely to occur as a subsystem of power based unit cost 
is rarely followed by a unit area cost based subsystem fn the photovoltaic 
power system chain. 
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Introducing this expression together with cq.IV.12 Into eq. IV. 3 |ivu the 
energy cost deterninntor T in terms of subsystem cost and performance data, and 
constants, only: 
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Following the approach used by Rodfield in his "cost/Watt optimization" 
(2), the parameters of a single subsystem or subsystem group k of interest can 
be isolated in eq. IV. 16: 


or: 
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The expressions 7. and Z stand for the sum or product, respectively, 
l#k i*k 

over all values of i from l to N, except for tin- value k. ThlR form of r 
permits the evaluation of various design options for the same subsystem, or 
group of subsystems, with differing costs and efficiencies, with respect to 
their influence on the cost of the energy produced. Such an evaluation is 
particularly simple, if only and are variables of the design options. 

Then, a first order Taylor expansion yields: 
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where AC^ and Ar^ arc - positive or negative - differences *ea«nst the base 
case in subsystem cost and efficiency, respecti”"l v , which ret from the 
change in design of subsystem k. The constant I*, in cq.IV.18 is the product of 
the first two of the three terms on the right hand .side of eq.IV.17a. A negative 
AT indicates a reduction in energy cost, and consoi.u ^nlly a design improvement. 
It is readily apparent from rq. IV. 18 that cost reductions and efficiency 
decreases counteract rich other. 


The condition imposed for the drr ivat ion of oq . IV. 18, that only and n'^ 
are variables of the design options, is in apparent conflict with several state- 
ments made in the preced^rp discussion. Thus, load factor can he affected 

by changes in tfie system efficiency, particularly • changes in tfie storage 
transfer factor t , To make the evaluations tractable, it Is practical to 

i 

proceed iteratively by considering tfie reduced load factor f ^ as temporarily 
constant and re-evaluating it only after several changes in the efficiencies. 

Tills procedure illuminates the need for the definition of a "reduced load factor" 
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according to eq. IV. 13a which contains only second order effects of the 
efficiencies and quasi-efficiencies. The iteration is frequently speeded 
by reinforcing properties of the second order effects. For instance, efficiency 
improvements tend, at constant E^, to result in /i. creased load factors. 

The condition for the validity of eq.IV.18 further requires that 
and are independent of the designs of the other subsystems, and particularly 
that the design choice of subsystem k does not influence costs and efficiencies 
of the remaining subsystems. Th*s condition can, in principle, always be 
fulfilled by judicious choice of the designation "subsystem k", so that inter- 
dependent parts of the system are included in the same subsystem. 

A change in the efficiency of one subsystem affects, however, the system 
as a whole. While the resulting change in output energy E^ ^ is appropriately 
accounted for in T , one or more of the subsystems subsequent to the changed 
subsystem in the chain may now be over- or underdimensioned, and the load may 
no longer be supplied as desired. This problem requires considering the system 
of concern in somewhat more detail. 

The majority of the functional subsystems of a photovoltaic solar energy 
conversion system are basically modular and thus essentially without economics 
of scale, at least within the range of concern in an individual design trade- 
off study. The costs of these subsystems can therefore be expressed as a unit 
cost times a quantity factor. Such quantity factors are the array area 
the power handling capacity P of some subsystems, and, for some energy storage 
related subsystems, the energy capacity E. Generalizing the usage in ref (1) 
and (2), the system cost can then he expressed as: 
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The area-based unit costs C A k apply to the array related subsystems, Ln- 


f eluding its installation and land costs. The power-based unit costs Cp^ 

are connected with the power conditioning and other power handling equipment, 
although a part of the costs of the energy storage subsystems can also be 
proportional to power, for instance through the charge or discharge rates. 

The energy-based unit costs c_ are concentrated in the storage subsystem. 

The remaining costs, including the system— status sensors and the control logic, 
represent the "fixed" costs, Cp 

Using this expression IV. 19 for the capital costs in the energy cost 

determinator eq. IV. 16, and simultaneously extracting the iteratively constant 

reduced load factor f' , from the efficiency product H n. , yields the form: 
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It will be observed that, in general, for every index i in the sum, 
only one of the unit cost factors c ., c ., , nr C . will be unequal 

to zero. An exception to this rule is known to exist in certain advanced 
storage batteries whose price is based on a combination of energy and power 
rating. Also, power conditioning subsystem groups may contain fixed cost sub- 
systems, such as control elements. 

In considering the quantity factors A^, P., and E. , several possible 
simplifications are immediately noticeable. First, all area based unit costs 
are commonly related either to the array area or to the solar cell area. The 
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latter is connected to the array area through the packing factor f pg < 1: 
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Second, cq. IV. 11 relates to P Sy through the equivalent storage time, and 
thus permits combined treatment of the second and third terras of eq. IV. 20. 
Third, the dimensioning of each subsystem i of power dependent cost in the 
chain is determined by the system output specifications and the efficiencies 
of the subsystems subsequent to i in the direction of energy flow, so that: 
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This permits expressing eq. IV. 20, under use of eq. IV. 6, IV. 10, IV. 11, IV. 13a 
and IV. 21, and under application of the subscript Ce to the cell area based costs, 
Ar to the array area based costs, as: 
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where t , is zero or t , depending on the existence of an energy based 
cost contribution in subsystem i. The fixed costs, shown in the third term 
in the brackets of ccj. IV. 23, contribute to the energy costs independently of 
the subsystem's or the system's performance. They are also the only ones ex- 
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hibiting any direct economics of scale. 


f The first term in the bracket of eq. IV. 23 can be evaluated for the 

I 

| impact of design options for an individual subsystem k in complete analogy 

i to eq. IV. 17 to IV. 18. The second term, however, requires a slightly different 
treatment: 
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Consequently, first order Taylor expansion of eq. IV. 23 yields the total cost- 
effectiveness criterion AF^ for a design change of subsystem k: 



where: 
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are the respective "investment per (unit energy per year)" ratios for all 
subsystems, except subsystem k, with area based unit costs, combined; for all 
subsystems preceding subsystem k in the chain, with power or energy based unit 
costs, combined; and for all subsystems, except subsystem k, with fixed sub- 
system costs, combined, correspondingly defined are the total subsystem in- 
vestments: 
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which represent the combined normalized costs of all subsystems, except sub- 
system k, with area based unit costs; of all subsystems preceeding subsystem 
k in the chain, with power or energy based unit costs; and of all subsystems, 
except subsystem k, with fixed subsystem costs; respectively. Examples of the 
application of eq, IV. 25 are shown in the section entitled: "Examples of 

Application of the Methodology." 

It is interesting to note that the three terms in the "cost effectiveness 
criterion" AF^ contain the "investment per (energy per year)" ratios for the 
remainder of the system, multiplied by the difference between two terms which 
are based on the relative cost change and the relative efficiency change, 
respectively. It is to be noted, however, that the relative cost change 
is based on the cost of the remainder of the system, while the relaHvp 
efficiency change is based on the efficiency of t lie subsystem under evaluation. 
The expression "remainder of the system" refers here to the subsystems with 
equally based unit costs, and, in the case of power or energy based unit costs, 
only to the subsystems preceding in the chain the subsystem being evaluated. 

For the "fixed cost subsystems", there is no efficiency influence. 

It may also be noted that the cost-effectiveness criterion (eq. IV. 25) 
contains the terms 
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where the relative efficiency change of subsystem k can refer to a subsystem 
of power based unit cost, but influence the cost-effectiveness through the 
subsystems of area based cost structure, or vice versa. The latter, inverse 
case is, however, not likely to occur as a subsystem of power based unit cost 
is rarely followed by a unit area cost based subsystem fn the photovoltaic 
power system chain. 
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The "cost effectiveness criterion" AI^ permits the evaluation of various 
subsystem design options both with respect to their benefit (or harm) in 
comparison to a baseline design, through the sign of and with respect 

to the relative merits of the different options, through the magnitude of 

Ar k* 

"Optimizations", that is a search for AI^ * 0 as discussed in ref. (2), will, 
with very few exceptions, not be possible, since the relationships between 
cost and performance are usually not available in functional form and, more- 
over, seem always to be limited by the contemporary, and often rapidly changing 
status of technology. "Relative evaluations", as discussed here, applied to 
specific subsystem design options, are, however, readily performed. 

The method is oasy to apply, since for the subsystem to be evaluated, 
only the cost and performance differences against a baseline design have to be 
known, and since the other needed inputs involve only a few summary data on the 
remainder of the system. While it may be, in some cases, difficult to obtain 
exact data for the remainder of the system, intelligent estimates will frequently 
suffice. When such estimates are used for the cost of the remainder of the 
system, error estimates should he made, as mls-est imat ion of the cost could 
shift the relative impact of the competing terms involving the subsystem cost - 
and efficiency - changes. 


4. Evaluation of Cost-Effectiveness of Manufacturing Process Options 

While many of the subsystems in a photovoltaic solar energy system are 
assembled of standard components by common methods, the solar cells, their 
assembly into arrays, and at a later time perhaps also the energy storage 
device, are specially manufactured items which represent a significant part 
of the total system cost. Since producing these devices with their highest 
possible performance at the lowest price is the fundamental condition for 
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success in large scale introduction of photovoltaic solar energy systems, 
comparative evaluations of the various available options for each step of 
the manufacturing process sequence need to be performed. A methodology very 
similar to that outlined for evaluation of the subsystem design options can 
be applied for this purpose. 

Evaluation methodologies for the solar cell and the module manufacturing 
processes are of greatest current interest. Both of these "subsystems" have 
an area based unit cost structure, and can therefore be treated by the same 
approach. The quantity to be reduced as far as possible is the "investment 
per (unit energy per year)" T (eq. IV. 23) which can be expressed as the sum 
of various sub-gammas for the different subsystems: 

N -1 
r - l (r A i + l p i + r Fji ); [$ kWh A yl (IV. 28) 
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where that for the subsystems of area based unit costs has the form: 
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As the solar cells and the modules are among the first subsystems in the chain, 

and are not preceded by power based subsystems, only the T . terms need to 

A, 1 

be considered for an evaluation of the manufacturing processes for one of these 
two subsystems. Thus, for the solar cells as subsystem k, it is: 
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using eq. IV. 26a for simplification. The product in the parenthesis of the 
first term is independent of subsystem k. 

The fabrication process sequence for subsystem k, the solar cells, shall 
be composed of I* process steps, with the individual step p costing c„ . on 
the basis of unit area of good work-in-process (partly processed solar cells) 
leaving the process station. The subsystem cost c Cg however, is based on 
the area of the finished, good cells leaving the end of the solar cell production 
line. Since each process step is afflicted with a certain yield y^, the amount 
of solar cell area to be processed through step p has to be increased above the 
finished cell area to make up for the yield losses of the subsequent process 
steps. Consequently, the unit cell area cost of the subsystem k can be expressed 
ns: 

P c 

c Ce,k B r - ~ p 6|k,P 5 [$ m " 21 (IV. 31) 

p-1 n y. 
n-p+l 

For solar cells, it has been long-standing practice^ to calculate an ideal- 
ized, theoretical "limit efficiency" Lim and to gauge the success in 
design and fabrication of the "real" solar cells by determining the various 
"loss factors" <{> which describe the degree of approach to ideality for the 
identified, efficiency influencing parameters. In variation of this practice, 
Redfield (2) assigned a loss factor to each of the process steps to facilitate 
his "cost/Watt" evaluation. Adapting this practice, the efficiency of the subsyntem 
k can be expressed as a limit efficiency times a product of loss factors. 

P 

n ■ n n l t l j 
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Each of the loss factors 4, is attributed to n step in the serial sequence 

k,p 

of process steps, and it expresses, by being norn.il lv loss than unltv, the 
degree to which the individual process step causes t he subsystem performance 
to deviate from Ideality. Different competing process options can usually be 
expected to cause different degrees of deviation from ideality. While for 
solar cells, a limit efficiency near 0.25 is usually discussed, for the module 
or panel assembly, a limit efficiency of unity will be practical to assume. 

Making use of eq. IV. 27a, IV. 29, IV. 31, and IV. 32 permits expressing eq. IV. 30 
in a form more conducive to derivation of the coat-effectlvenesa criterion: 
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In this form, the three characteristic attributes <! ,. _,Y_, and c„ , 

K,n* n* t.e,k,n 

of process step n which is the step to be evaluated, have been isolated. 

% 

Applying again a first order Taylor expansion to the investment per (energy 
per year) ratio, this time based on eq. IV. 28 and IV. 33, yields the cost 
effectiveness criterion AT for the individual solar cell manufacturing 
process step n : 
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where I. ... and C ... are used as before (eq. IV. 21a and IV. 27a), and where: 
A§ ifK A , IfR 
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expresses the fully yielded cost of the work-in-process required as input for 

step n In order to fabricate a unit area of output work-in-process from 

P 

this step. The factor fi y is the product o f the yields of the process steps 

^■n+i 

subsequent to step n. The inverse of this product gives the area of work-in- 
process to be processed through step n in order to obtain a unit area of fin- 
ished product (subsystem k). The application of eq.IV.34 is demonstrated on 
hand of an example in the next section. 

Similar to the subsystem cost-ef fectivonoss criterion AT , the manufacturing 

process cost-effectiveness criterion AT is the product of n variable 

R y 

factor and the "investment per (energy per year)" ratio for the remainder of 
the system, in this case, however, limited to the part of the system which is 
based on unit area costs. The variable factor contains three terms. The first 
describes the influence of the difference in cost t r of the subject 

» /i* i k 

process options against the baseline case, or against another option, taken 
relative to the total cost of all other subsystems of unit area based cost. 

The impact of this relative cost difference is magnilied by the inverse of 
the product of the yields of all process steps which follow the step under 
evaluation (n) in the process sequence up to the finished subsystem k. The 
second term describes the Impact <>f the iclative change in the yield of process 
step n which would he Incurred by switching to the option being evaluated. 

Tills relative yield change is multiplied by the cost of the input work-in- 
process to step n, divided by the total cost of all other subsystems of unit 
area based costs. Again, the impart of this term is increased through 
the yields of all subsequent process steps. The third term finally Is principally 
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the relative solar cell efficiency change resu 1 t i tig from Introduction of 
the subject process option. The impart of this relative efficiency change 
Is raised above unity by the ratio of the cost (per unit area) of the subsystem 
considered to the sum of the unit area costs of all other subsystems of area 
based cost structure. 

Examination of eq. IV. 34 shows that the knowledge of the "Investment per 
(energy per year)" ratio for the remainder of the system is not needed for 
comparative evaluation of different process options, as this ratio is n 
constant factor in the cost-effectiveness criterion. This leaves only four 
data required as constant inputs for the evaluation: the cost of the input 

work-in-process; the cost of the finished subsystem; the total coHt ot the 

remaining subsystems of area based costs; and the product of the yields of 
the subsequent process steps. The variable inputs are the relative changes in 
the three key attributes of the option for the process step to be evaluated: 
cost, yield, and efficiency coitrihution. Since exact data for the four 
constant inputs may be difficult to obtain, intelligent estimates will some- 
times be substituted. This procedure appears, at first look, appropriate as 
these quantities form constant multipliers. However, tills approach has to he 
applied with caution since significant mis-est i mat ion could shift the relative 
Impacts of the cost, yield, and efficiency terms. This caul ion will be necessary in 
the common cases, where the cost of the finished subsystem under evaluation 's 
small compared to the total cost of the remaining subsystems of area based 
coat, so that the multiplier on the relative efficiency change would not be 
large compared to unity. 

It is clear, that the method outlined here for the solar cell manufacturing 
process, and expressed in eq. IV. 34 , applies equally well to the array assembly 
processes, except for the omission of the packing factor f in that case, and 
the replacement of the subscript (e hv subscript Ar. 
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3. Examples of Applications of tha Methodology 

Two example* will demonstrate the application of eq. IV. 23 in evaluation 
of different design options for subsystem k. 

The subsystem under consideration shall be the solar cell. The base case 
is a cell with a conversion efficiency of 17.5% on the basis of the solar cell 
area. The following relevant data for the base case are known: 


Table IV. 1 



Item 

Symbol 

Data 

Units 

Basic 

1 . 

Solar cell price 

C Ce,k 

61.36 

$/m 2 

cell area 

2. 

Packing factor 

f Pg 

0.90 

— 

— 


Solar cell price 

c A,k 

55.24 

$/m 2 

module area 

1 

Module assembly 
add-on price 

C A,k+l 

23.50 

S/m 2 

module area 

1 

Foundation, array 
assembly, Installation, 
etc, add-on price 

C A,k+2 

50.00 

$/m 2 

module area 


Total area based costs 

C A 

128.74 


module area 

B 

Total area based coots 
except for subsystem k 

C A.iJ»k 

73.50 


module area 

D 

Module efficiency 


15.75 

■ 

module area 


t 

Problem 1. 

A process is anticipated by which the efficiency of the solar cells 
could be raised to 20%. How much more could the solar cells cost to pro- 
vide an at least equally cost-effective system? 
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Answer : 

a) Since the subsystem of concern is of area based coses only, the 
second and third terms of eq. IV. 25 are zero. 

b) The subsystem k contains only cell-area based costs, designated by 
subscript Ce, and no array-area based costs, designated by sub- 
scripts Ar. Thus: 


Ac 


Ar,k 


0 


c Ar,k 


0 


c) 


Since the packing factor fp^ does not change with the change of 
cell efficiency: 


A(f_ • c _ ) ■ f_ ■ Ac. , . 

Pg Ce,k Pg Ce, k 


In this case, also, It is immaterial cither module efficiences or 
cell efficiencies are used, as they are related through a constant, 
proportionality factor. 


d) Wanted is knowledge of Ac^ ^ for 


A!'. 


0 . 


A,i*k 


Transforming eq. IV. 25, after applying points a) to c) above; yields 
then: 
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e) The efficiency difference An^ going from the base case to the 
new subsystem option is 2.5%. All other numbers entering into 
the relationship given in point d) relate to the base case. Thus: 


» < +0.025 ,73.50 . ,, 

A W “0375 ( ~oX 


< 2 
* + 20.44 $/m cell area. 


A 14% cel; efficiency increase thus justifies a 33% cell cost increase 
for equal energy cost effectiveness, and any lower cost increase 
yields a more cost-effective system. 


The maximum price is thus: 


Base price: 


61.38 $/m 2 


cell area 


Maximum increase + 


Apply f pg = 0.90: 

Module add-on cost: 
Module cost 


2 

20.44 $/m cell area 
81.82 $/m 2 cell area 

2 

73.64 $/m module area 

2 

23.50 $/m module area 
97.14 $/m^ module area 


At 


180 W . / m 
pk 


output, this corresponds to 0.54 $/W 


pk' 


Problem 2 

In lieu of Czochralski grown wafers assumed to be used in the base case 
given above, the use of ribbon silicon is anticipated, resulting in a reduced 
cell efficiency of 14%, but an increased packing factor of 0.92. How much 
lower would the cell cost have to be to provide an at least equally cost 
effective system? 
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Answer: 


a) Points a) and b) of answer ] still apply. 

b) As the packing factor changes, 


% 


A(f_ • c_ , ) * f_ Ac_ , + c„ • Af_ 
Pg Ce,k Pg Ce,k Ce,k Pg 


will have to be used. 


c) Because of the change of packing factor, and since the energy 
cost determination is ultimately based on the array (or module) 
area related costs and efficiencies, the evaluation will have to 
use these latter efficiencies. For the base case, the module 
efficiency was 15.75%. For the option, it is 14*0.92 = 12.88%. 
Thus, An k = 2.87%. 

d) Under consideration of points 2a) and 2b) above, and solving 
for 



= 0 , 


as in Answer 1, eq. (25) transforms into: 
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e) The difference in packing factor is +0.02, compared to the base 
case. Outside of the efficiency difierenee, only data from the 
base case are needed: 
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2 

Ac Ce k “ “ 26,07 " 1,36 = ' 27,Ai cell area 

The maximum cell price for equal cost effectiveness is thus: 

61.38 $/m^ cell area 

-27.43 $/m 2 cell area 
33.95 $/m 2 cell area 

and the corresponding module price: 

Cells: 33.95 $/m^ * 0.92 = 31.23 $/m 2 module area 

2 

Module add-on cost +23.50 $/m module area 

54.73 $/m 2 module area 

2 

At 128.8 Wp^/m output, this corresponds to $0,425/1^^ for the 
module. 


Checks to Problems 1 and 2: 

Try 100 system: 

Base case: 

Area needed: 10"* W , : 157.5 W /m 2 = 632.9 m^ 

pk pk 

Module price: 0.50 -► 50,000 $ 

Installation etc.: 50 $/m 2 -* 31,64 5 $ 

Total 81,645 $ 


Option 1: 

Module efficiency: 18% 

Area needed: 10^ W : 180 W , /m* - = 555.6 m^ 

pk pk 

Module price 0.54 $/W , - 54,000 $ 

pk 

2 

Installtion cost etc. 50 $/m -*■ 27 , 780 $ 

81,780 $ 
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Option 2: 


Module efficiency 12.887' 
Area needed: 

Module price: 0.425 $/W 


10 5 W , : 128.8 W , /m 2 - 776.4 m 2 

pk pk 


pk 


42,500 $ 


Installation cost etc. $50/m > 38,820 $ 

81,320 $ 


Problem 3 

A process sequence for solar cell fabrication has been proposed by 
Motorola for 1986, which includes two diffusions for pn-junction and BSF 
layer formation. Starting with a texture-etched, cleaned wafer, a total 
of 5 process steps (spin-on silica front; BC 1 ^ diffusion; spin-on silica 
back, PH-j diffusion; strip oxide both surfaces) 's needed to produce a clean 
wafer ready for the next process step (AR coating). 

RCA has proposed a completely different process sequence for cell 
fabrication for 1986 which includes ion implantation for both pn-junction 
and BSF layer formation. The conditions of the wafer before and after the 
2-step process (ion- implantation, activation anneal) are equivalent to those 
before and after the 5-step Motorola diffusion process, except for possible 
differences in efficiency resulting from the two processes. Since the 
Motorola overall process sequence seems to he the less costly one, it will 
be used as the base case. Thus, in lieu of the diffusion process, ion im- 
plantation could be inserted into the base case process sequence. 


Question : 

One would like to know the relative cost eifect iveness of the 2 competing 
process options for pn-junction and BSF layer formation. 
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Answer: 


The costs and yields for the 2 process options are known, as well as 
tlit* costs and yields for all the other solar cell manufacturing process steps 
in the base case. The cost data from the 2 companies have been normalized 
to the same economic base through application of the SAM1CS standardized 
cost structure. No information is, however, available on the efficiency 
contributions of the 2 options. The evaluation will therefore be carried 
out by determining the efficiency difference which would make the 2 options 
equally cost-effective. Equation XV. 34 is therefore to be solved for 
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The information displayed in Table IV. 2 is available for the base process: 
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Based on the unit area of vork-in process leaving the respective process step. 

Based on the- unit area of vork-in-process leaving group of process steps, or sub-process. 

Based on rhe unit area of finished product. 


Table IV. 2 contains all the information needed for solving eq. IV. 38, which 




is summarised in Table IV. 3. 


Table IV. 3 



From Table II 
Column Line 

1 

Base 1 

Case Option 


C Ce,k,n 

5 10 

5 7a 

i 

9.40 | 

j 9.86 

^ c Ce,k,n 

— 

1 +0.46 

i 

y n 

6 6 

6 6a 

0.958 j 

j 0.980 

A 

yn 

— 

1 +0.022 

1 

C Ce,k 

C Ce,k,WPn 

P 

n y £ 

£**n+l 

^A.ij^k 

f Pg 

8 

. 5 5 

(divided by yield 
shown in column 6, 
Line 6) 

6 10 

from Probl. 1 
from Probl. 1 

73.95 | 

43.41 1 

1 

0.905 1 

1 

73.50 1 

1 

0.90 | 

1 


Thus: 


A< ^k , n 

^'k.n 


[0.46 


0.02 2 

0.958 


43.41 j 


(L9 

(73.5 +0.9 • 73.95) .0.905 
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- 0.0038 




Result: 


The RCA ion implantation process option thus could have an efficiency 
contribution 0.4% lower than that of the Motorola diffusion option, to 
achieve equal cost effectiveness in energy generation. The ion implantation 
process would thus, at equal efficiency contributions, be very slightly more 
cost-effective than the diffusion option, but the difference is so small 
that the two options really ought to be considered as equivalent. 

It may also be noted that experimental results obtained at various 
laboratories indicate that the expectation of equal efficiency contributions 
from the two process options considered is justified. Thus, the result of 
economic equivalence of the two particular options analyzed is realistic, 
as far as the projections to 1986 for the various cost contributions and yields 
can be considered realistic. 

Check : 

Since the efficiency contributions are considered equal for the two 
competing processes, the check can be performed on the cost and yield basis 
alone. 
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Table IV. 4 



Base Case 

1 

Option 

Units 

Input work in process 
on unit area basis 

41.59 

"T~ 

1 

41.59 

$/« 2 

Yield in process step 

95.8 

1 

98.0 

% 

Needed input work-in-process 
for unit output work-in- 
process 

1.044 

1 

1 

1.02 

m 2 /m 2 


«•— — — 

4* 

— — - 

- — — - 

Cost of input work-in- 
process 

43.41 

1 

42.44 

$/m 2 

Cost of process step per 
unit output work-in-process 

9.44 

1 

1 

9.86 

$/« 2 

Cost of output work-in-process 

52.85 

1 

_L 

52.30 

$/«" 


The option output work-in-process is thus 1% less costly. With its 0.4% 
lower permitted efficiency contribution, it becomes equivalent at the sys- 
tem level. At exact efficiency equality, it would be the (slightly) pre- 
ferable process. 
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6. Conclusions 


0< r&4t p 


A quantitative comparative evaluation is frequently needed of the different 
design options for a particular subsystem in a photovoltaic solar energy con- 
version system, or of the different options for a process step in the manu- 
facturing process sequence for such a subsystem. Such an evaluation lias to be 
functional, which means, based on the cost of the electrical energy produced 
by such a system. 


It is seen that such evaluations can he rather easily performed on the basis 
of knowledge of the quantitative differences of the key attributes of the 
particular option under consideration for a subsystem or a process step against 
the attributes of a baseline case or of a different option. The key attributes 
are cost and efficiency for the subsystem, assuming reliability and service life 
to be comparable, and cost, yield, and efficiency contribution for the process 
step. The other needed inputs are relatively few and of a rather fundamental 
nature, such as the investment needed for the whole system per unit of energy 
delivered annually, the total cost of the system exelusive of the subsystem 
being evaluated; or the cost of the input work-in-process to the particular 
process step being evaluated. In many instances, adequate evaluations can he 
performed by substituting estimated values for real data of these quantities. 


It is also noteworthy that, particularly for the manufacturing process 
step evaluation, an analysis on the "cost per peak Watt" basis will often be 
adequate as a first order approximation, since the load factor which is the 
principal variable in the conversion to the "cost per kWh delivered" basis, is 
•affected by the evaluation variables only through second order influences. 
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V. Energy Analysis 


1, Introduction 

Energy consumption and corresponding payback times were 
reviewed for the current production process and for pro- 
cesses which may be used on manufacturing lines in 1982 
and 1986. For a proper assessment of the payback times, 
the entire manufacturing process sequence from the mining of the 
ore, in this case quartzite, and its reduction to Si, up to 
the completion of the system, fully installed and ready for 
operation, should be analyzed. So far, we have only accumulated 
energy data up through module manufacture. Even for this part of 
the process sequence, we have so far performed only detailed energy 
analyses for those processes for which we have completed thourough 
technology c.nd cost studies. These processes include the 
Si0 2 reduction in the arc furnace, Cz crystal pulling, slicing, 
junction formation, and the energy content of encapsulation 
materials. For the remaining parts of the process sequence, we 
have inserted data from other studies, such as lies' 1974 com- 
pilation of the solar cell plant energy consumption^, and 
Bickler's design data for a $2/W(pk) and a $0 . 50/W(pk) Straw- 
man process. 

These energy consumption studies are summarized into 
two quantities of interest for energy source systems: 1. 

the total energy consumed in creating the energy source, and 
2. the relation of this "invested" energy to the useful energy 
output of the system Which, for solar energy utilization sys- 
tems, is often called the "payback time". While the former 
is a reasonably well defined quantity, which depends primarily 
on the methods used for creating the energy source, the pay- 
back time depends also on the use of the system and is 
thus less uniquely defined. 
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2. Energy Payback Time* 

To calculate payback times for photovoltaic systems, 
the annual energy output of a unit module area must be known. 

This power output, for photovoltaic systems, is 
dependent upon the encapsulated cell efficiency, the module 
packing factor, the solar insolation, the efficiency of 
the power conditioning and storage subsystems, and the mismatch 
of the energy availability and demand statistics. This list of 
dependencies makes it clear that the output of the module will vary 
according to climate at the locality of installation, and to the indi- 
vidual load to be satisfied. The factors of influence on the useful 
system output, outside of the cell efficiency, are combined into 
the "capacity factor". As a reasonably representative, not 
too optimistic value for this capacity factor, the number 0.11 

has been used, in consequence of the results of several system 
(4) 

studies. The capacity factor is essentially independent 
of the solar cell or module manufacturing process. Consequently, 
the energy payback time is only partly a function 
of the energy consumption for the solar module fabrication pro- 
cess, since the other system parameters strongly influence 
its absolute magnitude. In the calculations of energy pay- 
bark times, encapsulated cell efficiencies n of 11.5% and 
15% were employed for 1978 and 1982 and for 19Rf, , respectively, 

assuming the use of FFG ribbons, of 12%. To oLt.iin the 
energy payback times, the annual module output was calculated 

as E <■ n x 8766 x 0.11xHp k , H pk being the standardized peak 
solar irradiance, used as lkw(pk)/m 2 . The factor 0.11 repre- 
sents the "capacity factor". 
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Thus , annual energy outputs from a square meter of solar cells 
of 110.9, 144.6, and 115.7 kWh/m 2 , y are obtained for 1978, 

1982, and 1986 respectively. 

It is also to be observed that, since the energy recovery 
from solar cells is in the form of electrical energy, it is 
appropriate to convert those energy expenditures which occur 
in the form of heat of combustion, to equivalent electrical 
energy by applying the average efficiency of 0.30 experienced 
by the electric utilities in the conversion from heat of 

combustion of fuels to electric power delivered to the con- 

( 5 ) 
suroer' 
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i. Energy Consumption in Photovoltaic Solar Array 
Manufacturing Process Sequences 
3.1 Data Sources. 

The major sources of information for this energy analysis 
were data accumulated from LSA project reports and industry . 
interviews. The energy consumed through materials, both 
direct and indirect materials, was obtained by converting 
the material consumption to energy units through multi- 
plication with the material energy contents shown in Table V-l. 

Where specific information to the contrary was not available, the material 
energy values were assumed to be in the form of thermal energy. The 
equipment costs, were converted to energy units expressed on the 
basis of unit cell area of throughput, by assuming a lifetime 
of seven years, and an energy content of the equipment which corres- 
ponds in value to 2% of the equipment cost.^ 6 ^ This energy cost 

has then been converted to a thermal energy using an energy 

17 ) 

price of $0 . 003/kWh (th) ' . Finally, the general energy 

usage for operating the facility was derived from the machine floor 
area by using the annual SAMICS utility cost of $3.74/ft^ and assu- 
ming that all of this "indirect" energy, since it is used primarily for 
lighting , air conditioning, and ventilating, in view of the 

high heat load in the building, is in the form of electricity 

( 8 ) 

at a cost of $0. 0319/kWh . The conversion factors for 

the equipment and facility energies thus become 0.9523 kwh(th)/$y 

and 1262 kWh/m 2 y, respectively. 
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the energy payback time is only partly a function 
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cess, since the other system parameters strongly influence 
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Thus , annual energy outputs from a square meter of solar cells 
of 110.9, 144.6, and 115.7 kWh/m 2 , y are obtained for 1978, 

1982, and 1986 respectively. 

It is also to be observed that, since the energy recovery 
from solar cells is in the form of electrical energy, it is 
appropriate to convert those energy expenditures which occur 
in the form of heat of combustion, to equivalent electrical 
energy by applying the average efficiency of 0.30 experienced 
by the electric utilities in the conversion from heat of 
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( 5 ) 
suroer' 
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The major sources of information for this energy analysis 
were data accumulated from LSA project reports and industry . 
interviews. The energy consumed through materials, both 
direct and indirect materials, was obtained by converting 
the material consumption to energy units through multi- 
plication with the material energy contents shown in Table V-l. 
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Table V-l 

Energy Contents of Selected Materials 




Material 

Original Units 

Convenient Units 
(Thermal kWh's! 

Reference 

1. Acetic acid 

3.28 lb o' hi-press. 
ateaa jIus 
0.07 Wh(*v)/lb 

3.18 kWh/g 

1 

2. Aluminum (At) 

520x10® kWh(th)/ton 

0.17 kWh/g 

2 

3. Ammonia gaa (NH^) 

8.0Sxl0 3 Btu/lb 

0.00534 kWh/g 

3 

A. A— onlum Hydroxide 

• 

1.32 kWh/i 

4 

30% (hh 4 oh> 




S. Argon gas 

1,100 Btu/lb 

1.27xl0 -3 kWh/l 

5 

6. Butyl acetate 

4.32 lb of low-press, 
steam plus 
0.082 kWh(ei)/lb 

3.21 kWh/g 

6 

7. Copper (Cu> 

16.2x10* Btu/ton 

5.23 kWh/g 

7 

8. Energy for exhausting 
waste fuses without 

0.46 kW/1000 ft 3 /sin 

2.56x10" 2 kWh/ 
1000 ft 3 

8 

scrubbing. 




9. Freon-14 gas (CF ft ) 


2xlO -3 kWh/f 

9 

10. Hydrogen gas 

43,300 Btu/lb 

2. 51x1 0~ 3 kWh/f 

5 

11. HF (48X) 

7,000 Btu/lb 

5.22 kWh/f 

5 

12. Nitrogen (liquid) 

• 1,330 Btu/lb 

0.69 KWh/f. 

5 

13. Nitrogen (gas from 

- 

1.44xl0~ 3 kWh/i 

10 

liquid) 




14. Nitric Acid 

14,500 Btu/lb 

13.12 kWh H 

5 

(HN0 3 , 67X) 




15. Oxygen gas (0 2 ) 

830 Btu/lb 

7.66 kWh/* 

5 

16. Phosphorus (solid) 

23,790 Btu/lb 

1.54xl0" 2 kWh/g 

5 

17. Phosphine gas (PH^) 

- 

0.18 kWh/l 

11 

18. Phosphorous 


0.14 kWh/g 

12 

oxychloride (POCij) 

- 

19. Plating resist 

8,000 Btu/lb 

5.20 kWh/i 

13 

j 20. Silver 

1260xl0 6 Btu/ton 

0.406 kWh/g 

8 

21. Sodium Hydroxide 

by-product 

0 

5 

(NaOH) 




22. Toluene 

0.05 lb of low-press, 
steam plus 
0.025 kWh(ei)/lb 

0.0349 kWh/nf 

14 

23. Vacuua pump 
oil 

3,000 kWh(th)/ 
barrel 

18.87 kWh/* 

i 

15 

24. Wax 

2,000 kWh/barrel 

1.3xl0 -2 kWh/g 

16 
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has been used, in consequence of the results of several system 
(4) 

studies. The capacity factor is essentially independent 
of the solar cell or module manufacturing process. Consequently, 
the energy payback time is only partly a function 
of the energy consumption for the solar module fabrication pro- 
cess, since the other system parameters strongly influence 
its absolute magnitude. In the calculations of energy pay- 
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sents the "capacity factor". 

V-2 



Thus , annual energy outputs from a square meter of solar cells 
of 110.9, 144.6, and 115.7 kWh/m 2 , y are obtained for 1978, 

1982, and 1986 respectively. 

It is also to be observed that, since the energy recovery 
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combustion of fuels to electric power delivered to the con- 

( 5 ) 
suroer' 
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3. Energy Consumption in Photovoltaic Solar Array 
Manufacturing Process Sequences 
3.1 Data Sources. 

The major sources of information for this energy analysis 
were data accumulated from LSA project reports and industry . 
interviews. The energy consumed through materials, both 
direct and indirect materials, was obtained by converting 
the material consumption to energy units through multi- 
plication with the material energy contents shown in Table V-l. 

Where specific information to the contrary was not available, the material 
energy values were assumed to be in the form of thermal energy. The 
equipment costs, were converted to energy units expressed on the 
basis of unit cell area of throughput, by assuming a lifetime 
of seven years, and an energy content of the equipment which corres- 
ponds in value to 2% of the equipment cost.^ 6 ^ This energy cost 

has then been converted to a thermal energy using an energy 

17 ) 

price of $0 . 003/kWh (th) ' . Finally, the general energy 

usage for operating the facility was derived from the machine floor 
area by using the annual SAMICS utility cost of $3.74/ft^ and assu- 
ming that all of this "indirect" energy, since it is used primarily for 
lighting , air conditioning, and ventilating, in view of the 

high heat load in the building, is in the form of electricity 

( 8 ) 

at a cost of $0. 0319/kWh . The conversion factors for 

the equipment and facility energies thus become 0.9523 kwh(th)/$y 

and 1262 kWh/m 2 y, respectively. 
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Table V-l 

Energy Contents of Selected Materials 




Material 

Original Units 

Convenient Units 
(Thermal kWh's! 

Reference 

1. Acetic acid 

3.28 lb o' hi-press. 
ateaa jIus 
0.07 Wh(*v)/lb 

3.18 kWh/g 

1 

2. Aluminum (At) 

520x10® kWh(th)/ton 

0.17 kWh/g 

2 

3. Ammonia gaa (NH^) 

8.0Sxl0 3 Btu/lb 

0.00534 kWh/g 

3 

A. A— onlum Hydroxide 

• 

1.32 kWh/i 

4 

30% (hh 4 oh> 




S. Argon gas 

1,100 Btu/lb 

1.27xl0 -3 kWh/l 

5 

6. Butyl acetate 

4.32 lb of low-press, 
steam plus 
0.082 kWh(ei)/lb 

3.21 kWh/g 

6 

7. Copper (Cu> 

16.2x10* Btu/ton 

5.23 kWh/g 

7 

8. Energy for exhausting 
waste fuses without 

0.46 kW/1000 ft 3 /sin 

2.56x10" 2 kWh/ 
1000 ft 3 

8 

scrubbing. 




9. Freon-14 gas (CF ft ) 


2xlO -3 kWh/f 

9 

10. Hydrogen gas 

43,300 Btu/lb 

2. 51x1 0~ 3 kWh/f 

5 

11. HF (48X) 

7,000 Btu/lb 

5.22 kWh/f 

5 

12. Nitrogen (liquid) 

• 1,330 Btu/lb 

0.69 KWh/f. 

5 

13. Nitrogen (gas from 

- 

1.44xl0~ 3 kWh/i 

10 

liquid) 




14. Nitric Acid 

14,500 Btu/lb 

13.12 kWh H 

5 

(HN0 3 , 67X) 




15. Oxygen gas (0 2 ) 

830 Btu/lb 

7.66 kWh/* 

5 

16. Phosphorus (solid) 

23,790 Btu/lb 

1.54xl0" 2 kWh/g 

5 

17. Phosphine gas (PH^) 

- 

0.18 kWh/l 

11 

18. Phosphorous 


0.14 kWh/g 

12 

oxychloride (POCij) 

- 

19. Plating resist 

8,000 Btu/lb 

5.20 kWh/i 

13 

j 20. Silver 

1260xl0 6 Btu/ton 

0.406 kWh/g 

8 

21. Sodium Hydroxide 

by-product 

0 

5 

(NaOH) 




22. Toluene 

0.05 lb of low-press, 
steam plus 
0.025 kWh(ei)/lb 

0.0349 kWh/nf 

14 

23. Vacuua pump 
oil 

3,000 kWh(th)/ 
barrel 

18.87 kWh/* 

i 

15 

24. Wax 

2,000 kWh/barrel 

1.3xl0 -2 kWh/g 

16 
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References for Table V-l 


1. Assumed to be the same as that of acetone, which was 
taken from "Battelle Columbus Laboratories, Draft Target 
and Support Document for Developing a Maximum Energy 
Efficient Improvement Target for SIC 28',' Federal Energy 
Administration, Washington, DC (1976). 

2. J.T. Reding and B.P. Shepherd, "Energy Consumption, 

Report EPA-G5Q/2-75-032bV US Environmental Protection 
Agency, Washington, DC (4/75). 

3. Federal Energy Administration, "Project Independence 
Blueprint," Vol. 3, Federal Energy Administration, 

Washington, DC (1974). 

4. Calculated from the energy content for NH3 and using 
a density of 0.824 g/mil for the NH4GH solution. 

5. Battelle Columbus Laboratories, "Draft Target and Support 
Document for Developing a Maximum Energy Efficient 
Improvement Target for SIC 28',' Federal Energy Administra- 
tion, Washington, DC (1976). 

6. Assumed to be the same as that of butyl alcohol, as given in Ref. 5. 

7. H.W. Lownie, et al (Battelle Columbus Laboratories), 

"Draft Target Report on Development and Establishment of 
Energy Efficiency Improvement Targets for Primary Metal 
Industries," Federal Energy Administration, Washington, 

DC (9/76). 

8. M.G. Coleman, et al.. Motorola Final Report, DOE /JPL- 
954847-78/4, 183(11/78). 

9. Estimated as approximately 50% more than the average 
energy content of ccntnercial gases, as given in Ref. 5. 

10. Calculated from the energy content for LN2 assuming 
480 l of gas can be obtained from 18. of LN2 . 

11. Estimated as twice the energy content for P(s) plus 
3/2H 2 (g). 

12. Estimated as twice the energy content of P(s), plus ^02^)» 
plus |ci2(g)- The energy content of CI2 was 

taken from ref. 5. 

13. Assumed to be the same as that of wood rosin, as given in Ref. 5. 

1 

14. Taken as that of benzene, as given in ref. 5. 
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15. Taken as approximately twice the energy content of mid 
die oil distillates. 

16. From M. Sittig, Practical Techniques for Saving Energy 
in the Chemical Petroleum and Metals Industries, 

Noyes Data Corporation, Park Ridge, NJ (1977), and 
using a density of 0.97 g/ml for the wax. 


V-7 


3.2 Energy Consumption in Si reduction and purification. 

This process group starts with the reduction of Si0 2 in an arc 
furnace . This is a quite efficient and cost-effective pro- 
cess which is not very likely to be replaced by another 
approach. In contrast, the following step of Si purification 
is very inefficient and likely to be replaced by one of 
several alternate methods under development. Furthest advanced 

among these is the SiH purification process being prepared 

4 

for pilot line operations by Union Carbide Corp. It can be 

expected to be a production process by 1986. 

a) Theoretical Material Balance for the Arc Furnace 
Process. 


The chemical reaction of this process is: 


Atomic Weights 


Masses, normal- 
ized to 1 kg Si 
out: 


= f 28 + 

i { 2 . 


Si0 2 + 2C 


2(16) + 2(12) 
60 + 24 


Si + 2 CO 

28 + 2(12 + 16) 
28 + 56 


14 kg + 0.86 kg -*■ 1 kg + 2 kg 


According to Dow-Corning O) . r the industrially experienced 
and actual conversion efficiency of SiO^ to Si is 80%. 

The required input is thus: 2.68 kg Si0 2 for 1 kg Si out 

(Dow-Corning shows 2.71 kg Si0 2 ) . 

It is assumed that half of the lost SiO^ input, or 0.2 7 kg, 
is used in the reaction: 

Si0 2 + C -*■ Si*"* + CO 

Atomic Weights: 60 + 12 -*■ 44 + 28 

0.27kg + 0.05 kg * 0.2 kg + 0.12 kg 


where the SiO is lost at the top of the furnace. 

Thus, the total theoretical carbon input would be 0.91 kg C 
per 1 kg MG-Si output. 
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2. Energy Payback Time* 

To calculate payback times for photovoltaic systems, 
the annual energy output of a unit module area must be known. 

This power output, for photovoltaic systems, is 
dependent upon the encapsulated cell efficiency, the module 
packing factor, the solar insolation, the efficiency of 
the power conditioning and storage subsystems, and the mismatch 
of the energy availability and demand statistics. This list of 
dependencies makes it clear that the output of the module will vary 
according to climate at the locality of installation, and to the indi- 
vidual load to be satisfied. The factors of influence on the useful 
system output, outside of the cell efficiency, are combined into 
the "capacity factor". As a reasonably representative, not 
too optimistic value for this capacity factor, the number 0.11 

has been used, in consequence of the results of several system 
(4) 

studies. The capacity factor is essentially independent 
of the solar cell or module manufacturing process. Consequently, 
the energy payback time is only partly a function 
of the energy consumption for the solar module fabrication pro- 
cess, since the other system parameters strongly influence 
its absolute magnitude. In the calculations of energy pay- 
bark times, encapsulated cell efficiencies n of 11.5% and 
15% were employed for 1978 and 1982 and for 19Rf, , respectively, 

assuming the use of FFG ribbons, of 12%. To oLt.iin the 
energy payback times, the annual module output was calculated 

as E <■ n x 8766 x 0.11xHp k , H pk being the standardized peak 
solar irradiance, used as lkw(pk)/m 2 . The factor 0.11 repre- 
sents the "capacity factor". 
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Thus , annual energy outputs from a square meter of solar cells 
of 110.9, 144.6, and 115.7 kWh/m 2 , y are obtained for 1978, 

1982, and 1986 respectively. 

It is also to be observed that, since the energy recovery 
from solar cells is in the form of electrical energy, it is 
appropriate to convert those energy expenditures which occur 
in the form of heat of combustion, to equivalent electrical 
energy by applying the average efficiency of 0.30 experienced 
by the electric utilities in the conversion from heat of 

combustion of fuels to electric power delivered to the con- 

( 5 ) 
suroer' 
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3. Energy Consumption in Photovoltaic Solar Array 
Manufacturing Process Sequences 
3.1 Data Sources. 

The major sources of information for this energy analysis 
were data accumulated from LSA project reports and industry . 
interviews. The energy consumed through materials, both 
direct and indirect materials, was obtained by converting 
the material consumption to energy units through multi- 
plication with the material energy contents shown in Table V-l. 

Where specific information to the contrary was not available, the material 
energy values were assumed to be in the form of thermal energy. The 
equipment costs, were converted to energy units expressed on the 
basis of unit cell area of throughput, by assuming a lifetime 
of seven years, and an energy content of the equipment which corres- 
ponds in value to 2% of the equipment cost.^ 6 ^ This energy cost 

has then been converted to a thermal energy using an energy 

17 ) 

price of $0 . 003/kWh (th) ' . Finally, the general energy 

usage for operating the facility was derived from the machine floor 
area by using the annual SAMICS utility cost of $3.74/ft^ and assu- 
ming that all of this "indirect" energy, since it is used primarily for 
lighting , air conditioning, and ventilating, in view of the 

high heat load in the building, is in the form of electricity 

( 8 ) 

at a cost of $0. 0319/kWh . The conversion factors for 

the equipment and facility energies thus become 0.9523 kwh(th)/$y 

and 1262 kWh/m 2 y, respectively. 
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Table V-l 

Energy Contents of Selected Materials 




Material 

Original Units 

Convenient Units 
(Thermal kWh's! 

Reference 

1. Acetic acid 

3.28 lb o' hi-press. 
ateaa jIus 
0.07 Wh(*v)/lb 

3.18 kWh/g 

1 

2. Aluminum (At) 

520x10® kWh(th)/ton 

0.17 kWh/g 
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3. Ammonia gaa (NH^) 

8.0Sxl0 3 Btu/lb 

0.00534 kWh/g 

3 

A. A— onlum Hydroxide 
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1.32 kWh/i 
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30% (hh 4 oh> 




S. Argon gas 

1,100 Btu/lb 

1.27xl0 -3 kWh/l 

5 

6. Butyl acetate 

4.32 lb of low-press, 
steam plus 
0.082 kWh(ei)/lb 

3.21 kWh/g 
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7. Copper (Cu> 

16.2x10* Btu/ton 

5.23 kWh/g 

7 

8. Energy for exhausting 
waste fuses without 

0.46 kW/1000 ft 3 /sin 

2.56x10" 2 kWh/ 
1000 ft 3 

8 

scrubbing. 




9. Freon-14 gas (CF ft ) 


2xlO -3 kWh/f 

9 

10. Hydrogen gas 

43,300 Btu/lb 

2. 51x1 0~ 3 kWh/f 
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11. HF (48X) 

7,000 Btu/lb 

5.22 kWh/f 

5 

12. Nitrogen (liquid) 

• 1,330 Btu/lb 

0.69 KWh/f. 

5 

13. Nitrogen (gas from 

- 

1.44xl0~ 3 kWh/i 

10 

liquid) 




14. Nitric Acid 

14,500 Btu/lb 

13.12 kWh H 
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(HN0 3 , 67X) 




15. Oxygen gas (0 2 ) 

830 Btu/lb 

7.66 kWh/* 

5 

16. Phosphorus (solid) 

23,790 Btu/lb 

1.54xl0" 2 kWh/g 

5 

17. Phosphine gas (PH^) 

- 

0.18 kWh/l 

11 

18. Phosphorous 


0.14 kWh/g 

12 

oxychloride (POCij) 

- 

19. Plating resist 

8,000 Btu/lb 

5.20 kWh/i 

13 

j 20. Silver 

1260xl0 6 Btu/ton 

0.406 kWh/g 

8 

21. Sodium Hydroxide 

by-product 

0 

5 

(NaOH) 




22. Toluene 

0.05 lb of low-press, 
steam plus 
0.025 kWh(ei)/lb 

0.0349 kWh/nf 

14 

23. Vacuua pump 
oil 

3,000 kWh(th)/ 
barrel 

18.87 kWh/* 
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15 

24. Wax 

2,000 kWh/barrel 

1.3xl0 -2 kWh/g 

16 
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References for Table V-l 


1. Assumed to be the same as that of acetone, which was 
taken from "Battelle Columbus Laboratories, Draft Target 
and Support Document for Developing a Maximum Energy 
Efficient Improvement Target for SIC 28',' Federal Energy 
Administration, Washington, DC (1976). 

2. J.T. Reding and B.P. Shepherd, "Energy Consumption, 

Report EPA-G5Q/2-75-032bV US Environmental Protection 
Agency, Washington, DC (4/75). 

3. Federal Energy Administration, "Project Independence 
Blueprint," Vol. 3, Federal Energy Administration, 

Washington, DC (1974). 

4. Calculated from the energy content for NH3 and using 
a density of 0.824 g/mil for the NH4GH solution. 

5. Battelle Columbus Laboratories, "Draft Target and Support 
Document for Developing a Maximum Energy Efficient 
Improvement Target for SIC 28',' Federal Energy Administra- 
tion, Washington, DC (1976). 

6. Assumed to be the same as that of butyl alcohol, as given in Ref. 5. 

7. H.W. Lownie, et al (Battelle Columbus Laboratories), 

"Draft Target Report on Development and Establishment of 
Energy Efficiency Improvement Targets for Primary Metal 
Industries," Federal Energy Administration, Washington, 

DC (9/76). 

8. M.G. Coleman, et al.. Motorola Final Report, DOE /JPL- 
954847-78/4, 183(11/78). 

9. Estimated as approximately 50% more than the average 
energy content of ccntnercial gases, as given in Ref. 5. 

10. Calculated from the energy content for LN2 assuming 
480 l of gas can be obtained from 18. of LN2 . 

11. Estimated as twice the energy content for P(s) plus 
3/2H 2 (g). 

12. Estimated as twice the energy content of P(s), plus ^02^)» 
plus |ci2(g)- The energy content of CI2 was 

taken from ref. 5. 

13. Assumed to be the same as that of wood rosin, as given in Ref. 5. 

1 

14. Taken as that of benzene, as given in ref. 5. 
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2. Energy Payback Time* 

To calculate payback times for photovoltaic systems, 
the annual energy output of a unit module area must be known. 

This power output, for photovoltaic systems, is 
dependent upon the encapsulated cell efficiency, the module 
packing factor, the solar insolation, the efficiency of 
the power conditioning and storage subsystems, and the mismatch 
of the energy availability and demand statistics. This list of 
dependencies makes it clear that the output of the module will vary 
according to climate at the locality of installation, and to the indi- 
vidual load to be satisfied. The factors of influence on the useful 
system output, outside of the cell efficiency, are combined into 
the "capacity factor". As a reasonably representative, not 
too optimistic value for this capacity factor, the number 0.11 

has been used, in consequence of the results of several system 
(4) 

studies. The capacity factor is essentially independent 
of the solar cell or module manufacturing process. Consequently, 
the energy payback time is only partly a function 
of the energy consumption for the solar module fabrication pro- 
cess, since the other system parameters strongly influence 
its absolute magnitude. In the calculations of energy pay- 
bark times, encapsulated cell efficiencies n of 11.5% and 
15% were employed for 1978 and 1982 and for 19Rf, , respectively, 

assuming the use of FFG ribbons, of 12%. To oLt.iin the 
energy payback times, the annual module output was calculated 

as E <■ n x 8766 x 0.11xHp k , H pk being the standardized peak 
solar irradiance, used as lkw(pk)/m 2 . The factor 0.11 repre- 
sents the "capacity factor". 
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Thus , annual energy outputs from a square meter of solar cells 
of 110.9, 144.6, and 115.7 kWh/m 2 , y are obtained for 1978, 

1982, and 1986 respectively. 

It is also to be observed that, since the energy recovery 

from solar cells is in the form of electrical energy, it is 

appropriate to convert those energy expenditures which occur 

in the form of heat of combustion, to equivalent electrical 

energy by applying the average efficiency of 0.30 experienced 

by the electric utilities in the conversion from heat of 

combustion of fuels to electric power delivered to the con- 
(5 ) 


sumer 



3. Energy Consumption in Photovoltaic Solar Array 
Manufacturing Process Sequences 
3.1 Data Sources. 

The major sources of information for this energy analysis 
were data accumulated from LSA project reports and industry . 
interviews. The energy consumed through materials, both 
direct and indirect materials, was obtained by converting 
the material consumption to energy units through multi- 
plication with the material energy contents shown in Table V-l. 

Where specific information to the contrary was not available, the material 
energy values were assumed to be in the form of thermal energy. The 
equipment costs, were converted to energy units expressed on the 
basis of unit cell area of throughput, by assuming a lifetime 
of seven years, and an energy content of the equipment which corres- 
ponds in value to 2% of the equipment cost.^ 6 ^ This energy cost 

has then been converted to a thermal energy using an energy 

17 ) 

price of $0 . 003/kWh (th) ' . Finally, the general energy 

usage for operating the facility was derived from the machine floor 
area by using the annual SAMICS utility cost of $3.74/ft^ and assu- 
ming that all of this "indirect" energy, since it is used primarily for 
lighting , air conditioning, and ventilating, in view of the 

high heat load in the building, is in the form of electricity 

( 8 ) 

at a cost of $0. 0319/kWh . The conversion factors for 

the equipment and facility energies thus become 0.9523 kwh(th)/$y 

and 1262 kWh/m 2 y, respectively. 
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Energy Contents of Selected Materials 






Material 

Original Units 
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1. Acetic acid 

3.28 lb o' hi-press. 
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13 
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21. Sodium Hydroxide 

by-product 
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0.05 lb of low-press, 
steam plus 
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11. Estimated as twice the energy content for P(s) plus 
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plus |ci2(g)- The energy content of CI2 was 

taken from ref. 5. 

13. Assumed to be the same as that of wood rosin, as given in Ref. 5. 
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15. Taken as approximately twice the energy content of mid 
die oil distillates. 

16. From M. Sittig, Practical Techniques for Saving Energy 
in the Chemical Petroleum and Metals Industries, 

Noyes Data Corporation, Park Ridge, NJ (1977), and 
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3.2 Energy Consumption in Si reduction and purification. 

This process group starts with the reduction of Si0 2 in an arc 
furnace . This is a quite efficient and cost-effective pro- 
cess which is not very likely to be replaced by another 
approach. In contrast, the following step of Si purification 
is very inefficient and likely to be replaced by one of 
several alternate methods under development. Furthest advanced 

among these is the SiH purification process being prepared 

4 

for pilot line operations by Union Carbide Corp. It can be 

expected to be a production process by 1986. 

a) Theoretical Material Balance for the Arc Furnace 
Process. 


The chemical reaction of this process is: 


Atomic Weights 


Masses, normal- 
ized to 1 kg Si 
out: 


, f 28t 

i { 2 . 


Si0 2 + 2C 


2(16) + 2(12) 
60 + 24 


Si + 2 CO 

28 + 2(12 + 16) 
28 + 56 


14 kg + 0.86 kg -*■ 1 kg + 2 kg 


According to Dow-Corning O) . r the industrially experienced 
and actual conversion efficiency of SiO^ to Si is 80%. 

The required input is thus: 2.68 kg Si0 2 for 1 kg Si out 

(Dow-Corning shows 2.71 kg SiOj) . 

It is assumed that half of the lost SiO^ input, or 0.2 7 kg, 
is used in the reaction: 

SiO 


Atomic Weights: 


2 

60 + 


C 

12 


Si o 

44 


+ 

+ 


CO 

28 


0.27kg + 0.05 kg * 0.2 kg + 0.12 kg 

where the SiO is lost at the top of the furnace. 

Thus, the total theoretical carbon input would be 0.91 kg C 
per 1 kg MG-Si output. 
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2. Energy Payback Time* 

To calculate payback times for photovoltaic systems, 
the annual energy output of a unit module area must be known. 

This power output, for photovoltaic systems, is 
dependent upon the encapsulated cell efficiency, the module 
packing factor, the solar insolation, the efficiency of 
the power conditioning and storage subsystems, and the mismatch 
of the energy availability and demand statistics. This list of 
dependencies makes it clear that the output of the module will vary 
according to climate at the locality of installation, and to the indi- 
vidual load to be satisfied. The factors of influence on the useful 
system output, outside of the cell efficiency, are combined into 
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Thus , annual energy outputs from a square meter of solar cells 
of 110.9, 144.6, and 115.7 kWh/m 2 , y are obtained for 1978, 

1982, and 1986 respectively. 

It is also to be observed that, since the energy recovery 

from solar cells is in the form of electrical energy, it is 

appropriate to convert those energy expenditures which occur 

in the form of heat of combustion, to equivalent electrical 

energy by applying the average efficiency of 0.30 experienced 

by the electric utilities in the conversion from heat of 

combustion of fuels to electric power delivered to the con- 
(5 ) 
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3. Energy Consumption in Photovoltaic Solar Array 
Manufacturing Process Sequences 
3.1 Data Sources. 

The major sources of information for this energy analysis 
were data accumulated from LSA project reports and industry . 
interviews. The energy consumed through materials, both 
direct and indirect materials, was obtained by converting 
the material consumption to energy units through multi- 
plication with the material energy contents shown in Table V-l. 

Where specific information to the contrary was not available, the material 
energy values were assumed to be in the form of thermal energy. The 
equipment costs, were converted to energy units expressed on the 
basis of unit cell area of throughput, by assuming a lifetime 
of seven years, and an energy content of the equipment which corres- 
ponds in value to 2% of the equipment cost.^ 6 ^ This energy cost 

has then been converted to a thermal energy using an energy 

17 ) 

price of $0 . 003/kWh (th) ' . Finally, the general energy 

usage for operating the facility was derived from the machine floor 
area by using the annual SAMICS utility cost of $3.74/ft^ and assu- 
ming that all of this "indirect" energy, since it is used primarily for 
lighting , air conditioning, and ventilating, in view of the 

high heat load in the building, is in the form of electricity 

( 8 ) 

at a cost of $0. 0319/kWh . The conversion factors for 

the equipment and facility energies thus become 0.9523 kwh(th)/$y 

and 1262 kWh/m 2 y, respectively. 
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Table V-l 

Energy Contents of Selected Materials 






Material 

Original Units 

Convenient Units 
(Thermal kWh's! 

Reference 

1. Acetic acid 

3.28 lb o' hi-press. 
ateaa jIus 
0.07 Wh(*v)/lb 

3.18 kWh/g 

1 

2. Aluminum (At) 

520x10® kWh(th)/ton 

0.17 kWh/g 

2 

3. Ammonia gaa (NH^) 

8.0Sxl0 3 Btu/lb 

0.00534 kWh/g 

3 

A. A— onlum Hydroxide 

• 

1.32 kWh/i 

4 

30% (hh 4 oh> 




S. Argon gas 

1,100 Btu/lb 

1.27xl0 -3 kWh/l 

5 

6. Butyl acetate 

4.32 lb of low-press, 
steam plus 
0.082 kWh(ei)/lb 

3.21 kWh/g 

6 

7. Copper (Cu) 

16.2x10* Btu/ton 

5.23 kVh/g 

7 

8. Energy for exhausting 
waste fuses without 

0.46 kW/1000 ft 3 /sln 

2.56x10" 2 kWh/ 
1000 ft 3 

8 

scrubbing. 




9. Freon-14 gas (CF ft ) 


2xlO -3 kWh/f 

9 

10. Hydrogen gas 

43,300 Btu/lb 

2. 51x1 0~ 3 kWh/* 

5 

11. HF (48X) 

7,000 Btu/lb 

5.22 kWh/f 

5 

12. Nitrogen (liquid) 

• 1,330 Btu/lb 

0.69 KWh/i 

5 

13. Nitrogen (gas from 

- 

1.44xl0~ 3 kWh/* 

10 

liquid) 




14. Nitric Acid 

14,500 Btu/lb 

13.12 kWh /t 

5 

(HN0 3 , 67X) 




15. Oxygen gas (0 2 ) 

830 Btu/lb 

7.66 kWh/* 

5 

16. Phosphorus (solid) 

23,790 Btu/lb 

1.54xl0" 2 kWh/g 

5 

17. Phosphine gas (PH^) 

- 

0.18 kWh/l 

11 

18. Phosphorous 


0.14 kWh/g 

12 

oxychloride (POCij) 

- 

19. Plating resist 

8,000 Btu/lb 

5.20 kWh/i 

13 

j 20. Silver 

1260xl0 6 Btu/ton 

0.406 kWh/g 

8 

21. Sodium Hydroxide 

by-product 

0 

5 

(NaOH) 




22. Toluene 

0.05 lb of low-press, 
steam plus 
0.025 kWh(ei)/lb 

0.0349 kWh/m* 

14 

23. Vacuua pump 
oil 

3,000 kWh(th)/ 
barrel 

18.87 kUh/l 

i 

15 

24. Wax 

2,000 kWh/barrel 

1.3xl0 -2 kWh/g 

16 
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References for Table V-l 


1. Assumed to be the same as that of acetone, which was 
taken from "Battelle Columbus Laboratories, Draft Target 
and Support Document for Developing a Maximum Energy 
Efficient Improvement Target for SIC 28',' Federal Energy 
Administration, Washington, DC (1976). 

2. J.T. Reding and B.P. Shepherd, "Energy Consumption, 

Report EPA-G5Q/2-75-032bV US Environmental Protection 
Agency, Washington, DC (4/75). 

3. Federal Energy Administration, "Project Independence 
Blueprint," Vol. 3, Federal Energy Administration, 

Washington, DC (1974). 

4. Calculated from the energy content for NH3 and using 
a density of 0.824 q/ml for the NH4GH solution. 

5. Battelle Columbus Laboratories, "Draft Target and Support 
Document for Developing a Maximum Energy Efficient 
Improvement Target for SIC 28',' Federal Energy Administra- 
tion, Washington, DC (1976). 

6. Assumed to be the same as that of butyl alcohol, as given in Ref. 5. 

7. H.W. Lownie, et al (Battelle Columbus Laboratories), 

"Draft Target Report on Development and Establishment of 
Energy Efficiency Improvement Targets for Primary Metal 
Industries," Federal Energy Administration, Washington, 

DC (9/76). 

8. M.G. Coleman, et al.. Motorola Final Report, DOE /JPL- 
954847-78/4, 183(11/78). 

9. Estimated as approximately 50% more than the average 
energy content of ccntnercial gases, as given in Ref. 5. 

10. Calculated from the energy content for LN2 assuming 
480 l of gas can be obtained from 18. of LN2 . 

11. Estimated as twice the energy content for P(s) plus 
3/2H 2 (g). 

12. Estimated as twice the energy content of P(s), plus ^02^)» 
plus |ci2(g)- The energy content of CI2 was 

taken from ref. 5. 

13. Assumed to be the same as that of wood rosin, as given in Ref. 5. 

1 

14. Taken as that of benzene, as given in ref. 5. 
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15. Taken as approximately twice the energy content of mid 
die oil distillates. 

16. From M. Sittig, Practical Techniques for Saving Energy 
in the Chemical Petroleum and Metals Industries, 

Noyes Data Corporation, Park Ridge, NJ (1977), and 
using a density of 0.97 g/ml for the wax. 
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3.2 Energy Consumption in Si reduction and purification. 

This process group starts with the reduction of Si0 2 in an arc 
furnace . This is a quite efficient and cost-effective pro- 
cess which is not very likely to be replaced by another 
approach. In contrast, the following step of Si purification 
is very inefficient and likely to be replaced by one of 
several alternate methods under development. Furthest advanced 

among these is the SiH purification process being prepared 

4 

for pilot line operations by Union Carbide Corp. It can be 

expected to be a production process by 1986. 

a) Theoretical Material Balance for the Arc Furnace 
Process. 


The chemical reaction of this process is: 


Atomic Weights 


Masses, normal- 
ized to 1 kg Si 
out: 


, f 28t 

i { 2 . 


Si0 2 + 2C 


2(16) + 2(12) 
60 + 24 


Si + 2 CO 

28 + 2(12 + 16) 
28 + 56 


14 kg + 0.86 kg -*■ 1 kg + 2 kg 


According to Dow-Corning O) . r the industrially experienced 
and actual conversion efficiency of SiO^ to Si is 80%. 

The required input is thus: 2.68 kg Si0 2 for 1 kg Si out 

(Dow-Corning shows 2.71 kg SiOj) . 

It is assumed that half of the lost SiO^ input, or 0.2 7 kg, 
is used in the reaction: 

SiO 


Atomic Weights: 


2 

60 + 


C 

12 


Si o 

44 


+ 

+ 


CO 

28 


0.27kg + 0.05 kg * 0.2 kg + 0.12 kg 

where the SiO is lost at the top of the furnace. 

Thus, the total theoretical carbon input would be 0.91 kg C 
per 1 kg MG-Si output. 
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2. Energy Payback Time* 

To calculate payback times for photovoltaic systems, 
the annual energy output of a unit module area must be known. 

This power output, for photovoltaic systems, is 
dependent upon the encapsulated cell efficiency, the module 
packing factor, the solar insolation, the efficiency of 
the power conditioning and storage subsystems, and the mismatch 
of the energy availability and demand statistics. This list of 
dependencies makes it clear that the output of the module will vary 
according to climate at the locality of installation, and to the indi- 
vidual load to be satisfied. The factors of influence on the useful 
system output, outside of the cell efficiency, are combined into 
the "capacity factor". As a reasonably representative, not 
too optimistic value for this capacity factor, the number 0.11 

has been used, in consequence of the results of several system 
(4) 

studies. The capacity factor is essentially independent 
of the solar cell or module manufacturing process. Consequently, 
the energy payback time is only partly a function 
of the energy consumption for the solar module fabrication pro- 
cess, since the other system parameters strongly influence 
its absolute magnitude. In the calculations of energy pay- 
bark times, encapsulated cell efficiencies n of 11.5% and 
15% were employed for 1978 and 1982 and for 1980, respectively, 

assuming the use nl' FFG ribbons, of 12%. To obtain the 
energy payback times, the annual module output was calculated 

as E <■ n x 8766 x O.llxHp^, Hp k being the standardized peak 
solar irradiance, used as lkw(pk)/m 2 . The factor 0.11 repre- 
sents the "capacity factor". 
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Thus, annual energy outputs from a square meter of solar cells 
of 110.9, 144.6, and 115.7 kWh/m 2 , y are obtained for 1978, 

1982, and 1986 respectively. 

It is also to be observed that, since the energy recovery 

from solar cells is in the form of electrical energy, it is 

appropriate to convert those energy expenditures which occur 

in the form of heat of combustion, to equivalent electrical 

energy by applying the average efficiency of 0.30 experienced 

by the electric utilities in the conversion from heat of 

combustion of fuels to electric power delivered to the con- 
(5 ) 


sumer 


3. Energy Consumption in Photovoltaic Solar Array 
Manufacturing Process Sequences 
3.1 Data Sources. 

The major sources of information for this energy analysis 
were data accumulated from LSA project reports and industry . 
interviews. The energy consumed through materials, both 
direct and indirect materials, was obtained by converting 
the material consumption to energy units through multi- 
plication with the material energy contents shown in Table V-l. 

Where specific information to the contrary was not available, the material 
energy values were assumed to be in the form of thermal energy. The 
equipment costs, were converted to energy units expressed on the 
basis of unit cell area of throughput, by assuming a lifetime 
of seven years, and an energy content of the equipment which corres- 
ponds in value to 2% of the equipment cost.^ 6 ^ This energy cost 

has then been converted to a thermal energy using an energy 

17 ) 

price of $0 . 003/kWh (th) ' . Finally, the general energy 

usage for operating the facility was derived from the machine floor 
area by using the annual SAMICS utility cost of $3.74/ft^ and assu- 
ming that all of this "indirect" energy, since it is used primarily for 
lighting , air conditioning, and ventilating, in view of the 

high heat load in the building, is in the form of electricity 

( 8 ) 

at a cost of $0. 0319/kWh . The conversion factors for 

the equipment and facility energies thus become 0.9523 kwh(th)/$y 

and 1262 kWh/m 2 y, respectively. 
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of the energy availability and demand statistics. This list of 
dependencies makes it clear that the output of the module will vary 
according to climate at the locality of installation, and to the indi- 
vidual load to be satisfied. The factors of influence on the useful 
system output, outside of the cell efficiency, are combined into 
the "capacity factor". As a reasonably representative, not 
too optimistic value for this capacity factor, the number 0.11 

has been used, in consequence of the results of several system 
(4) 

studies. The capacity factor is essentially independent 
of the solar cell or module manufacturing process. Consequently, 
the energy payback time is only partly a function 
of the energy consumption for the solar module fabrication pro- 
cess, since the other system parameters strongly influence 
its absolute magnitude. In the calculations of energy pay- 
bark times, encapsulated cell efficiencies n of 11.5% and 
15% were employed for 1978 and 1982 and for 1980, respectively, 

assuming the use nl' FFG ribbons, of 12%. To obtain the 
energy payback times, the annual module output was calculated 

as E <■ n x 8766 x O.llxHp^, Hp k being the standardized peak 
solar irradiance, used as lkw(pk)/m 2 . The factor 0.11 repre- 
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Table V-l 

Energy Contents of Selected Materials 






Material 

Original Units 

Convenient Units 
(Thermal kWh 'si 

Reference 

1. Acetic ecld 

3.28 lb o' hi-press. 
steam ^lus 
0.07 Wh(*v)/lb 

3.18 kWh/g 

1 

2. Aluminum (At) 

520x10® kWh(th)/ton 

0.17 kWh/g 

2 

3. Ammonia gaa (NH^) 

8.0Sxl0 3 Btu/lb 

0.00534 kWh/g 

3 

A. Ammonium Hydroxide 

• 

1.32 kWh/1 

4 

30% (nh 4 oh> 




S. Argon gas 

1,100 Btu/lb 

1.27xl0 -3 kWh/l 

5 

6. Butyl acetate 

4.32 lb of low-press, 
steam plus 
0.082 kWh(el)/lb 

3.21 kWh/g 

6 

7. Copper (Cu) 

16.2x10* Btu/ton 

5.23 kWh/g 

7 

8. Energy for exhausting 
waste fuses without 

0.46 kW/1000 ft 3 /min 

2.56x10" 2 kWh/ 
1000 ft 3 

8 

scrubbing. 




9. Freon-14 gas (CF ft ) 


2xlO -3 kWh/f 

9 

10. Hydrogen gas 

43,300 Btu/lb 

2. 51x1 0~ 3 kWh/l 

5 

11. HF (48X) 

7,000 Btu/lb 

5.22 kWh/f 

5 

12. Nitrogen (liquid) 

• 1,330 Btu/lb 

0.69 KWh/f. 

5 

13. Nitrogen (gas from 

- 

1.44xl0~ 3 kWh/i 

10 

liquid) 




14. Nitric Acid 

14,500 Btu/lb 

13.12 kWh It 

5 

(HN0 3 , 67X) 




15. Oxygen gas (0 2 ) 

830 Btu/lb 

7.66 kWh/* 

5 

16. Phosphorus (solid) 

23,790 Btu/lb 

1.54xl0" 2 kWh/g 

5 

17. Phosphine gas (PH^) 

- 

0.18 kWh/t 

11 

18. Phosphorous 


0.14 kWh/g 

12 

oxychloride (POCtj) 

- 

19. Plating resist 

8,000 Btu/lb 

5.20 kMh/1 

13 

j 20. Silver 

1260xl0 6 Btu/ton 

0.406 kWh/g 

8 

21. Sodium Hydroxide 

by-product 

0 

5 

(NaOH) 




22. Toluene 

0.05 lb of low-press, 
steam plus 
0.025 kWh(ei)/lb 

0.0349 kWh/nf 

14 

23. Vacuua pump 
oil 

3,000 kWh(th)/ 
barrel 

18.87 kWh/t 

i 

15 

24. Wax 

2,000 kWh/barrel 

1.3xl0 -2 kWh/g 

16 
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Agency, Washington, DC (4/75). 

3. Federal Energy Administration, "Project Independence 
Blueprint," Vol. 3, Federal Energy Administration, 

Washington, DC (1974). 

4. Calculated from the energy content for NH3 and using 
a density of 0.824 g/mil for the NH4GH solution. 

5. Battelle Columbus Laboratories, "Draft Target and Support 
Document for Developing a Maximum Energy Efficient 
Improvement Target for SIC 28" Federal Energy Administra- 
tion, Washington, DC (1976). 

6. Assumed to be the same as that of butyl alcohol, as given in Ref. 5. 

7. H.W. Lownie, et al (Battelle Columbus Laboratories), 

"Draft Target Report on Development and Establishment of 
Energy Efficiency Improvement Targets for Primary Metal 
Industries," Federal Energy Administration, Washington, 

DC (9/76). 

8. M.G. Coleman, et al.. Motorola Final Report, DOE /JPL- 
954847-78/4, 183(11/78). 

9. Estimated as approximately 50% more than the average 
energy content of ccntnercial gases, as given in Ref. 5. 

10. Calculated from the energy content for LN2 assuming 
480 l of gas can be obtained from 18. of LN2 . 

11. Estimated as twice the energy content for P(s) plus 
3/2H 2 (g). 

12. Estimated as twice the energy content of P(s), plus ^02^)» 
plus |ci2(9)- The energy content of CI2 was 

taken from ref. 5. 

13. Assumed to be the same as that of wood rosin, as given in Ref. 5. 

\ 

14. Taken as that of benzene, as given in ref. 5. 
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15. Taken as approximately twice the energy content of mid 
die oil distillates. 

16. From M. Sittig, Practical Techniques for Saving Energy 
in the Chemical Petroleum and Metals Industries, 

Noyes Data Corporation, Park Ridge, NJ (1977), and 
using a density of 0.97 g/ml for the wax. 
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3.2 Energy Consumption in Si reduction and purification. 

This process group starts with the reduction of Si0 2 in an arc 
furnace . This is a quite efficient and cost-effective pro- 
cess which is not very likely to be replaced by another 
approach. In contrast, the following step of Si purification 
is very inefficient and likely to be replaced by one of 
several alternate methods under development. Furthest advanced 

among these is the SiH purification process being prepared 

4 

for pilot line operations by Union Carbide Corp. It can be 

expected to be a production process by 1986. 

a) Theoretical Material Balance for the Arc Furnace 
Process. 


The chemical reaction of this process is: 


Atomic Weights 


Masses, normal- 
ized to 1 kg Si 
out: 


, f 28t 

i { 2 . 


Si0 2 + 2C 


2(16) + 2(12) 
60 + 24 


Si + 2 CO 

28 + 2(12 + 16) 
28 + 56 


14 kg + 0.86 kg -*■ 1 kg + 2 kg 


According to Dow-Corning O) . r the industrially experienced 
and actual conversion efficiency of SiO^ to Si is 80%. 

The required input is thus: 2.68 kg Si0 2 for 1 kg Si out 

(Dow-Corning shows 2.71 kg SiOj) . 

It is assumed that half of the lost SiO^ input, or 0.2 7 kg, 
is used in the reaction: 

SiO 


Atomic Weights: 


2 

60 + 


C 

12 


Si o 

44 


+ 

+ 


CO 

28 


0.27kg + 0.05 kg * 0.2 kg + 0.12 kg 

where the SiO is lost at the top of the furnace. 

Thus, the total theoretical carbon input would be 0.91 kg C 
per 1 kg MG-Si output. 
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2. Energy Payback Time* 

To calculate payback times for photovoltaic systems, 
the annual energy output of a unit module area must be known. 

This power output, for photovoltaic systems, is 
dependent upon the encapsulated cell efficiency, the module 
packing factor, the solar insolation, the efficiency of 
the power conditioning and storage subsystems, and the mismatch 
of the energy availability and demand statistics. This list of 
dependencies makes it clear that the output of the module will vary 
according to climate at the locality of installation, and to the indi- 
vidual load to be satisfied. The factors of influence on the useful 
system output, outside of the cell efficiency, are combined into 
the "capacity factor". As a reasonably representative, not 
too optimistic value for this capacity factor, the number 0.11 

has been used, in consequence of the results of several system 
(4) 

studies. The capacity factor is essentially independent 
of the solar cell or module manufacturing process. Consequently, 
the energy payback time is only partly a function 
of the energy consumption for the solar module fabrication pro- 
cess, since the other system parameters strongly influence 
its absolute magnitude. In the calculations of energy pay- 
bark times, encapsulated cell efficiencies n of 11.5% and 
15% were employed for 1978 and 1982 and for 1980, respectively, 

assuming the use nl' FFG ribbons, of 12%. To obtain the 
energy payback times, the annual module output was calculated 

as E <■ n x 8766 x O.llxHp^, Hp k being the standardized peak 
solar irradiance, used as lkw(pk)/m 2 . The factor 0.11 repre- 
sents the "capacity factor". 

V-2 


Thus, annual energy outputs from a square meter of solar cells 
of 110.9, 144.6, and 115.7 kWh/m 2 , y are obtained for 1978, 

1982, and 1986 respectively. 

It is also to be observed that, since the energy recovery 

from solar cells is in the form of electrical energy, it is 

appropriate to convert those energy expenditures which occur 

in the form of heat of combustion, to equivalent electrical 

energy by applying the average efficiency of 0.30 experienced 

by the electric utilities in the conversion from heat of 

combustion of fuels to electric power delivered to the con- 
(5 ) 


sumer 


3. Energy Consumption in Photovoltaic Solar Array 
Manufacturing Process Sequences 
3.1 Data Sources. 

The major sources of information for this energy analysis 
were data accumulated from LSA project reports and industry . 
interviews. The energy consumed through materials, both 
direct and indirect materials, was obtained by converting 
the material consumption to energy units through multi- 
plication with the material energy contents shown in Table V-l. 

Where specific information to the contrary was not available, the material 
energy values were assumed to be in the form of thermal energy. The 
equipment costs, were converted to energy units expressed on the 
basis of unit cell area of throughput, by assuming a lifetime 
of seven years, and an energy content of the equipment which corres- 
ponds in value to 2% of the equipment cost.^ 6 ^ This energy cost 

has then been converted to a thermal energy using an energy 

17 ) 

price of $0 . 003/kWh (th) ' . Finally, the general energy 

usage for operating the facility was derived from the machine floor 
area by using the annual SAMICS utility cost of $3.74/ft^ and assu- 
ming that all of this "indirect" energy, since it is used primarily for 
lighting , air conditioning, and ventilating, in view of the 

high heat load in the building, is in the form of electricity 

( 8 ) 

at a cost of $0. 0319/kWh . The conversion factors for 

the equipment and facility energies thus become 0.9523 kwh(th)/$y 

and 1262 kWh/m 2 y, respectively. 
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Table V-l 

Energy Contents of Selected Materials 






Material 

Original Units 

Convenient Units 
(Thermal kWh 'si 

Reference 

1. Acetic ecld 

3.28 lb o' hi-press. 
steam ^lus 
0.07 Wh(*v)/lb 

3.18 kWh/g 

1 

2. Aluminum (At) 

520x10® kWh(th)/ton 

0.17 kWh/g 

2 

3. Ammonia gaa (NH^) 

8.0Sxl0 3 Btu/lb 

0.00534 kWh/g 

3 

A. Ammonium Hydroxide 

• 

1.32 kWh/1 

4 

30% (nh 4 oh> 




S. Argon gas 

1,100 Btu/lb 

1.27xl0 -3 kWh/l 

5 

6. Butyl acetate 

4.32 lb of low-press, 
steam plus 
0.082 kWh(el)/lb 

3.21 kWh/g 

6 

7. Copper (Cu) 

16.2x10* Btu/ton 

5.23 kWh/g 

7 

8. Energy for exhausting 
waste fuses without 

0.46 kW/1000 ft 3 /min 

2.56x10" 2 kWh/ 
1000 ft 3 

8 

scrubbing. 




9. Freon-14 gas (CF ft ) 


2xlO -3 kWh/f 

9 

10. Hydrogen gas 

43,300 Btu/lb 

2. 51x1 0~ 3 kWh/l 

5 

11. HF (48X) 

7,000 Btu/lb 

5.22 kWh/f 

5 

12. Nitrogen (liquid) 

• 1,330 Btu/lb 

0.69 KWh/f. 

5 

13. Nitrogen (gas from 

- 

1.44xl0~ 3 kWh/i 

10 

liquid) 




14. Nitric Acid 

14,500 Btu/lb 

13.12 kWh It 

5 

(HN0 3 , 67X) 




15. Oxygen gas (0 2 ) 

830 Btu/lb 

7.66 kWh/* 

5 

16. Phosphorus (solid) 

23,790 Btu/lb 

1.54xl0" 2 kWh/g 

5 

17. Phosphine gas (PH^) 

- 

0.18 kWh/t 

11 

18. Phosphorous 


0.14 kWh/g 

12 

oxychloride (POCtj) 

- 

19. Plating resist 

8,000 Btu/lb 

5.20 kMh/1 

13 

j 20. Silver 

1260xl0 6 Btu/ton 

0.406 kWh/g 

8 

21. Sodium Hydroxide 

by-product 

0 

5 

(NaOH) 




22. Toluene 

0.05 lb of low-press, 
steam plus 
0.025 kWh(ei)/lb 

0.0349 kWh/nf 

14 

23. Vacuua pump 
oil 

3,000 kWh(th)/ 
barrel 

18.87 kWh/t 

i 

15 

24. Wax 

2,000 kWh/barrel 

1.3xl0 -2 kWh/g 

16 


2. Energy Payback Time* 

To calculate payback times for photovoltaic systems, 
the annual energy output of a unit module area must be known. 

This power output, for photovoltaic systems, is 
dependent upon the encapsulated cell efficiency, the module 
packing factor, the solar insolation, the efficiency of 
the power conditioning and storage subsystems, and the mismatch 
of the energy availability and demand statistics. This list of 
dependencies makes it clear that the output of the module will vary 
according to climate at the locality of installation, and to the indi- 
vidual load to be satisfied. The factors of influence on the useful 
system output, outside of the cell efficiency, are combined into 
the "capacity factor". As a reasonably representative, not 
too optimistic value for this capacity factor, the number 0.11 

has been used, in consequence of the results of several system 
(4) 

studies. The capacity factor is essentially independent 
of the solar cell or module manufacturing process. Consequently, 
the energy payback time is only partly a function 
of the energy consumption for the solar module fabrication pro- 
cess, since the other system parameters strongly influence 
its absolute magnitude. In the calculations of energy pay- 
bark times, encapsulated cell efficiencies n of 11.5% and 
15% were employed for 1978 and 1982 and for 1980, respectively, 

assuming the use nl' FFG ribbons, of 12%. To obtain the 
energy payback times, the annual module output was calculated 

as E <■ n x 8766 x O.llxHp^, Hp k being the standardized peak 
solar irradiance, used as lkw(pk)/m 2 . The factor 0.11 repre- 
sents the "capacity factor". 
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Thus, annual energy outputs from a square meter of solar cells 
of 110.9, 144.6, and 115.7 kWh/m 2 , y are obtained for 1978, 

1982, and 1986 respectively. 

It is also to be observed that, since the energy recovery 

from solar cells is in the form of electrical energy, it is 

appropriate to convert those energy expenditures which occur 

in the form of heat of combustion, to equivalent electrical 

energy by applying the average efficiency of 0.30 experienced 

by the electric utilities in the conversion from heat of 

combustion of fuels to electric power delivered to the con- 
(5 ) 
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3. Energy Consumption in Photovoltaic Solar Array 
Manufacturing Process Sequences 
3.1 Data Sources. 

The major sources of information for this energy analysis 
were data accumulated from LSA project reports and industry . 
interviews. The energy consumed through materials, both 
direct and indirect materials, was obtained by converting 
the material consumption to energy units through multi- 
plication with the material energy contents shown in Table V-l. 

Where specific information to the contrary was not available, the material 
energy values were assumed to be in the form of thermal energy. The 
equipment costs, were converted to energy units expressed on the 
basis of unit cell area of throughput, by assuming a lifetime 
of seven years, and an energy content of the equipment which corres- 
ponds in value to 2% of the equipment cost.^ 6 ^ This energy cost 

has then been converted to a thermal energy using an energy 

17 ) 

price of $0 . 003/kWh (th) ' . Finally, the general energy 

usage for operating the facility was derived from the machine floor 
area by using the annual SAMICS utility cost of $3.74/ft^ and assu- 
ming that all of this "indirect" energy, since it is used primarily for 
lighting , air conditioning, and ventilating, in view of the 

high heat load in the building, is in the form of electricity 
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at a cost of $0. 0319/kWh . The conversion factors for 

the equipment and facility energies thus become 0.9523 kwh(th)/$y 

and 1262 kWh/m 2 y, respectively. 
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References for Table V-l 


1. Assumed to be the same as that of acetone, which was 
taken from "Battelle Columbus Laboratories, Draft Target 
and Support Document for Developing a Maximum Energy 
Efficient Improvement Target for SIC 28',' Federal Energy 
Administration, Washington, DC (1976). 

2. J.T. Reding and B.P. Shepherd, "Energy Consumption, 

Report EPA-G5Q/2-75-032bV US Environmental Protection 
Agency, Washington, DC (4/75). 

3. Federal Energy Administration, "Project Independence 
Blueprint," Vol. 3, Federal Energy Administration, 

Washington, DC (1974). 

4. Calculated from the energy content for NH3 and using 
a density of 0.824 g/mil for the NH4GH solution. 

5. Battelle Columbus Laboratories, "Draft Target and Support 
Document for Developing a Maximum Energy Efficient 
Improvement Target for SIC 28" Federal Energy Administra- 
tion, Washington, DC (1976). 

6. Assumed to be the same as that of butyl alcohol, as given in Ref. 5. 

7. H.W. Lownie, et al (Battelle Columbus Laboratories), 

"Draft Target Report on Development and Establishment of 
Energy Efficiency Improvement Targets for Primary Metal 
Industries," Federal Energy Administration, Washington, 

DC (9/76). 

8. M.G. Coleman, et al.. Motorola Final Report, DOE /JPL- 
954847-78/4, 183(11/78). 

9. Estimated as approximately 50% more than the average 
energy content of ccntnercial gases, as given in Ref. 5. 

10. Calculated from the energy content for LN2 assuming 
480 l of gas can be obtained from 18. of LN2 . 

11. Estimated as twice the energy content for P(s) plus 
3/2H 2 (g). 

12. Estimated as twice the energy content of P(s), plus ^02^)» 
plus |ci2(9)- The energy content of CI2 was 

taken from ref. 5. 

13. Assumed to be the same as that of wood rosin, as given in Ref. 5. 

\ 

14. Taken as that of benzene, as given in ref. 5. 
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15. Taken as approximately twice the energy content of mid 
die oil distillates. 

16. From M. Sittig, Practical Techniques for Saving Energy 
in the Chemical Petroleum and Metals Industries, 

Noyes Data Corporation, Park Ridge, NJ (1977), and 
using a density of 0.97 g/ml for the wax. 
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3.2 Energy Consumption in Si reduction and purification. 

This process group starts with the reduction of Si0 2 in an arc 
furnace . This is a quite efficient and cost-effective pro- 
cess which is not very likely to be replaced by another 
approach. In contrast, the following step of Si purification 
is very inefficient and likely to be replaced by one of 
several alternate methods under development. Furthest advanced 

among these is the SiH purification process being prepared 

4 

for pilot line operations by Union Carbide Corp. It can be 

expected to be a production process by 1986. 

a) Theoretical Material Balance for the Arc Furnace 
Process. 


The chemical reaction of this process is: 


Atomic Weights 


Masses, normal- 
ized to 1 kg Si 
out: 


, f 28t 

i { 2 . 


Si0 2 + 2C 


2(16) + 2(12) 
60 + 24 


Si + 2 CO 

28 + 2(12 + 16) 
28 + 56 


14 kg + 0.86 kg -*■ 1 kg + 2 kg 


According to Dow-Corning O) . r the industrially experienced 
and actual conversion efficiency of SiO^ to Si is 80%. 

The required input is thus: 2.68 kg Si0 2 for 1 kg Si out 

(Dow-Corning shows 2.71 kg SiOj) . 

It is assumed that half of the lost SiO^ input, or 0.2 7 kg, 
is used in the reaction: 

SiO 


Atomic Weights: 


2 

60 + 


C 

12 


Si o 

44 


+ 

+ 


CO 

28 


0.27kg + 0.05 kg * 0.2 kg + 0.12 kg 

where the SiO is lost at the top of the furnace. 

Thus, the total theoretical carbon input would be 0.91 kg C 
per 1 kg MG-Si output. 
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2. Energy Payback Time* 

To calculate payback times for photovoltaic systems, 
the annual energy output of a unit module area must be known. 

This power output, for photovoltaic systems, is 
dependent upon the encapsulated cell efficiency, the module 
packing factor, the solar insolation, the efficiency of 
the power conditioning and storage subsystems, and the mismatch 
of the energy availability and demand statistics. This list of 
dependencies makes it clear that the output of the module will vary 
according to climate at the locality of installation, and to the indi- 
vidual load to be satisfied. The factors of influence on the useful 
system output, outside of the cell efficiency, are combined into 
the "capacity factor". As a reasonably representative, not 
too optimistic value for this capacity factor, the number 0.11 

has been used, in consequence of the results of several system 
(4) 

studies. The capacity factor is essentially independent 
of the solar cell or module manufacturing process. Consequently, 
the energy payback time is only partly a function 
of the energy consumption for the solar module fabrication pro- 
cess, since the other system parameters strongly influence 
its absolute magnitude. In the calculations of energy pay- 
bark times, encapsulated cell efficiencies n of 11.5% and 
15% were employed for 1978 and 1982 and for 1980, respectively, 

assuming the use nl' FFG ribbons, of 12%. To obtain the 
energy payback times, the annual module output was calculated 

as E <■ n x 8766 x O.llxHp^, Hp k being the standardized peak 
solar irradiance, used as lkw(pk)/m 2 . The factor 0.11 repre- 
sents the "capacity factor". 
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Thus, annual energy outputs from a square meter of solar cells 
of 110.9, 144.6, and 115.7 kWh/m 2 , y are obtained for 1978, 

1982, and 1986 respectively. 

It is also to be observed that, since the energy recovery 

from solar cells is in the form of electrical energy, it is 

appropriate to convert those energy expenditures which occur 

in the form of heat of combustion, to equivalent electrical 

energy by applying the average efficiency of 0.30 experienced 

by the electric utilities in the conversion from heat of 

combustion of fuels to electric power delivered to the con- 
(5 ) 
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3. Energy Consumption in Photovoltaic Solar Array 
Manufacturing Process Sequences 
3.1 Data Sources. 

The major sources of information for this energy analysis 
were data accumulated from LSA project reports and industry . 
interviews. The energy consumed through materials, both 
direct and indirect materials, was obtained by converting 
the material consumption to energy units through multi- 
plication with the material energy contents shown in Table V-l. 

Where specific information to the contrary was not available, the material 
energy values were assumed to be in the form of thermal energy. The 
equipment costs, were converted to energy units expressed on the 
basis of unit cell area of throughput, by assuming a lifetime 
of seven years, and an energy content of the equipment which corres- 
ponds in value to 2% of the equipment cost.^ 6 ^ This energy cost 

has then been converted to a thermal energy using an energy 

17 ) 

price of $0 . 003/kWh (th) ' . Finally, the general energy 

usage for operating the facility was derived from the machine floor 
area by using the annual SAMICS utility cost of $3.74/ft^ and assu- 
ming that all of this "indirect" energy, since it is used primarily for 
lighting , air conditioning, and ventilating, in view of the 

high heat load in the building, is in the form of electricity 

( 8 ) 

at a cost of $0. 0319/kWh . The conversion factors for 

the equipment and facility energies thus become 0.9523 kwh(th)/$y 

and 1262 kWh/m 2 y, respectively. 
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Industries," Federal Energy Administration, Washington, 

DC (9/76). 

8. M.G. Coleman, et al.. Motorola Final Report, DOE /JPL- 
954847-78/4, 183(11/78). 

9. Estimated as approximately 50% more than the average 
energy content of ccntnercial gases, as given in Ref. 5. 

10. Calculated from the energy content for LN2 assuming 
480 l of gas can be obtained from 18. of LN2 . 

11. Estimated as twice the energy content for P(s) plus 
3/2H 2 (g). 

12. Estimated as twice the energy content of P(s), plus ^02^)» 
plus |ci2(9)- The energy content of CI2 was 

taken from ref. 5. 

13. Assumed to be the same as that of wood rosin, as given in Ref. 5. 

\ 

14. Taken as that of benzene, as given in ref. 5. 
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References for Table V-l 


1. Assumed to be the same as that of acetone, which was 
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Agency, Washington, DC (4/75). 
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Washington, DC (1974). 
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a density of 0.824 g/mil for the NH4GH solution. 
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\ 
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15. Taken as approximately twice the energy content of mid 
die oil distillates. 

16. From M. Sittig, Practical Techniques for Saving Energy 
in the Chemical Petroleum and Metals Industries, 

Noyes Data Corporation, Park Ridge, NJ (1977), and 
using a density of 0.97 g/ml for the wax. 
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3.2 Energy Consumption in Si reduction and purification. 

This process group starts with the reduction of Si0 2 in an arc 
furnace . This is a quite efficient and cost-effective pro- 
cess which is not very likely to be replaced by another 
approach. In contrast, the following step of Si purification 
is very inefficient and likely to be replaced by one of 
several alternate methods under development. Furthest advanced 

among these is the SiH purification process being prepared 

4 

for pilot line operations by Union Carbide Corp. It can be 

expected to be a production process by 1986. 

a) Theoretical Material Balance for the Arc Furnace 
Process. 


The chemical reaction of this process is: 


Atomic Weights 


Masses, normal- 
ized to 1 kg Si 
out: 


, f 28t 

i { 2 . 


Si0 2 + 2C 


2(16) + 2(12) 
60 + 24 


Si + 2 CO 

28 + 2(12 + 16) 
28 + 56 


14 kg + 0.86 kg -*■ 1 kg + 2 kg 


According to Dow-Corning O) . r the industrially experienced 
and actual conversion efficiency of SiO^ to Si is 80%. 

The required input is thus: 2.68 kg Si0 2 for 1 kg Si out 

(Dow-Corning shows 2.71 kg SiOj) . 

It is assumed that half of the lost SiO^ input, or 0.2 7 kg, 
is used in the reaction: 

SiO 


Atomic Weights: 


2 

60 + 


C 

12 


Si o 

44 


+ 

+ 


CO 

28 


0.27kg + 0.05 kg * 0.2 kg + 0.12 kg 

where the SiO is lost at the top of the furnace. 

Thus, the total theoretical carbon input would be 0.91 kg C 
per 1 kg MG-Si output. 
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b) Theoretical Energy Balance. 

The energy required from the electric arc for the reduction of 

Si0 o to Si would firstly be expected to equal the heat of formation of 

o 

SiOg at the reduction temperature ( -AHf 2Qoo = 210kc:al , or 8.7 kWh for 
1 kg of Si formed. However, the carbon charged into the furnace also 
participates in the reaction, and, in its oxidation to CO, supplies 
approximately 25% of the required energy. Thus, the theoretical 
minimum energy required from the electric arc would be 6.6 kWh/kg Si. 

In addition, energy is required to heat the furnace charge to 
reaction temperature, an energy which is not recovered. Heating the 
Si0 2 charge (2.7 kg) to the reaction temperature will require approximately 
1.2 kWh, plus 0.1 kWh, for a total minimum theoretical energy requirement 
of 10.8 kWh/Si, outside of the heat loss of the furnace. 

c) Experienced Material /Energy Balances. 

Table II compares the experienced material and energy balances 
with the theoretical one. Experience data have been taken from ref. 

9 which lists "data from a major manufacturer of MG-Si", ref. 10 
which cites data from a not widely distributed Battelle report, and 
ref. 11 which gives data on an experimental arc furnace run at 
El kern. Several observations have resulted from the study of these 
data: 

i. The quality of the date does not justify better than 2-digit 
precision. 

ii. Mining and transportation energy expenditures for quartz 
and coal are quite variable, depending on source, method, 
location of user, etc., but are in all cases so small as 
to be negligible. (It is good that the referenced authors, 
as well as these investigators, have checked the data, so 
that this statement can be made with confidence.) 
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Balance Data 
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iii. Energy expenditures for plant operation (lighting, 
ventilation, crane operation, etc.) and energy con- 
tent of equipment have evidently not yet been analyzed. 

We added 1 kWh (el) per kg MG- Si produced as 
an estimate for these energy expenditures (Tables V-4A, 

4B, - 4C, each Line 1 (p.p. V-15, Y-l'., V-17). 

iv • The electrical energy consumption in the industrial 
processes seems only about a factor of 2 higher than 
theoretically required. 

v. The use of "fixed carbon" (F.C.) is also approximately 
a factor of 2 higher than theoretically required. 

i 

vi. The use of thermal energy is an order of magnitude larger 
than theoretical. The following reasons prevail: 

1. The oxidation of carbon goes only to CO, with only 
about 1/3 of the heat of combustion of C to CO^ 
utilized. 

2. The carbon sources used contain al«- combustible 
volatile components of high heat of combustion 
(hydrocarbons) which are not utilized. 

vii. Both ref.’s 9 and 10 add a considerable energy expenditure 

for the coking process. The petro coke is, however, a 
byproduct of gasoline refining, obtained by coking the 
heavy residues from the distillation process. This coking 
process provides a higher yield of gasoline. It seems, 
therefore, that the energy expenditure for the coking 
process should be attributed to the gasoline production and 
not to the "co’-.c". 
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viii. For the conversion from heat of combustion to electri- 
cal energy , the factor 0.33 has been used, corresponding 
the average efficiency of electric power generation in 
thermal plants. In this, the transmission/distribution 
losses have been omitted since large power users, 
such as arc furnace operators, are usually located 
close to generating plants. 

ix* The data given in ref. (9) and Table 13 of ref. (11) 
do not seem to agree with the standard specifications 
Q.2) for "low volatile coal" shown in ref. (9). 

x. The energy content of "wood chips" used in ref. (9) 
is approximately a factor of 2 lower than that of 
ref. (10) and (13), with credance given to the latter. 

The energy content of the volatile component of 
the woodchips on the unit mass basis, was originally 
assumed to equal that of metallurgical coal. However, 
comparing the total heat of combustion for wood thus 
obtained with that given in standard tables, after 
adjustment for moistrre content, showed that the 
energy content of the volatiles of wood, as the only 
unknown, had to be lower by approximately a factor 
of 2 than that of the volatiles of coal. There 
probably is better information on this subject 
available in the literature, but it was not deemed 
beneficial to the project to expend more effort on 
this subject, particularly since this energy does 
not participate in the reaction, but exits via the 
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off-gas. Thus, both input and output of the energy 
balance are reduced by an equal amount. 

xi . Ref. (9) lists the conversion efficiency of Si02 
to Si as 80%, but ret. (10) uses only 66%. Never- 
theless, ref. (10) arrives at approximately the same 
total energy consumption. 

xi i. . The July 1977 Dow-Corning quarter l'> report 

contains very useful detailed data on this subject, 

which have been augmented by additional information 

(15) 

obtained directly from Interlake. 

The final energy balance, which was capable of accounting 
for 90% of the input energy (Table V- 3) , contains the interesting 
finding that nearly half of the energy input to the process 
leaves the furnace as chemical energy in the off-gas. 

About 60% of this energy is contained in volatile gases, in 
good part originating from the wood chips. The other 40% is con- 
tained in the carbon monoxide which results from the fixed car- 
bon used to reduce the SiC^. This off-gas is not utilized at 
present, but is mostly burned off at the top of the furnace. A 
utilization of this off-gas could reduce the net energy 
consumption of the process to approximately one half its 
current value, and such an improvement has been assumed 
accomplished for the energy balance projections to 1982 and 
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Table V-3 
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For the purification of silicon by the standard 
SiHCl 3 process and the generation of the "Polylog", pub- 
lished data (16) and corroboration were used. The 

energy consumption data computed on this basis are shown in 
Table v-4A, lines 2, 3, and 4 for the SiHCl 3 generation, the 
SiHCl 3 distillation, and the SiHCl 3 reduction to purified 
silicon, respectively, and correspondingly in Table V-4B. 

As an alternate purification process, which can be ex- 
pected to have replaced the SiHCl 3 purification process by 1986, 
the projected energy consumption data for the SiH 4 process under 
development at Union Carbide Corp. ^ are detailed in lines 
2 and 3 of Table V-4C. 
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PROCESS ENERGY CONSUMPTION IN THE CURRENT (lt7«) MANUFACTURING SEQUENCE FOR SILICON 




i 




• 

1 

^ #d W ^ 

•• fi « « J 3 


*4 

!. 

r 

*3 

11 

a 

2 

Y_ 

« * * 

s a a $ s § 

a 

ft *♦ 

Id or 

«n r- 

#N I-! 

r 

* e* ^ 

J « r» «•» 2 

SR* 35* 

I 

•• 

d * —* 

£ 

1 

J 

F 

O IN 

*i o l in « 


3 

2« 

tii ii i 


! 

a 

<" • _ _ *» 
5 ® ' 1 8 S 



. . 9 III 


1 

2 

2 

I 

a 

'* 

ff 

1 S m ^ 

« ' o’ 8 • 


a - 

* 

* * 

* • ^ 

8 R . 5 


1 

1 

a 

a 

a 

& 

• * 

IN <N ^ JJ 2 


a 

*. “i o 2 * 

i s s a 5 s 


i 

fc 

« 

a 

a 

j 

i 

* 

>s ^3 « ^ ■ 

§ A u> »< & ll« 

“3*5 -T “ - 8 3 &S 

pi p -3 . v,u g> v. in p pa 

S V M P «J£ Jt *9 i" •* C «< 

a -4 Jt m4 H M a ». MW 

r- r- ui u w e* u w <■» u * “J 

>» <a o s « X •* * si ? 

• Pi • p • •«• p • •«• P • P te • >■ 

M jf —• M rt (fl h •; < O 0 

Conversion 
0.3kWh(el) 
per kUh(th) 

dto. 


1 ^ t 

. 1 3 3 I i 

8 m •* *» E S 

•3 w -H 0 d» « £ 

| 1 i 1 3 |I 

v a XT' n w *n — . *> «•> d 

fi ^ id 9 (AtddtdCfp 

£ f JX* .r«K t< 3 5 

o*N «S “S a 8 t °£- 5 5 

SS •<: - * - o)i N « H ^ 

Crt— W *• W>— Ul *j — O — 019 

fa 

•ft 

u • A 

6 d d™ 

W V •* V ff 

_ * ■ — * • 0-N* 

3 i.". S t £3 
53* u £e3 

V > f h< ■* 

w 3 S y w| 
v v f* m **- » 5 
i w w eo o* 

u u «3 


I * * • • • * j. • 

|h in m v m * j r% m 



o» in • 

• . « « 
•2 * i * w 

O e c * o • 
«n d >m O 
fa — x 


£ 2 


I 


5 «l 5 * M * lU 

•J u t m • 



TAILt V-4B 
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3.3 Energy Consumption in Sheet Generation 




The current process group for sheet generation contains 
Czochralski crystal pulling followed by wafer slicing 
either by ID diamond blade sawing, multiblade or multi- 
wire slurry sawing. 

The Czochralski process can be expected to be improved 
both with respect to cost and energy consumption. Competing 
are the heat exchanger method (HEM) of single crystal growth, 
and semicrystal casting. All of these processes, however, re- 
quire sawing, with very substantial kerf losses. Consequently, 
the energy balance of the sawing processes is dominated by 
the energy content of the silicon material lost in the kerf. 

The primary possibility for energy savings lies therefore 
in the replacement of the bulk crystal growing methods, with 
subsequent sawing, by one of the ribbon growing methods under 
development. In the projections, improved Czochralski 
crystal growing and sawing methods are used for 1982, and 
a ribbon growing method for 1986. 

a) Czochralski crystal growth and slicing. 

The analysis of the Czochralski crystal pulling process, 

(19) 

was based on a review of an earlier analysis and the 

( 20 ) 

addition of newer data from Texas Instruments and Dow 

(21) 

Corning . The projections contain primarily three 
improvements of energy impact: a reduction of the crucible 

usage, an increase in the furnace productivity, and additional 
technology advances which include better furnace design for 
reduced energy consumption. 

The projected crucible usage reduction is based on 
the assumption that crucibles can bo used for the equivalent 


of If individual crystal pulls, either with re-seeding or 
(quasi-) continuous pulling, rather than the currently 
practiced usage for one crystal only. 

The second projected reduction of energy consumption 
is connected with a projected furnace productivity 
increase. Approximately half of this productivity increase 
is expected to result from an increase in crystal diameter 
from the presently common, nominally 75 mm (3") diameter 
to nominally 100 mm (4") by 1982 and to 150 mm (6") by 1986. 

The other half of the productivity increase, however, is 
expected to come from a higher linear pull rate, thus more 
closely approximating the thermodynamically computed theoreti- 
cal maximum pull rate. This prediction of a linear pull 
rate increase is more risky as two currently not adequately 
explored effects are involved. The first concerns crystal 
perfection which may decrease with increasing pulling speed 
and may possibly prevent attainment of the expected pull rates. 
The second unknown is based on an analysis by Rea (22) 
who found that the radiative energy transfer from the melt 
surface and the heater environment to the grown crystal 

prevents any close approach to the limit growth rates com- 

(23) 

puted in the earlier thermodynamic analyses . This 
spurious radiative heat transfer could, in principle, be 
reduced by introduction of appropriate heat shields. To 
what degree this can be achieved in practice, without inter- 
fering with other aspects of the crystal growing process, 
needs to be explored. 

Some of the possible reduction in energy consumption 
is not just related to the productivity increase 
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through a reduction of the time for unit mass crystal growth, 
during which heat losses occur at a constant rate, but is 
directly connected with the dependence of the heat losses 
on the system geometry. The estimates are based on the - 
only conditionally valid - assumption that the heat losses 
in unit time are directly related to the crystal 

geometry change, ard thus would increase proportionately to 
the crystal radius. The increased mass pull rate, which 

results from the radius increase, however, leads to a reduc- 
tion of the energy losses inversely proportional to the radius. 
Finally, the required decrease of the linear pull rate results 
in a reduction of the energy losses inversely proportional 
to the square root of the crystal radius. 

Twice as large a reduction in energy consumption than 
by the geometry change is, however, expected to result from 
technology improvements, consisting in the use of better heat 
shielding and insulation in future furnace designs. Another 
technology advance has been assumed in the substantial decrease 
in usage of replacement parts for the crystal pulling furnace, 
which reflects itself in reduced energy content of the replace- 
ment parts used per unit mass of crystal pulled. 

Prior to slicing, the ingots are brought to constant 
diameter in a grinding process, which, by current production 
experience, costs 8% of the mass of the crystals grown. With 
the use of larger diameter crystals, projected for 1982, this 
grinding loss has been assumed to be reduced to 6%. 

A similar analysis has been carried out for the slicing 
process, both as currently practiced and with projected tech- 
nology improvements. Data on the current production processes 

/ 2 A ) 

include the multi-blade slurry sawing process and the in- 

side diameter diamond blade sawing process Also, ox - 
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perimental and projected data on advanced slicing methods were 
examined for the multi-blade slurry saw ' , for the 

ID diamond blade aaw^^ , and for the Yasunaga and Crystal 
Systems multiwire slurry sawing processes. OO ,31) The 
slurry and the blade packs used in the multi-blade process 
constitute a substantial indirect materials consumption with 
significant energy content, as do the diamond saw blades in 
the ID sawing process. 

The slicing technology improvements projected for the 
1982 production lines as far as they concern the energy 
balance, are the results of current experimental runs. For 
the multi-blade saw , the primary advancement will be a 25% 

blade thickness reduction, in combination with a 37.5% wafer 
thickness decrease, to 250 ym wafer thickness and 200 ym 

kerf, while maintaining the present practical wafer yield of 
95%. This results in a 50% increase in the mass to area 
conversion, to 0.9 m^ per kg silicon crystal. A similar 
reduction in wafer thickness is anticipated in the ID 
sawing process, but without reduction of the kerf. The 
slurry multi-wire saw, which also could be on the 1982 
production lines, can yield wafers with similar thickenss 
and kerf as the multi-blade slurry saw. 

For the 1986 projections, energy consumption data given 
for the EFG ribbon growing process but not yet 

reviewed by us, were used. 

In consequence of the discussed process improvements, 
the energy content of the wafers or ribbons of silicon is 
expected to fall from the 1978 value of 1537 kWh (el) /m 2 
to 880 kWh(el)/m 2 in 1982 and 165 kWh(el)/m 2 in 1986 (Tables 
V-4A to C, pp. V-15 to V-17) . . 
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3.4 Enerqy Consumption in the Solar Cell and Module 

Fabrication Process Sequenco 

The energy analysis of the solar cell fabrication pro- 
cess sequence is, in some respects, simpler than that for 
the sequence up through sheet generation, and more compli- 
cated in others. It is simpler, because one is dealing only 
with areas of silicon wafers or ribbons, no longer converting 
from one material form to another, or from mass units 
to area units. It is more complicated, however, because there 
exist many more process options and potential sequences. 

In any such sequence, the accumulation of yields from 
the individual process steps is as important for the energy 
consumption as it is for the costs. In fact, for the entire 
process from SiC >2 to finished cells, the energy content of 
the silicon lost in the various conversions and due to 
yields of the process steps far exceeds the energy going in 
direct line into the finished product. 

Consequently, the total energy expended for producing 
a unit of work-in-process (or finished product) leaving 
a given process step n is described by E n _i/A n + AE n ' 
where E n _^ is the total energy expended for producing a unit 
of work-in-process entering the respective step > AE n is 
the total energy needed to process a unit of work-in-process 
through the step, and > n is the yield of the process step. 

The total energy E T , N expended in a process sequence up to 
step N (inclusive) is then: 


N n 





where E Q is the energy content of a unit of work-in-process 
or direct material (e.g. wafer or ribbon) entering into step 
1 of the sequence. E is, for instance, the energy shown 

T $ 

in the right hand columns headed "Total Energy" in Tables V-4 

A to C and V-8A to C, or "Subtotal" and "Total" in Table V-6. 

These total energies cannot be summed, but can be used for 

entry as E^ ^ to another step or sub-sequence. Thus, the 

data in the left-hand columns headed "Direct Material Energy" 

or "Input Material" in Table V-6, and in Tables V-4A to C and 

V-8A to C, represent the values of 

The total enerqy "content" of the input work-in-process required 

E 

to produce a unit of good output from step n is n-1 , so that 

*n 

the energy content E of the material lost in the step n is 

L,n 

given by: 


E 


L,n 



E 

n-1 


(V-2) 


This quantity includes the energy content of all material lost 


in the preceding steps. The energy given by eq. (V-2) is an important 
ingredient of the "add-on-energy" E. for a unit of good work-in-process 

leaving process step n: 


E A,n 


E 


L,n 


+ AE ; 
n 


(V-3) 


Data generated by use of eq. (V-3) are included, e.g., in 
the right-hand column headed, "Total Add-On Energy" of 
Table V-5, while the left-hand column headed "Energy Content 
of Lost Silicon" includes data obtained by use of eq. (V-2) , 
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with 


E . = E _ + E. „ , 
n-1 n-2 A,n-1 


(V-4) 


Neither the add-on energies, nor the energy values for 
material lost in a process step, can be summed directly to 
obtain the total energy expenditures in a process sequence 
up to step N inclusively, because of the yields of the process 
steps subsequent to a given step n. Thus, for a unit area of 
product to leave step N, the total energy content E L N n 
of the material lost in step n is : 


E L,N,n 


E 


L,n 


N 


If 

k*n+l 




(V-5) 


and the total add-on energy for the sequence from step 1 to 
step N, inclusive is given by: 


N 



n=l 


N + 1 

IT » k 


(V-6) 


k=n+l 

with * 

Eq. (V-5) has been applied to obtain the "Totals" for the 
"Energy of Lost Silicon" in Table (V-7) , and an equivalent 
relationship for the other energy component totals. The 
farthest right-hand column of Table (V-7) represents data 
resulting from, application of eq. (V-6) . The cnerqy content 
of the total material lost is then simply: 
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(V-7) 


N 

E L,T ■ E A/T "j^ En • 


The quantities E L n and E A n , and correspondingly E L N n 
and E a t , include in the energy content that of all material 
lost in processing the good, remaining work-in-process. 

They are therefore not suitable for the determination of the 
total energy content of the material lost in a sequence of 
process steps, or of the "total net energy content" of 
the good work-in-process leaving such a sequence. This 
total net energy content does not include the energy content 
of the lost material. Thus its computation has to be based 
on the net energy contents of the input work-in-process, 
including consideration of the yields of the subsequent process 
steps. In analogy to eq. (V-2) , the net energy lost in 
step n for a unit of good work- in-process leaving this step is 


E L,n,n 



E n,n-1 


(V-8) 


where : 


J n, n -1 


= E, 


n-1 

k=l 


AE, 


(V— 9 ) 


is the net energy content of the input work-in-process. 
Similarly, the total net energy lost in step n for a unit 
of good work-in-process to leave a process sequence after 
step N is given by: 


E L,n,N,n 


E t 

L,n ,n 
N 

TT*. 

k=n+l 


(V-10) 
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and the total lost energy in all process steps by: 


N 


E L,n,T 



E L f n,n 
N+l 

I I *k 

k=n+l 


(V-ll) 


with X N+ i » 1. 


Corresponding to eq. (V-9) , the total net energy content of 
the good work-in-process, or finished product, leaving 
step N is: 

N 




k=l 


(V-12 ) 


This total is represented in the farthest right hand column labeled 
"Total" in Tables V-9A to C, while the entries in the columns labeled 
"Energy content in lost Si" of those tables represent 
data according to eq. (V-10) , and the totals of those columns 
correspond to eq. (V-ll) as well as eq. (V-?). 

In the solar cell processing sequence, we have accu- 
mulated the energy consumption data for the texture etching 
and the junction formation processes, and we have examined 
the direct material content of the encapsulation materials. 

The metallization and antireflection coating processes are 
presently being analyzed, and their energy data will be 
presented in a later report. In addition to these detailed 
analyses, we have accumulated the energy data for the 1982 
$2 . 90/W (pk) and the 1986 $G.!iO/W(pk) JPI. St rawman process 
sequences. Prom our detail analyses for t tie individual 
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processes and the available data for the remainder of the 
sequence, such as the Strawman data, we have synthesized 
overall process sequence energy consumption and payback time 
data. 

a. Energy Consumption of Junction Formation Processes 

The energy consumption of the various present and projected 
junction formation processes was studied in detail in 
connection with the analysis of their cost-effectiveness. 

In the evaluations, SAMICS overhead energy standards were 
used throughout, except for the 1974 and 1977 experience data 
shown in Table V-5, lines 1 and 2. 

The wagers or ribbons resulting from the Cz/slicing (1978 
and 1982), or EFG ribbon (1986) processes were tentatively 
assumed to be jubjected to the same texture etching process, with 
the same yield, independent of the junction formation processes j 

used or the time frame. This unified wafer preparation consists 
of a 3-step process sub-sequence , which includes the appli- 
cation of an etch-stop (wax) to the buck surface of the cell, tex- 
turing of the front surface by a hot NaOH etch, and removal of 
the etch-stop (Table V-6) . The wafers resulting from this 
process sub-sequence are assumed to form the input material for 
the various junction formation processes. 

The add-on energies for the junction formation processes 
are summarized in Table V-5. These processes have been grouped 
into present (Table V-5, lines 2 - 3) , near term (Table 
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V-5, line 4) and long range projected junction formation 
processes (Table V-5, lines 5-13). For comparison purposes, 
the 1974 data ^ ^ ^ for the direct energy consumed in a 
diffusion process have been included in Table V-5, line 1. 

Since the data of this reference include indirect energy 
consumption only for the entire plant, including wafer 
generation, cell processing, and environmental testing which 
is probably more connected with the space power cells pro- 
duced predominantly in 1974, the indirect energy consumption 
was, for the purposes of this report, allocated to the various 
process areas in proportion to their direct energy consumption. 

The 1978 Spectrolab phosphine diffusion process (Table V-5, 

line 3) , for which detailed experience data had been made 
(34) 

available , shows approximately an order of magnitude lower 
total equivalent electrical energy consumption than lies' 
numbers. This energy consumption is, in many respects, 
comparable to that for the POCI3 diffusion process contained 
in the JPL 1982 Strawman process (Table V-5, line 5) ^5) . It 
noteworthy that the "equipment energy," derived from the 
equipment price as outlined on p. V-4 of this report, is the 
item of highest energy consumption in the Spectrolab process 
(Table V-5, line 3) . 

Motorola's diffusion process data also » • >hosphine 

diffusion, arc used as a projection t cj l'lHG, with . .to much 
lower energy content, per m , of the silicon ribbon expected 
to be used then. This difference expresses itself in the 
energy consumption for lost silicon. Essentially counter- 
balancing this change, however, is Motorola's assumption ol 
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a lower yield. Another significant difference is in the 
material energy consumption, which is based on Spectrolab's 
substantial use of acids for quartzware cleaning, an item 
which is not included in Motorola's data. 

Also included in Table V-5, line 7 , is a 5-step process 
sub-sequence (Motorola) for the formation of both the pn junc- 
tion and the BSF layer. This 5-step process includes the 
phosphine diffusion step just discussed. Its five process 
steps and their energy consumption are detailed in Table V-7. 

Besides diffusion, ion implantation is presently used 
for junction formation. But in its present practice^^ , 
ion implantation is much more energy intensive than the diffu- 
sion process (Table V-5, line 4). The process energy of ion 
implantation is, however, projected to decrease significantly 
through improvements in equipment design. These designs 
stress higher projected throughput rates relative to material 
usage, and capital equipment and facility requirements. In 
contrast, for the diffusion process, the add-on energy is not 
expected to change substantia llv since significant future in- 
creases in furnace efficiency are not anticipated . 

It is also necessarv to consider associated processes 
in the sequence, and the yields in such seouences. Thus, 
the ion implantation process should never be considered 
without adding the needed annealing step (Table V-5, 
line 13) . Consequently, the 5-step pn junction and BSF 
layer diffusion process and the 2-sided ion implantation 
process with annealing (Table V-5, lines 2 and 9, respectively) 
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constitute comparable processes. Differing yields, or assump- 
tions of yields, can have significant influence on the energy 
consumption. Thus, the lower overall yield (96.4%) of the 
5-step Motorola sequence compared to that of the competing ion 
implantation process (98.4%) accounts for nearly 50% of the 
difference seen in the energy consumption of the 2 processes. 

This difference is in favor of the ion implantation for 
this particular example. Because of the importance of high 
yields for achieving processes with low energy consumption, 
most of the 1936 solar module fabrication sequence has been 
projected to be composed of individual processes of very high 
yields. Modifying the projected Motorola PH 3 diffusion 
process to be consistent with Spectrolab's experience of 99.9% 
yield and, in addition, a 17.5% "capacity factor" for the 
furnace energy consumption, would reduce its payback time 
from 24 to 9 days, and that of the 5-step sequence to 
33 days. This would then be slightly less than the energy 
payback time of the competing 2 -sided ion implantation 
process with annealing. 

The data of Table V-5 are illustrated in Figure V-l. 

It shows the projected ion implantation processes to require 
less add-on process energies than the projected diffusion 
processes. However, after adding the process energy for 
activation annealing to that of the implantation process leads 
to comparable energy expenditures and payback times for 
the future ion implantation processes and the diffusion processes. 
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Figure 1: 


Junction Formation 


It may also be noted that the energy payback times for the 
projected diffusion processes are more "firm" than those for 
ion implantation, since the latter are based on equipment which 
has not yet been built. In contrast, the projected diffusion 
processes represent relatively small extrapolations from current 
production equipment and practices, primarily modifying 
throughputs to meet the LSA-JPL output goals, 
b. Module assembly (encapsulation) 

The major energy contributions in the module assembly 
step seem to come from thj energy content of the encapsulation 
materials. Consequently, time has not yet been spent on determining 
the direct energy consumption in the encapsulation process 
(interconnect attachment, matrix connection, encapsulation 
material layup, potting material curing, junction box assembly) , 
and the facility energy consumption for the respective part 
of the plant. 

For the 1978 module assembly process, the encapsula- 
tion materials chosen were two 3 mm (1/8" thick) glass plates, 
along with a 1 mm thick layer of potting material. The cell 
packing density of 80% requires 1.2d m^ of encapsulation 
material for each square meter of cell area. Cell yield in 
this area was tentatively taken as 100%. With a glass energy 
content of 46.1 kWh/m^ ,and a potting material enerqy content of 
about 12 kWh/kg, both assumed to be predominantly thermal 
energy, the energy content of the encapsulation materials 
alone is about 130 kWh(th)/m of cell area. For the present 
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process, lies ( ^ gave no energy consumption data for the encapsu- 
lation sequence, but Solarex presented data (39) which are 
tentatively used in Table V-8A. 

For 1982, an 85% packing factor has been used for the 

solar cells in the module, reducing the encapsulation area 

to 1.18 m^/m^ of cells, and the energy content of the 

encapsulation materials to 120 kWh(th)/m^ of cells. At 

the same time, the energy consumption data from the JPL 

$2/W(pk) Strawman process have been used for the direct 
energy use and for the facility energy. (35) 

For 1986, the packing factor has been assumed to be 
further improved to 95%, and the back glass layer replaced with 
a 0.25 mm thick Mylar film, or another material of similar energy 
content. Use of the improved packing factor alone without back 
glass replacement, results in an encapsulation material energy content 
of 107 kWh (th) /m2 of cells , while use of the Mylar film backing 
leads to 62 kWh/m of cells. In addition, the relatively 
small direct and indirect energy consumption data of the 
1982 Strawman process for the encapsulation process group have 
also been used for 1986. 

c. Full solar module process sequences. 

Using the energy content of the input work- in-process 
to the solar module fabrication process sequence which was 
summarized in Tables V-4A to C for the current silicon wafers 
as well as projected wafers (1982) and ribbons (1'86), the 
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energy content of the completed modules was estimated or 
projected, based on available data. It turned out that the 
differences in these data for current processes are so 
large that they cannot be significantly improved by intro- 
ducing the data resulting from our texturizing and junction 
formation process analysis into these other sequences. 

A first analysis of the electrical power requirements 
of a solar cell production plant had been performed by Iles^. 
lies estimated the contemporary (1974) solar cell fabrication 
plant energy consumption from the installed equipment 
power ratings, and extrapolated the plant's productive 
capacity by assuming that a five times larger solar cell 
area could be handled annually in the same plant, 
with approximately the same energy consumption, by going 
from the then-prevalent 2x2 cm space power cells to 2" 
diameter cells. Deducting lies' direct energy consumption 
values for crystal growing and slicing leaves a direct 
electrical energy consumption of 120 kWh/m 2 of cell area pro- 
duced. Allocating the indirect energy consumption in pro- 
portion to the direct energy, leads to 125 kWh/m 2 of cell area 
for the solar cell process plant, These data contain, e.g. , 
a direct energy consumption for diffusion which is an order 
of magnitude larger than Spectrolab's experience data (see 

Table V- 5 , lines 1 and 2). On the other hand, lies' data 
do not contain the energy content of the indirect materials 
consumed, nor that of the equipment installed (see Table V-8A) . 
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The Solarex data ^ 39 ^ amount to about half the energy 
consumption found by lies, but include the energy contents 
of the indirect materials consumed and the equipment. Com- 
paring these data with those resulting from the analysis of 
texturizing and junction formation shows the Solarex energy 
consumption to be much higher than that given by Motorola for 
the texturizing process group# and by Spectrolab for diffusion. 

In contrast, the indirect energy consumption of Solarex (plant 
lighting etc.) amounts to only 1.7 kWh/m^ of cells processed 
and is much lower than any other data given for this item. 

We consequently increased this consumption to SOkWh/m^ of cells, 
or nearly 100% of the direct energy consumption, about in line 
with the other available data. One reason for this change 
is also the apparent omission by Solarex of the indirect 
energy consumption for the common areas, offices, etc. 

In the encapsulation area, we used the energy content of 
the direct materials from section V-3.4bof this report, amounting 
to 130 kWh/m^ of cells, for 2 sheets of glass and potting mater- 
ial. This number omits the energy content of the interconnectors, 
a junction box, or a frame possibly applied. Considering 
these facts, as well as a possibly lower packing factor for 
Solarex, Solarex' s value of 205 kWh/m of cells for this 
item appears quite compatible. For the direct and the indirect 
energy consumption in this process step, the Solarex data 
have been used as the only ones so far available. It msv 
also be noted that 100% cell yield has been assumed for this 
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encapsulation process group, which may be slightly optimistic. 

The consequence of these entries is a total energy 
consumption of 2179 kWh equivalent electrical energy consumed 
in the production of 1 m^ of encapsulated cell area. At 
80% cell yield from wafers to finished modules, it turns 
out that 1920 kWh, or 88% of the total module energy content, 
was already contained in the wafers entering the solar cell 
processing line. Of the 259 kWh added in the solar cell/ 
module process line, over one quarter is attributable to 
module assembly and encapsulation. 

For the 1982 projections, the summary numbers given for 
the JPL-Task IV Strawman process ( 35 ^ have been reviewed. 

These data include a 93% overall cell processing yield, and 
an 85% module packing factor. The energy data for the diffusion 
step in this Strawman process sequence agree quite well 
with those of the current Spectrolab diffusion process, 
except for the equipment and facility energy values, which 
are considerably lower in the Strawman process because of 
higher assumed throughput rates. This comparison on one 
significant process step gives a degree of credibility to the 
remainder of the data. Again, we used our energy content data 
for the double glass encapsulation. As Table V-8B shows, 
the energy content of the completed module has been reduced to 
about half of that of the 1978 module, but the energy content 
of the input wafers now constitutes 93% of the module energy 
content. Also, module assembly and encapsulation now consumes 
59% of the cell and module process energy. These shifts are 
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due primarily to the considerable reduction in energy con- 
sumption on the solar ceil process line, and are rein- 
forced by the yield improvement on this line from 80% 
to 90%. 

The projections to 1986 were similarly based on the JPL 
"Candidate Process" data^ 3 ' , which lead to a $0.50/W(pk) 
encapsulated module price. The input material is assumed 
to be silane purified, EFG grown ribbon, according to Table 
V-4C. Use of the "Candidate Process" leads to a total energy 
consumption of 206 kWh/m 2 of equivalent electrical energy 
(Table V-8C) , of which 172 kWh/m 2 , or 84%, is represented 
in the input ribbon material. Also, 57% of the solar cell 
and module process energy is added in module assembly, 
predominantly in the energy content of the encapsulation 
material . 

Tables V-9A to C present additional data to augment those 
of Tables V-8A to C. They give the mass fiow of silicon to 
the unit of finished encapsulated cell area, the net energy 
content per unit mass of the work-in-process at the key pro- 
cess stages, as well as the energy content of the silicon 
lost in the major process groups, and the net process energy 
of the material appearing in the good finished product, as 
contributed by these process groups. Again, these data 
are provided for the c r, omporary processes and for the 
projections to 1792 and 1986. 
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4. 


Conclusions 


The results of the energy consumption analysis are 
summarized pictorial form in Fig. V-2. This figure 
clearly demonstrates three points: 

a. most of the current high energy content is associated 
with the losses incurred in material conversion and 
in process yields. 

b. the biggest reductions in the energy consumption 
will be connected with the introduction of new 
processes for silicon purification and sheet 
generation. 

c. much of the reduction in cell processing energy comes 
from higher throughput rates; this effect is already 
observable now. 

The numbers on energy consumption are to be considered 
as rough approximations, since the data for the current 
process practice show a large spread, and since the future 
data represent projections. But in toto, the energy payback 
times can be expected to decrease rapidly from their recent 
value near 20 years to below 10 years by 1982 and to less 
than 2 years by 1986. This last prediction is somewhat obvious 
since the modules are expected to be close to cost-effective 
by then as replacement supplies for energy generally available 
from other sources. 
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Figure 2: Processing Energy Carponents For Module Manufacturing 
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Chapter VI 


Plans 

In progress are the analyses of metallization processes 
and of antireflection coatings, which will be reported in 
the next quarterly reports. The effort will then shift to 
masking and material removal processes, back surface 
field layer formation processes, and the evaluation of ingot 
versus ribbon technologies. 
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